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ABSIHACT 


Time-  and  depth-varying  cross-shore  currents  forced  by  gravity  and 
inffagravity  waves  in  the  nearshore  are  numerically  modeled.  Model  evaluation  is 
made  with  analytical  solutions  and  laboratory  data. 

Cross-shore  flow  is  simulated  with  two  linked  models:  a  depth- 
integrated  model  and  an  eigenfunction  solution  for  velocity  distribution  over  depth. 
Water  surface  gradients  generated  by  the  dq)th-integrated  model  are  used  to  drive 
the  profile  model.  Both  models  are  time  dependent.  Two  numerical  schemes  are 
used  for  the  dqpth-integrated  model,  Predictor-Corrector  and  Lax-Wendroff.  The 
Predictor-Corrector  scheme  exhibited  water  surface  and  velocity  oscillations  in 
space  when  executed  with  only  linear  terms.  Two  disjoint  solutions  in  space  were 
generated. 

The  depth-integrated  and  profile  models  were  applied  to  an  analytical, 
dq)th-integrated  surf  beat  solution  (Schaffer  and  Svendsen  1988  and  Schaffer  1990, 
1993).  The  solution  includes  both  forced  infragravity  waves  propagating  toward 
the  shoreline  and  free  infragravity  waves  released  in  the  breaking  process  and 
reflected  from  shore.  Mean  radiation  stress  forcing  was  added  to  the  model.  The 
d^th-integiated  model  results  compared  well  with  the  analytical  solution.  The 


X 


eigenfunction  profile  model  results  were  compared  to  the  periodic  analytical  profile 
model  of  Putrevu  and  Svendsen  (1995)  for  the  time-varying  portion  of  the  surf  beat 
solution.  The  results  were  identical. 

Laboratory  data  to  evaluate  the  depth-integrated  and  profile  models 
were  collected  at  the  SUPERTANK  Laboratory  Project.  Wave  and  current  data 
were  collected  in  a  large-scale  wave  channel  with  a  sand  beach.  The  models  did 
not  represent  mean  cross-shore  currents  well  due  to  limitations  in  the  radiation 
stress  calculation.  Time-varying  flows  were  represented  fairly  well. 

Characteristics  exhibited  by  the  profile  model  differed  outside  and 
inside  tire  surf  zone.  Outside  the  surf  zone,  most  vertical  variation  in  the  cross¬ 
shore  current  occurs  near  the  water  surface,  velocities  near  the  surface  and  bottom 
are  out  of  phase,  and  profile  variation  over  depth  is  fairly  symmetric.  Inside  the 
surf  zone,  the  mean  undertow  profile  is  the  dominant  feature,  time  variation  of  the 
velocity  is  fairly  constant  over  depth,  and  in  the  inner  surf  zone,  velocity  profiles 
are  uniform  over  depth. 
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Chapter  1 
INTRODUCTION 

As  gravity  waves  approach  a  beach,  they  steepen  and  break  in  water 
depths  approximately  equal  to  the  wave  height.  Breaking  waves  dissipate  energy 
and  induce  nearshore  currents  (longshore  and  cross-shore)  and  changes  in  the  water 
level.  The  surf  zone  is  the  region  from  the  shoreline  to  the  seaward  boundary  of 
wave  breaking.  This  zone  is  the  most  dynamic  coastal  region  with  sediment 
transport  and  bathymetry  change  driven  by  breaking  waves  and  wave-induced 
currents.  The  surf  zone  also  serves  as  the  interface  between  humans  and  the 
ocean,  so  surf  zone  processes  greatly  impact  the  human-ocean  inter-relationship. 

By  the  year  2(X)0,  45  percent  of  the  U.S.  population  will  live  in  coastal 
communities  (Senate  Coastal  Coalition  1995)  that  are  vulnerable  to  coastal  erosion, 
flooding,  and  pollution  problems.  U.S.  beaches  contribute  substantially  to  the 
U.S.  economy  with  tourism  producing  $50  billion  in  tax  revenue,  contributing 
$746  billion  to  the  Gross  Domestic  Product  in  1995,  and  supporting  more  than  14 
million  jobs.  Eighty-five  percent  of  all  tourism-generated  revenue  is  the  U.S. 
comes  from  coastal  states  (Houston  1995a,  b).  U.S.  beaches  are  the  key  element 
protecting  these  coastal  communities  and  generating  tourism  dollars.  An 


understanding  of  nearshore  coastal  processes  is  the  basis  for  engineers  and 
scientists  to  make  wise  decisions  when  endeavoring  to  "solve"  problems  of 
shoreline  and  bathymetry  change  (erosion/accretion),  coastal  flooding,  and 
pollutant  transport.  Engineering  tools  used  in  decision  making  take  the  form  of 
empirical  relationships,  physical  models,  and  numerical  models.  In  this  thesis,  all 
three  types  of  tools  are  exploited,  although  the  major  emphasis  is  numerical 
modeling.  Numerical  modeling  is  increasingly  the  tool  of  choice  due  to  its 
flexibility  and  efficiency. 

Although  the  most  dramatic  nearshore  feature  observed  by  someone  on 
the  beach  is  the  breaking  gravity  waves,  the  related  processes  of  nearshore 
circulation  and  infragravity  waves  exhibit  higher  correlation  to  sediment  and 
pollutant  transport.  Infragravity  waves,  oscillations  with  periods  of  tens  to 
hundreds  of  seconds,  increase  in  amplitude  near  the  shoreline,  in  the  region  where 
gravity  waves  (short  waves)  have  been  dissipated.  The  energy  in  infragravity 
frequencies  can  be  substantial,  even  exceeding  energy  in  short-wave  frequencies 
(Wright  et  al.  1982,  Guza  and  Thornton  1982). 

Nearshore  currents  are  composed  of  motions  at  many  scales,  forced  by 
several  processes.  Currents  are  composed  of  wave-driven  currents,  tidal  currents, 
wind-driven  currents,  and  gravity  and  infragravity  wave  oscillations.  This  study 
focuses  on  the  wave-driven  and  infragravity  wave  components  of  nearshore 
currents.  The  time-  and  depth-dependent  cross-shore  flow  is  modeled  under  the 
forcing  of  gravity  and  infragravity  waves. 


The  following  section  gives  background  information  on  nearshore 
circulation  and  infragravity  wave  modeling.  Section  1.2  provides  an  overview  of 
the  material  covered  in  the  following  chapters. 

1.1  Background 

The  study  of  nearshore  circulation  was  revolutionized  by  the 
introduction  of  the  concept  of  radiation  stress.  Radiation  stress  was  applied  to 
calculate  mean  water  elevations  and  currents,  undertow,  and  wave-group  forced 
infragravity  waves. 

In  the  early  1960s,  Longuet-Higgins  and  Stewart  (1962,  1963,  1964) 
introduced  the  concept  of  radiation  stress,  a  stress  tensor  associated  with  gravity 
waves.  Gradients  in  radiation  stress,  due  to  non-uniform  wave  amplitudes,  induce 
currents  and  changes  in  water  level.  Radiation  stress  forcing  has  been  applied  to 
problems  of  wave  setup  and  setdown,  longshore  currents,  undertow,  and  wave- 
group  forced  inffagravity  waves.  Radiation  stress  is  an  additional  term  in  the 
momentum  equation  for  the  solution  of  water  surface  elevation  and  velocity. 

Nearshore  circulation  is  generally  modeled  with  the  depth-integrated 
equations  of  mass  and  momentum.  The  equations  are  short-wave  averaged,  and 
formulated  much  like  tidal  circulation  models  with  additional  short-wave  forcing. 
Radiation  stress  is  calculated  externally  and  then  imposed  as  an  external  force. 
Generally  the  models  are  formulated  to  be  time-dependent,  but  in  practice  are  used 
with  steady-state  forcing.  The  steady-state  limitation  stems  from  the  use  of  steady- 
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state  short-wave  models  that  provide  the  radiation  stress  forcing.  Depth-integrated 
circulation  models  provide  water-levels  and  depth-integrated  currents  over  a  one- 
or  two-dimensional-horizontal  domain. 

Nearshore  circulation  models  are  used  to  study  longshore  current 
distributions,  wave  setup,  two-dimensional  current  patterns  (including  rip  currents 
and  cell  circulation),  and  ultimately  predict  sediment  transport  and  bathymetry 
change.  Examples  of  nearshore  circulation  models  include  steady  longshore 
current  models,  e.g.,  Bowen  (1969),  Longuet-Higgins  (1970a,  b),  Thornton 
(1970),  Thornton  and  Guza  (1986),  and  Larson  and  Kraus  (1991);  and  steady  and 
time-dependent  two-dimensional-horizontal  circulation  models,  e.g,,  Noda  (1974), 
Birkemeier  and  Dalrymple  (1975),  Keely  and  Bowen  (1977),  Ebersole  and 
Dalrymple  (1979),  Vemulakonda  (1984),  Wind  and  Vreugdenhill  (1986),  and 
Cialone  et  al.  (1994).  Examples  of  linked  wave,  current,  and  sediment  transport 
models  are  discussed  by  de  Vriend  et  al.  (1993). 

Longshore  currents  are  generally  relatively  constant  over  depth.  Wave 
setup,  the  vertical  nonuniformity  of  radiation  stress,  and  mass  transport  combine  to 
generate  a  depth-varying  cross-shore  flow  or  undertow.  The  mean  undertow 
profile  is  parabolic  with  a  maximum  offshore  flow  near  the  bottom  and  an  onshore 
flow  above  trough  level.  Undertow  may  be  relatively  strong,  generally  eight  to  ten 
percent  ofVJh  near  the  bottom,  where  g  is  gravitation  acceleration  and  h  is  the 
water  dqpth.  The  vertical  profile  of  the  undertow  is  determined  as  a  balance 
between  radiation  stress,  the  pressure  gradient  from  the  sloping  mean  water 


surface,  and  vertical  mixing.  The  first  quantitative  analysis  of  undertow  was  given 
by  Dyhr-Nielsen  and  Sorensen  (1970).  Solutions  for  the  undertow  profile  are 
given  by  Dally  and  Dean  (1984),  Hansen  and  Svendsen  (1984),  Stive  and  Wind 
(1986),  and  Svendsen  et  al.  (1987). 

The  quasi-three-dimensional  modeling  approach  has  been  used 
successfully  to  extract  steady  undertow  and  longshore  current  profiles  based  on 
local  wave  and  water  level  forcing  (de  Vriend  and  Stive  1987,  Svendsen  and 
Lorenz  1989,  Svendsen  and  Putrevu  1990,  Sanchez-Arcilla  et  al.  1992,  and 
Svendsen  and  Putrevu  1994).  The  quasi-3-D  approach  provides  information  on 
current  variation  through  the  vertical  without  the  complication  of  a  full  3-D  model. 
The  modeling  approach  has  more  recently  been  applied  to  situations  with  time- 
varying  forcing.  Van  Dongeren  et  al.  (1994)  calculated  the  development  of 
undertow  and  longshore  current  profiles  from  zero  to  steady-state  forcing  using  a 
quasi-steady  approach.  Steady-state  profiles  are  solved  analytically  using  the 
method  of  Svendsen  and  Putrevu  (1994).  A  quasi-3-D  approach  has  been  used  to 
extract  time-varying  profiles  in  large-scale  circulation  models  (e.g.,  Davies  1987, 
1980).  Putrevu  and  Svendsen  (1995)  develop  an  analytic  solution  for  time-varying 
inftagravity  wave  velocity  profiles  in  the  surf  zone,  assuming  constant  periodicity 
in  all  forcing.  Time-varying  profiles  have  not  been  compared  to  laboratory  or  field 
measurements. 

Infragravity  waves  fall  into  three  categories:  (a)  forced  long  waves, 

(b)  edge  waves,  and  (c)  leaky  waves.  Forced  long  waves  are  generated  by 


giEdients  in  mdintion  stress  found  in  wave  groups  causing  a  lowering  of  the  mean 
water-level  under  high  waves  and  a  raising  under  low  waves  (Longuet-Higgins  and 
Stewart  1962).  The  forced  wave  travels  at  the  group  speed  of  the  short  waves, 
hence  is  bound  to  the  wave  group.  Edge  waves  are  freely  propagating  long  waves 
which  reflect  from  the  shoreline  and  are  trapped  along  the  shore  by  refraction. 
Long  waves  may  be  progressive  or  standing  along  the  shore.  Edge  waves  travel 
alongshore  with  an  antinode  at  the  shoreline,  and  the  amplitude  decays 
ejqwnentially  offshore.  Leaky  waves  are  also  freely  propagating  long  waves  or 
standing  waves,  but  they  reflect  from  the  shoreline  to  deep  water  and  are  not 
trapped  by  the  bathymetry.  Proposed  generation  mechanisms  for  propagating  long 
waves  include  time-varying  break  point  of  groupy  waves  (Symonds  et  al.  1982), 
release  of  forced  waves  through  wave  breaking  (Longuet-Higgins  and  Stewart 
1964),  and  nonlinear  wave-wave  interactions  (Gallagher  1971). 

Field  studies  have  clearly  identified  forced  long  waves  and  edge  waves 
in  the  nearshore  (see  Huntley  et  al.  1981  and  Oltman-Shay  and  Hathaway  1989). 
The  relative  amount  of  infragravity  energy  and  incident  gravity  wave  energy  is 
often  considered  a  function  of  the  surf  similarity  parameter  (Holman  and  Sallenger 
1985,  Holman  1986),  with  infragravity  energy  dominating  for  low  values  of  the 
surf  similarity  parameter  (<  1.5).  For  low  surf  similarity  values,  the  energy 
q)ectrum  at  incident  frequencies  is  generally  saturated  (the  spectral  energy  density 
is  independent  of  the  offshore  wave  height,  due  to  wave  brealdng),  but  at 
infragravity  frequencies,  the  energy  density  increases  linearly  with  increasing 
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offshore  wave  height  (Guza  and  Thornton  1982,  Mase  1988).  This  implies  that 
shortrwave  forcing  is  a  significant  source  of  infragravity  wave  forcing  (Svendsen 
and  Putrevu  1996).  Storm  conditions  with  high  steepness  waves  tend  to  have  low¬ 
valued  surf  similarity  parameters,  so  infragravity  waves  are  prevalent  in  storms. 
Svendsen  and  Putrevu  (1996)  suggest  that  a  substantial  amount  of  infragravity  wave 
generation  takes  place  in  the  surf  zone,  based  on  recent  measurements  of  free 
infragravity  waves  in  the  offshore  (e.g.,  Elgar  et  al.  1992  and  Berbers  et  al.  1992) 
and  proposed  surf  zone  generation  mechanisms. 

1.2  Thesis  Content 

The  purpose  of  this  study  is  to  numerically  model  the  time-  and  depth- 
varying  cross-shore  current  forced  by  gravity  waves  and  infragravity  waves  in  the 
nearshore.  The  model  results  are  compared  to  analytic  solutions  and  laboratory 
measurements. 

First,  the  governing  depth-integrated  equations  of  mass  and  momentum 
and  the  depth-dependent  momentum  equations  are  derived  in  Chapter  2.  These 
equations  serve  as  the  basis  for  model  development  in  the  following  chapters. 
Chapter  3  describes  the  formulation  of  a  depth-integrated  model,  including  gravity 
and  infragravity  wave  forcing.  The  finite-difference  model  calculates  water  surface 
elevations  and  depth-integrated  currents.  Chapter  4  describes  the  formulation  of  an 
eigenfunction  solution  for  calculation  of  the  local  velocity  profile.  The  profile 
model  is  forced  with  output  from  the  depth-integrated  model. 


The  depth-integrated  and  profile  models  are  applied  using  forcing  from 
the  analytical  surf  beat  model  of  Schaffer  and  Svendsen  (1988)  and  Schaffer  (1990, 
1993).  The  analytical  surf  beat  model  provides  the  offshore  long- wave  input  and 
the  time-varying  short-wave  forcing  required  to  execute  the  depth-integrated  and 
profile  models.  The  depth-integrated  model  is  compared  to  the  analytic  surf  beat 
solution  and  errors  are  evaluated.  The  profile  model,  with  the  surf  beat  forcing,  is 
compared  to  the  periodic  profile  model  of  Putrevu  and  Svendsen  (1995). 

Chapter  6  describes  data  collected  during  the  SUPERTANK  Laboratory 
Data  Collection  Project.  SUPERTANK  was  a  large-scale,  nearshore  hydrodynamic 
and  beach  profile  change  study  conducted  at  Oregon  State  University.  The  chapter 
describes  how  the  data  were  collected  and  analyzed.  Sample  data  are  also 
presented.  The  SUPERTANK  data  are  applied  in  Chapter  7  to  drive  the  depth- 
integrated  and  profile  models  and  to  evaluate  model  results.  Last,  Chapter  8 
summarizes  the  study  and  provides  conclusions,  limitations,  and  suggestions  for 
further  study. 


Chapter  2 

DERIVATION  OF  GOVERNING  EQUATIONS 

The  purpose  of  this  chapter  is  to  derive  and  discuss  the  equations  that 
govern  the  nearshore  circulation.  The  depth-integrated  continuity  and  momentum 
equations  and  equation  for  variation  of  flow  over  depth  are  derived.  These 
equations  will  be  used  in  Chapters  3  and  4  to  develop  models  for  depth-integrated 
and  depth-varying  flows,  respectively. 

2.1  Depth-Integrated  Equations  of  Continuity  and  Momentum 

The  equations  governing  the  nearshore  flow  are  the  wave-averaged 
continuity  and  momentum  equations. 

2.1.1  Continuity  equation 

The  continuity  equation,  assuming  incompressible  flow,  is  given  by 

^  .0  (2.1) 

dx  dy  dz 

where  u  is  the  cross-shore  velocity,  v  is  the  longshore  velocity,  w  is  the  vertical 
velocity,  p  is  the  water  density,  and  t  is  time.  The  right-handed  coordinate  system 
is  defined  with  x  in  the  on-shore  direction,  y  in  the  longshore  direction,  and  z 


vertical,  measured  from  the  still-water  surface  (see  Figure  2.1).  Integrating 
Equation  2. 1  from  the  bottom,  z  =  -ho,  to  the  water  surface,  z  =  f ,  and  using 
Leibnitz’s  Theorem  gives 


±]udz^±].dz-u,^-.,K 

dh,  dh, 
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cbc  ^  dy  ^  ■*« 


(2.2) 


The  subscripts  -ho  and  f  denote  that  the  variable  is  evaluated  at  these  z  elevations. 


Next,  the  kinematic  boundary  conditions  (no  flow  through  the  interface) 
are  applied.  The  kinematic  free  surface  boundary  condition  is 


The  bottom  boundary  condition,  assuming  the  bottom  position  does  not  change 
with  time,  is 
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where  u„,  v,^,  and  w„  are  short-wave  velocities,  and  U,  V,  and  W  are  short-wave- 
averaged  velocities.  Applying  the  boundary  conditions  (Equations  2.3  and  2.4), 
Equation  2.2  becomes 
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where  u\  v’,  and  w’  denote  turbulent  fluctuations. 

Velocity  is  divided  into  three  components:  a)  short-wave-averaged 
component  (17,  V,  0),  b)  short-wave  component  v„,  wj,  and  c)  turbulent 
fluctuation  («’,  v’,  w”).  Substituting  these  components  into  Equation  2.5  and  time 
averaging  (over  the  short-wave  period)  gives 
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where  is  the  short-wave  trough  elevation.  The  overline  denotes  short-wave- 
averaged  components.  The  time-averaged  short-wave  velocities  are  zero  below 
trough  level 


“w  “  ''h.  “  ®  below  trough  level 


(2.7) 


For  the  region  above  trough  level,  the  integral  of  the  short-wave  velocities  u„  and 
V,,  are  x  and  y  components  of  wave-induced  mass  transport,  which  will  be  denoted 
as  Qya  and  Q^,  respectively.  The  first  two  terms  in  Equation  2.6  are  the  volume 
fluxes  below  mean  water  level  and  are  defined  as 
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Applying  these  definitions  reduces  Equation  2.6  to 

(2.9) 

dx  By  dt 


2.1.2  Momentum  equations 

The  momentum  equations,  neglecting  viscous  stresses,  are  given  by 
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where  P  is  pressure  and  g  is  gravitational  acceleration.  Integrating  Equations  2. 10- 
12  over  depth  and  using  Leibnitz’s  Theorem  gives 
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Applying  the  kinematic  boundary  conditions  (Equations  2.3  and  2.4),  assuming  that 
bed  level  is  constant  in  time,  simplifies  Equations  2.13-15  to 
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(2.18) 


dt 


C  c  C  p 

j wdz  I vwdz  +  g(C  +  hf)  -  — — 


,  ,  dhn  ,  dhn 

-  w  ,(w'  +  u' — -  +  1=0 

dx  6y 


The  last  three  terms  on  the  left-hand  sides  of  Equation  2. 16-18  represent  turbulent 
stresses  at  the  top  of  the  bottom  boundary  layer.  The  bottom  boundary  condition 
(Equation  2.4)  applies  only  to  the  wave-averaged  and  short-wave  components  of 
velocity,  and  thus  the  turbulent  stresses  do  not  cancel  out  when  the  boundary 
condition  is  applied.  The  notation  r^p,  Tj/p,  and  Tj^p  will  be  used  to  represent 
these  bottom  stresses,  respectively. 

Time  averaging  Equations  2.16-18  results  in 
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The  time  average  of  the  vertical  velocity  is  zero,  so  the  first  term  in  Equation  2.21 


is  zero.  Equation  2,21  can  be  rewritten  as 
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The  second  and  third  terms  on  the  right-hand  side  of  Equation  2.22  are  zero  by 

linear  theory  due  to  the  time  integration.  These  terms  will  be  neglected.  The  last 

term  on  the  right-hand  side  of  Equation  2.22  is  much  smaller  than  the  horizontal 

bottom  stresses,  and  in  Equations  2.19  and  2.20  and  thus  will  be  neglected. 

Substituting  Equation  2.22  into  2. 19  yields, 
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and  Equation  2.23  becomes 
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Similarly,  Equation  2.20  becomes 
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The  wave-induced  radiation  stress  components  are  defined  as 
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(2.27) 
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Tnsprring  the  time-averaged,  short-wave,  and  turbulent  velocity  components  (U,  u„, 
u’,  V,  V,,,  v’,  w„,  w')  into  Equations  2.25  and  2.26,  using  the  definitions  given  by 
Equations  2.27-2.29,  gives 
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The  turbulent  shear  stresses  (Reynolds  stresses)  will  be  parameterized  by 
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where  v  is  the  eddy  viscosity,  and  U  and  V  are  depth-integrated  short-wave- 
averaged  velocities.  The  turbulent  normal  stresses  (Equations  2.32  and  2.34)  are 
typically  small  and  will  be  neglected.  The  x  and  y  momentum  equations  become 
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The  bottom  stress  terms  will  be  parameterized  by 


(2-37) 

where  ^  is  the  bottom  friction  coefficient,  and  and  are  the  short-wave 
velocity  amplitudes.  This  formulation  is  a  linearized  version  of  the  standard 
quadratic  bottom  friction  (e.g.,  Jonsson  1966)  which  assumes  a  weak  current,  i.e., 
u„  >>  U  and  >  >  V. 


2.2  Depth-Varying  jc-Momentum  Equation 

Returning  to  Equation  2.10,  prior  to  depth  integrating  the  x-momentum 
equation,  and  substituting  the  three  velocity  components  (short-wave  averaged, 
short  wave,  and  turbulent),  gives 

d(U+uJ  ^  ^  e[(t/+«J(V+v,,)]  ^ 

dt  *  dx  By  Bz  (2.39) 

_  _  1  ^  _  B(u')^  _  Bu'v'  _  Bu'w' 
p  Bx  Bx  By  Bz 

Integrating  the  z  momentum  equation  (Equation  2.12)  to  an  arbitrary  elevation  z, 
time  averaging,  and  applying  the  kinematic  free  surface  boundary  condition  gives 


—  =  5(C  -  z)  -  + 

P 


d 

-J'uwdz  +  — jvwdz 


(2.40) 


The  last  two  terms  on  the  right-hand  side  of  Equation  2.40  are  second  order  and 
are  neglected,  as  discussed  previously.  Time  averaging  Equation  2.39  and 
substituting  the  mean  pressure  from  Equation  2.40  into  2.39  yields 


# 


du  ,  du'^  .  dUV  ^  ^  du^w^ 

-  +  -  +  -  4*  -  +  -  +  - 

dt  etc  dy  _  dy  dz  (2.41) 

d(u'^  -  w'^)  _  du*\'  _  du'w' 

^  dx  dx  dy  dz 


Since  the  wave  motion  is  periodic,  the  short-wave  average  of  the  short-wave 
components  equal  zero  (Equation  2.7). 

We  are  considering  only  the  cross-shore  flow,  so  we  assume  the  bottom 
contours  are  straight  and  parallel  in  the  longshore  direction.  This  means  that  the 
terms  d/dy  =  0,  and  Equation  2.41  reduces  to 


W  ^  ^  ^  ^  du^w^ 

dt  dx  dx  dz 

6C  d(u'^  -  w'^)  _  du'w' 
^  dx  dx  dz 


(2.42) 


Stive  and  Wind  (1982)  showed  that  the  turbulent  apparent  normal  stress  term, 

(i/’2  -  is  approximately  5%  of  the  wave-induced  apparent  normal  Stress  term 
in  surf  zone.  Thus,  they  concluded  that  the  effects  of  breaking-induced  turbulence 
are  weak.  This  term  will  be  neglected.  The  term  represents  the  horizontal 
shear  stress  created  by  the  oscillatory  wave  component  and  is  also  typically 
neglected  in  nearshore  current  modeling.  See  Putrevu  and  Svendsen  (1993)  for 
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discussion.  The  term  is  identically  zero  for  horizontal-bottom  wave  theories  except 
in  the  bottom  boundary  layer  (Battjes  et  al.  1990). 

Next,  we  introduce  a  Reynolds  stress  closure  for  the  turbulent  shear 
stress  « V’,  where  the  turbulent  shear  stress  is  related  directly  to  the  mean  flow 
through  the  eddy  viscosity,  p„ 

u'w'  =  -V,  —  (2.43) 


The  eddy  viscosity  is  assumed  here  to  be  constant  over  depth.  Applying  this 
closure  and  neglecting  the  terms  (u’^  and  reduces  Equation  2.42  to 

m  -  ±L  -  K  -  ~  O  _  dUi  (2.44) 

dt  8!z\  ^  dz)  *  etc  etc  etc 

The  first  term  in  Equation  2.44  represents  the  time  variation  in  the  cross-shore 
flow.  Typically  undertow  has  only  been  analyzed  for  steady  situations  (dU/dt 
=  0).  The  purpose  of  this  study  is  to  extend  analysis  to  unsteady  flow,  where 
dU/dt  is  important.  The  second  term  in  Equation  2.44  represents  the  "stiffiiess"  of 
the  water  column.  The  value  of  the  eddy  viscosity  will  determine  to  what  extent 
surface  stresses  are  transferred  through  the  water  column.  The  terms  on  the  right- 
hand  side  of  Equation  2.44  will  be  treated  as  the  driving  forces  of  the  cross-shore 
flow.  The  first  term  on  the  right-hand  side  is  the  dominant  forcing  term.  This  is 
the  setup  term,  and  it  has  no  vertical  structure  (constant  over  depth).  The  second 
term  on  the  right-hand  side  is  the  wave  radiation  stress  term.  This  term  has  a 
significant  vertical  variation  with  a  large  concentration  above  the  wave  trough  level 


and  a  fairly  uniform  vertical  distribution  below  trough  level.  This  variation  over 
the  vertical,  to  a  large  part,  drives  the  undertow.  The  last  term  in  Equation  2.44  is 
the  momentum  flux  from  the  cross-shore  current,  and  it  is  typically  neglected  in 
cross-shore  current  modeling.  Svendsen  et  al.  (1987)  showed  that  this  term  can  be 
significant,  although  it  is  smaller  in  magnitude  than  the  setup  and  radiation  stress 
terms.  This  term  would  typically  be  included  on  the  left  side  of  the  equation  and 
not  as  a  driving  term  on  the  right-hand  side.  For  simplicity,  we  will  be  treating  it 
as  a  known  quantity  (from  the  depth-integrated  model)  that  is  averaged  over  depth. 


Chapter  3 

DEPTH-INTEGRATED  MODEL  FOR  INFRAGRAVITY  WAVES 

This  chapter  describes  the  first  of  two  model  components  used  to 
calculate  the  horizontal  and  vertical  variations  in  the  time-varying  infragravity 
wave  velocity  field.  The  approach  used  is  quasi-three-dimensional,  where  the 
depth-integrated  velocity  is  solved  for  the  domain  first,  and  then,  as  described  in 
Chapter  4,  the  horizontal  velocity  components  at  several  elevations  in  the  vertical 
are  solved  locally.  The  local  solution  over  the  vertical  is  more  efficient 
computationally  than  a  full  three-dimensional  model.  Quasi-three-dimensional 
models  provide  information  on  velocity  variation  over  the  vertical  with  only 
marginally  more  computational  effort  than  a  depth-integrated  model.  For  the 
depth-integrated  model,  this  chapter  provides  discussion  of  the  model  assumptions, 
solution  techniques,  input  and  output  parameters,  and  an  example  application. 

The  depth-integrated  model  is  based  on  the  equations  of  continuity  and 
momentum  given  in  Equations  2.9,  2.35,  and  2.36.  For  the  present  application, 
additional  simplifications  have  been  made  to  the  equations.  The  first  assumption  is 
longshore  homogeneity  (3/3y  =  0).  This  assumption  is  made  because  the  model 
will  be  applied  to  data  taken  in  a  long,  narrow  wave  channel.  The  second 


assumption  is  that  the  cross-shore  and  longshore  components  of  momentum  flux 
(dlfi/dx  and  dV^/by)  are  small.  In  the  SUPERTANK  laboratory  measurements, 
blP/dx  is  small.  The  effect  of  one  current  component  on  the  other  {UV  terms)  will 
also  be  neglected,  although  Svendsen  and  Putrevu  (1994)  foimd  that  the  nonlinear 
interaction  has  a  dispersion  effect  on  the  longshore  current.  The  third  assumption 
is  to  neglect  the  wave-current  interaction  which  in  the  momentum  equations  occur 
in  the  trough-to-crest  region.  The  U  and  V  velocities  are  undefined  in  this  region, 
and  presently,  no  guidance  exists  on  specifying  these  values.  Under  these 
assumptions.  Equations  2.9,  2.35,  and  2.36,  respectively,  become: 
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For  the  applications  in  this  thesis,  waves  are  normally  incident  and  =  0.  Thus, 
the  longshore  flux  (j2y)  and  velocity  (V)  are  also  zero. 


3.1  Numerical  Solution 

Two  numerical  schemes  have  been  used:  a  two-dimensional  (x-y) 
Predictor-Corrector  method  and  a  one-dimensional  (x)  Lax-Wendroff  scheme. 
Both  methods  are  second-order  accurate,  and  the  solutions  are  explicit. 


3.1.1  Predictor-Corrector  model 

The  Predictor-Corrector  model  SHORECIRC  (Van  Dongeren  et  al. 
1994)  was  used  as  the  baseline  model.  An  absorbing-generating  offshore  boundary 
condition  was  implemented  for  arbitrary,  time-varying  surface  elevation. 


The  model  uses  the  following  definitions: 
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The  spatial  derivatives  are  calculated  with  second-order,  central  finite  differences. 
At  the  offshore  grid  boundary,  second-order  forward  differences  are  applied,  and 
at  the  shore  boundary,  second-order  backward  differences  are  applied. 

The  predictor  step  is: 


(3.7) 
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Q  *  =  <2 »  +  —13F".  -  F":‘) 


%  ‘  %  *  ') 


where  i,  j  subscripts  represent  the  x,  y  grid  indices  where  the  variable  is  evaluated, 
the  superscripts  represent  the  time  step  (e.g.,  n  and  n-1)  when  the  variable  is 
evaluated,  the  *  superscript  represents  an  intermediate  time  step,  and  At  is  the  time 
step.  The  corrector  stq)  is: 


C  =  -  C‘) 


(3.10) 


=  <?/  +  —(5Fu  +  8F".  -  F":^) 


(3.11) 


n  =  O  "  +  — (5G.  *  +  8G.."  - 

^yu  ^Vij  12'  ' 


(3.12) 


The  boundary  conditions  applied  in  the  model  include  a  no-flux  condition  at  the 
beach  (no  flow  into  the  beach),  an  absorbing-generating  offshore  boundary 
condition,  and  periodicity  at  the  lateral  boundaries. 

Boundary  conditions 

Initially,  the  offshore  boundary  was  specified  as  the  measured  total 
infiragravity  wave  surface  elevation  f .  Although  this  condition  is  theoretically 
correct,  practically,  it  requires  that  the  reflected  wave  generated  in  the  model 
exactly  matches  the  measured  reflected  wave.  If  the  calculated  and  measured 
reflected  waves  differed,  the  difference  was  re-reflected  into  the  model  from  the 
boundary.  Also,  this  mode  of  operation  requires  that  the  "solution,"  the  total 
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surface  elevation  at  the  boundary,  be  known  a  priori.  Therefore,  an  absorbing¬ 
generating  offshore  boundary  condition  was  developed. 

The  offshore  boundary  generates  a  specified  incoming  wave  and 
absorbs  the  outgoing  reflected  wave  using  the  method  of  characteristics  (Kobayashi 
et  al.  1987,  Van  Dongeren  et  al.  1994).  Assuming  shallow-water  conditions  with 

c  =  (3.13) 

where  ft  =  /Jq  +  f ,  the  continuity  equation  (Equation  3.1)  can  be  rewritten  as 

+  U—  =  0  (3.14) 

dt  dx  dx 

and  the  x-momentum  equation  (Equation  3.2)  can  be  rewritten  as 

.  ci?£  =  ^  (315) 

dt  dx  dx  p  j  p 

Note  that  the  longshore  momentum  flux  term  has  been  re-introduced. 

Adding  and  subtracting  Equations  3.14  and  3.15  generates  the 
characteristic  form  of  the  equations.  Using  the  definitions  a  =  U  +  2c  and 
jS  =  U  -  2c,  the  characteristic  equations  are 

—  +  (U  +  c)—  =  ^  (3.16) 

dt  dx  dx[  p  )  p 

+  (7)  -  _ _ —  -  2^  (3.17) 

dt  ^  ^  dx  P  J  P 
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Equation  3.17,  the  j8-equation,  gives  the  relationship  between  V  and  h  on  the 
offshore  boundary  at  x  =  0. 


The  jS-equation  can  be  written  as 


P  =  (^  +  t/,  +  UJ  -  2c 


(3.18) 


where  Ui  is  the  velocity  associated  with  the  incoming  or  incident  long  wave,  is 
the  velocity  associated  with  the  outgoing  or  reflected  long  wave,  and  17„  is  the 
time-mean  velocity  at  x  =  0.  Again  applying  the  shallow-water  assumptions 


i.  C  U 


(3.19) 


Equation  3.18  becomes 


P  =  ^  I  (C,  -  y  *  t;„  -  2,/ift 


(3.20) 


The  first  two  terms  of  the  binomial  expansion  of  h,  assuming  f ,  +  f ,.  <  <  at 
X  =  0,  is 


(3.21) 


Applying  Equation  3.21,  the  jS-equation  becomes 


V  6"o  ; 


(3.22) 
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The  term  (1  +  (f,-  +  ^J/(2h^y^  can  be  expressed  as  the  geometric  series 
1  ■  (?i  +  neglecting  terms  of  order  and  higher.  The  /3-equation  then 

simplifies  to 


P  =  -2^ 
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(3.23) 


The  reflected  wave  amplitude  is  given  by  solving  Equation  3.23  for 
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where  C(t)  =  -0.5  U„.  The  total  depth,  h,  is  the  sum  of  ho  + 

where  f,.  is  given  by  Equation  3.24 
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and  the  velocity,  using  the  definition  of  /3,  is  given  by 


17  =  P  +  2vP 


(3.26) 


To  implement  the  absorbing-generating  boundary  condition,  the 
Predictor-Corrector  method  is  used  to  solve  for  using  Equation  3. 13.  Next, 
the  values  of  T,-  and  C(t)  at  n+1  are  specified.  Then,  Equation  3.25  is  solved  for 
and  Equation  3.26  is  solved  for  These  values  are  then  used  as  the 
offshore  boundary  conditions  in  the  Predictor-Corrector  scheme  illustrated  in 
Equations  3.4-3.12.  This  form  of  the  absorbing-generating  boundary  condition 
assumes  that  the  incoming  infragravity  wave  is  propagating  in  the  +x  direction.  A 


more  general  formulation  in  given  by  Van  Dongeren  et  al,  (1994).  A  periodicity 
condition  is  applied  to  the  lateral  boundaries.  This  assumes  the  longshore 
wavenumber  is  a  multiple  of  the  grid  width  in  the  y  direction.  For  the  problems 
studied  in  this  thesis,  no  longshore  forcing  is  applied  and  V  is  always  zero.  Also, 
since  we  are  simulating  a  closed  wave  tank,  C(t)  is  zero. 

Solution  oscillation 

For  application  of  the  Predictor-Corrector  model  with  surf  beat  forcing 
(Schafer  and  Svendsen  1988,  Schaffer  1993,  see  also  Chapter  5),  the  solutions  for 
f  and  oscillated  from  grid  point  to  grid  point  in  the  cross-shore  direction.  The 
Schaffer  and  Svendsen  analytic  solution  is  linear  {h  =  with  no  bottom  friction, 
and  the  model  was  run  under  these  same  conditions.  The  oscillations  did  not 
become  unbounded,  but  started  from  zero  at  the  cold  start  and  quickly  reached  a 
bounded  value  (up  to  approximately  30  percent  of  Q^.  The  amplitudes  of  the 
oscillations  were  reduced  slightly  by  reducing  the  grid  spacing,  but  were  not 
eliminated.  Figures  3.1a  and  b  show  Predictor-Corrector  results  for  a  plane  beach 
with  Schaffer  and  Svendsen  surf  beat  forcing. 
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a.  Sample  results  for  f 


Rgure  3.1.  Predictor-Corrector  results  with  surf  beat  forcing  (solid  line  — 
Schaffer-Svendsen  analytic  solution,  dashed  line  -  model  results).  Oscillations  are 
due  to  a  disjoint  numerical  solution  (see  page  31  for  discussion). 
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The  cause  of  this  spatial  oscillation  of  the  solution  is  analogous  to  the 
time  splitting  instability  found  in  the  "leap  frog"  numerical  scheme  (Roach  1982). 

In  the  time-splitting  problem,  the  solution  jumps  back  and  forth  between  two 
solutions  on  alternating  time  steps.  This  occurs  because  the  time  derivatives  use 
function  values  at  time  levels  n-1  and  n+1,  while  spatial  derivative  are  taken  at 
time  level  n,  and  thus  the  solutions  for  even  and  odd  time  steps  are  disjoint. 

The  second-order  Predictor-Corrector  scheme  applied  here  has  similar 
construction  in  the  horizontal  plane.  For  example  consider  the  solution  for  I*,  in 
both  the  predictor  and  corrector  steps,  depends  only  on  f "  and  evaluated 
at  i+1  and  i-1  for  time  levels  n-1,  n,  n+1,  and  *.  Thus,  at  even  spatial  steps  f  is 
disjoint  from  the  odd  spatial  step  solution  for  J ,  and  vice  versa.  Two  disjoint 
solutions  can  develop  and  the  values  of  J  and  Q,  oscillate  between  these  solutions. 
Figure  3.2  illustrates  the  finite-difference  molecule  for  the  Predictor-Corrector 
method.  When  the  depth  nonlinearity  is  included,  the  oscillation  disappears 
because  the  nonlinear  term  links  solutions  for  f  and  at  steps  i-1,  i,  and  i+L 

3.1.2  Lax-Wendroff  model 

As  an  alternate  solution  scheme,  the  two-step  Lax-Wendroff  method 
was  applied.  For  simplicity,  the  model  was  developed  in  only  one  horizontal 
direction  (Equations  3.1  and  3.2  only).  The  first  step  of  the  solution  is  given  by 


a.  Predictor  step  b.  Corrector  stq) 

Figure  3.2.  Predictor-Corrector  finite-difference  molecule 
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The  second  step  is: 
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Figure  3.3  shows  the  finite  difference  molecule  for  the  Lax-Wendroff  method  (for 
Equations  3.27  and  3.29).  Note  that  the  nonlinearity  in  water  depth  and  bottom 


firiction  have  been  neglected  in  Equations  3.27-3.30.  This  is  because  this  scheme 
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is  applied  only  to  the  Schaffer-Svendsen  surf  beat  solution,  which  also  neglects 
these  ^components. 
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b.  Step  2 


Hgure  3.3.  Lax-Wendroff  finite-difference  molecule 


The  boundary  conditions  for  the  Lax-Wendroff  model  are  the  same  as 
used  in  the  Predictor-Corrector  model:  an  absorbing-generating  offshore  boundary 
condition  and  a  no-flux  condition  at  the  beach.  The  Lax-Wendroff  model  is  used 
for  applications  using  the  Sch^fer-Svendsen  surf  beat  forcing  in  Chapter  4.  The 
Predictor-Corrector  model  is  used  for  applications  to  the  SXJPERTANK  data  in 
Chapter  5. 

3.2  Input/Output  Parameters 

Both  the  Predictor-Corrector  and  Lax-Wendroff  models  require  the 
same  input  and  produce  the  same  output.  The  input  includes  the  forcing  functions 
of  the  incident  infragravity  wave  elevations  and  the  radiation  stress  forcing 


throughout  the  grid.  A  bottom  friction  coefficient  and  some  short-wave 
information  are  also  required  to  calculate  bottom  friction,  when  bottom  friction  is 
included.  Parameters  must  also  be  specified  to  define  the  grid  domain.  The  output 
includes  the  infragravity  wave  surface  elevations  and  depth-integrated  velocities 
throughout  the  grid. 

3.2.1  Input  parameters 
Incident  infragravitv  surface  elevation 

The  absorbing-generating  boundary  condition  at  the  offshore  boundary 
requires  the  input  of  the  time  series  of  the  incident  infragravity  wave  elevation  f 
This  boundary  condition  can  be  specified  from  an  analytic  solution  or  from 
measurements.  If  measurements  are  used,  the  incident  and  reflected  wave 
components  must  be  separated  (see  Hughes  1993  and  Section  6.3.3). 

Radiation  stress 

A  time  series  of  the  radiation  stress  is  required  at  each  grid  point.  This 
is  the  most  difficult  input  to  supply  to  the  model  because  a  general  model  to 
predict  the  temporal  and  ^atial  radiation  stress  variation,  commensurate  with  the 
level  of  effort  employed  in  the  depth-integrated  model,  is  not  available.  For 
applications  in  this  thesis,  radiation  stress  was  calculated/interpolated  from  wave 
measurements  or  calculated  using  the  analytic  model  of  Schaffer  and  Svendsen 
(1988).  Radiation  stress  is  calculated  based  on  linear  theory  as 


5  =  _  l]  (3.31) 
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where  H  is  wave  height,  is  wave  group  celerity,  and  C  is  wave  celerity.  For 
SUPERTANK  applications,  shallow  water  is  also  assumed,  so  =  C.  The  wave 
height  varies  as  a  function  of  x  and  t. 

Bottom  friction 

The  linear  bottom  friction  formulation  is  given  in  Equation  2.37.  This 
formulation  requires  an  input  bottom  ftiction  coefficient  and  an  estimate  of  the 
short-wave  velocity  amplitude  at  the  bottom  u^.  The  short-wave  velocity  amplitude 
is  estimated  for  each  grid  point  based  on  linear  theory,  an  input  short-wave 
amplitude  and  period  (treated  as  a  constant),  and  the  local  stUl-water  depth.  The 
bottom  Motion  coefficient  is  a  function  of  bottom  roughness  and  bottom  velocity 
(laminar  or  smooth  turbulent  boundary  layer).  Typical  bottom  friction  coefficients 
used  in  previous  steady-state  current  studies  range  from  approximately  0.005  to 
0.05,  based  on  the  Equation  2.37  (e.g.,  Jonsson  1966,  Thornton  1970,  Noda  1974, 
Ebersole  and  Dalrymple  1979,  Thornton  and  Guza  1981,  and  Smith  et  al.  1992). 
Coefficient  values  of  0.01  and  0.05  were  applied  with  the  SUPERTANK  data  with 
minor  variations  in  the  results.  The  bottom  Motion  coefficient  in  both  the  depth- 
integrated  and  profile  models  is  treated  as  a  constant,  and  a  value  of  0.01  is  used. 
The  coefficient  could  be  used  as  a  calibration  parameter. 
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Grid  domain 


;  The  grid  domain  is  specified  by  the  number  of  cells  in  the  x  and  y 
coordinates,  nx  and  ny,  and  the  cell  spacings,  Ax  and  Ay.  The  number  of  time 
stq)s  Tit  and  the  step  increment  At  must  also  be  specified.  The  still-water  depth  ho 
must  be  specified  throughout  the  domain. 

3.2.2  Output  parameters 

The  model  output  parameters  are  the  time  series  of  f  and  U  for  each 
grid  point.  These  output  parameters  provide  the  forcing  for  the  profile  model  to 
extract  the  distribution  of  U  over  depth.  The  reflected  infragravity  wave  at  the 
offshore  boundary  can  also  be  estimated. 

3.3  Example 

As  an  example  model  application  illustrating  the  absorbing-generating 
boundary  condition,  a  standing  wave  is  generated  with  no  radiation  stress  forcing 
and  no  bottom  friction.  The  Predictor-Corrector  model  is  used  for  this  example. 
The  offshore  boundary  condition  is  given  by 

Cj  =  a.cos(kx  -  at)  (3.32) 

where  a,-  is  the  incident  infragravity  wave  amplitude,  k  is  wave  number,  and  a  is 
wave  frequency.  The  parameter  values  specified  are  a,-  =  0.05  m  and  wave  length 
L  =  50  m  (from  which  k  and  a  are  calculated  using  linear  theory).  The  domain  is 
100  m  in  length  and  is  specified  with  nx  =  101,  ny  =  5,  Ax  =  Ay  =  1  m,  and  At 


=  0.25  sec.  A  constant  depth  of  ho  =  0.5  m  is  applied  over  the  entire  grid.  The 
incident  infragravity  wave  input  was  linearly  ramped  up  over  the  first  20  time 
steps. 

The  analytic  solution  for  this  case  is  given  by 

C  =  aJcos(ifcc  -  of)  +  cos(Jbc  +  of)]  (3.33) 


and 

U  =  fl.-^[cos(foc  -  of)  -  cosQac  +  of)]  (3.34) 

kho 

The  solution  has  antinodes  in  f  at  each  of  the  boundaries.  The  maximum  water 
surface  amplitude  is  0.05  m  at  the  antinodes,  and  the  maximum  horizontal  velocity 
is  0.44  m/sec  at  the  nodes.  Example  plots  of  the  results  at  cell  11  and  51  are 
given  in  Figures  3.4  and  3.5.  The  plots  show  a  flat  portion  at  the  beginning  of  the 
record,  prior  to  the  wave  reaching  the  cell  from  a  cold  start.  Then,  the  incident 
wave  reaches  the  cell,  and  finally  the  wave  reflected  from  the  end  of  the  grid 
returns  to  the  cell.  A  close  match  to  the  analytic  solution  is  reached  after  about 
100  sec  at  cell  11  and  after  80  sec  at  cell  51.  The  small  errors  in  amplitude  and 
phase  are  due  to  a  small  re-reflection  at  the  absorbing-generating  boundary. 

Root-mean-square  (rms)  errors  in  the  calculated  solution  for  time  steps 
500-1000  (f  =  125-250  sec)  are  given  in  Table  3.1,  for  selected  grid  cells.  The 
rms  error  is  defined 


# 
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(3.35) 


2  («„  - 


where  a  represents  the  variable  (7  or  lo,  the  subscript  a  denotes  the  analytic 
solution,  n  denotes  the  numerical  solution,  and  N  is  the  number  of  points  in  the 
summation.  The  relative  error  is  expressed  as  the  rms  error  divided  by  the  local 
maximum  of  the  analytic  solution,  rJlmax  and  %  rJU^. 


Table  3.1.  rms  error  in  standing  wave  model  computations 


Grid  CeU 


1  2  11  26  41  51  61  76  91  100  101 


0.0  0.02  0.2  0.5  0.8  1.0  1.2  1.5  1.8  1.98  2.0 


Bl 


.0076  .0077  .0039  .0076  .0036  .0077  .0035  .0078  .0035  .0079  .0078 


Ev.m,  (m/sec) 


■ 

i 


.6  7,8  12.6  7,6  11.6 


11.3  7.8  11.3  7.8  7.8 


.0135  .0134  .0366  .0113  .0338  .0086  .0337  .0057  .0377  .0067  .0 


00  24.4  8.0  00  8.1  00  8.0  oo  8.0  12.2  0. 


Figure  3.4.  Example  standing  wave  results  at  cell  11  (dashed  line  --  analytic 
solution,  solid  line  --  model  results) 
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Figure  3.5.  Example  standing  wave  results  at  cell  51  (dashed  line  --  analytic 
solution,  solid  line  --  model  results) 
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Chapter  4 

PROFILE  MODEL  FOR  INFRAGRAVITY  WAVES 

This  chapter  describes  the  model  component,  called  the  profile  model, 
used  to  calculate  the  vertical  variations  in  the  time-varying  infiagravity  wave 
velocity  field.  The  water  column  is  divided  into  three  regions:  a  bottom  boundary 
layer,  a  central  layer  (from  the  bottom  boundary  layer  to  the  mean  water  surface), 
and  an  upper  layer  (above  the  mean  water  surface).  The  velocity  profiles  are 
solved  in  the  central  layer  with  the  boundary  conditions  of  bottom  and  surface 
stresses  representing  bottom  and  upper  layers.  This  chapter  provides  discussion  of 
the  model  assumptions,  solution  technique,  input  and  output  parameters,  and 
example  applications  for  the  profile  model. 

The  profile  model  is  based  on  the  cross-shore  momentum  equation 
given  in  Equation  2.44.  The  only  additional  assumption  applied  is  that  the  term 
dlf^/dx  is  neglected,  as  was  assumed  for  the  depth-integrated  model  (see 
Chapter  3).  Equation  2.44  then  becomes 

-  ±L  iE]  =  -aK  -  ~  (4.1) 

dt  *  dz)  dx  dx 
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The  right-hand  side  of  Equation  4.1  is  the  forcing  function,  defined  ^&f(t).  The 
bottom,  surface,  and  initial  boundary  conditions  are; 


=  0 


at  z  =  -K 


(4.2) 


iE  =  2^ 

dz  pv, 


at  z  =  C 


(4.3) 


and 

U  =  Ujiz)  at  t  =  0  (4.4) 

where  r/t)  is  the  time-varying  stress  at  the  mean  water  surface  and  UJz)  is  the 
initial  condition.  The  surface  stress  is  due  to  radiation  stress  above  mean  water 
level  and  wind  stress  (which  is  neglected  for  the  applications  in  this  thesis).  The 
eddy  viscosity  is  assumed  constant  over  depth. 

4.1  Eigenfunction  Expansion 

Equation  4.1  is  solved  using  separation  of  variables  and  the  Sturm- 
Liouville  approach.  This  approach  involves  determining  the  solution  in  two  parts: 
u  which  satisfies  the  boundary  conditions  in  z,  and  v  which  is  broken  into  a 
homogeneous  solution  (for  /  =  0)  and  a  particular  solution,  v^. 


4.1.1  Solution  satisfying  boundary  conditions  in  z 

•  The  first  step  is  to  choose  a  partial  solution  u(z,t)  which  satisfies  the 
boundary  conditions  in  z 


u(U)  =  — 


K*  z  + 


fc^^O, 


(4.5) 


The  solution  for  Equation  4.1  then  becomes 

J7(z,f)  =  M(z,f)  +  v(z/)  (^'^ 


4.1.2  Homogeneous  solution 

Combining  Equations  4.5  and  4.6  and  substituting  into  4.1,  gives  the 
governing  equation  for  v 


^  -  V 

at  ” 


= m 


(4.7) 


where 


m  -  AO 

ot 


J_£Z£ 

pv,  at 


fcw^o 


(4.8) 


The  boundary  conditions  for  Equation  4.7  are  given  by  substituting  Equations  4.5 
and  4.6  into  the  boundary  conditions  (Equations  4.2-4)  for  Equation  4.1.  The 
boundary  conditions  are  given  by 


(4.9) 


and 


—  =  0  at  z  =  C 
dz 

pv,-^  -  -fcw^o^  =  0  at  z  =  -/to  <4-10) 

OZ  ’  1C 


v(z,0)  =  u;(z)  t  =  0 


(4.11) 


The  homogeneous  equation  for  v 


0v  _  3^ 


=  0 


(4.12) 


is  solved  by  the  method  of  separation  of  variables,  in  the  form  v  =  Z(z)T(t),  where 
Z  and  T  are  to  be  determined.  Substituting  the  assumed  form  of  v  into  Equation 
4,12  gives 

=  0  (4.13) 

dt  *dz^ 


To  separate  the  variables,  Equation  4.13  is  divided  by  v,  ZT  and  the  first  term  is 
moved  to  the  right-hand  side 
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(4.14) 


dz^ 


M... 


Since  z  and  t  are  independent,  both  sides  of  Equation  4. 14  must  be  equal  to  a 
constant,  which  is  defined  as  V.  The  separation  procedure  produces  two  ordinary 
differential  equations 


—  +  TC^Z  =  0 

dz'^ 


(4.15) 


—  +  VjK^r  =  0 

dt  ' 


(4.16) 


The  general  solutions  for  Equations  4.15  and  4.16  are  given  by 


{IcosKZ  +  JSsinKz 

C  +  Dz 


(4.17) 


(4.18) 


where  A,  B,  C,  D,  E,  and  F  are  constants.  The  general  solution  for  v  is 


V  =  ZT  =  (iicosKz  +  BsmKz)e  ‘  +  C  +  Dz 


(4.19) 
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where  the  constant  coefficients  have  been  redefined  {E  has  been  absorbed  into  A 
and  B,  and  F  has  been  absorbed  into  C  and  D).  The  solutions  for  0  and 
1^  =  0  have  been  summed. 

To  simplify  Equation  4.19,  we  define  ?  =  z  -  T-  Then  Equation  4.19 

becomes 

V  =  ZT  =  (ilcosK^  +  BsinKOe'"''’^'  +  C  +  Dl 

The  boundary  conditions  are  applied  to  determine  the  coefficient  values.  At 

?  =0, 

=  0  =  +  D  (4.21) 

dl 


Since  the  first  and  second  term  on  the  right-hand  side  of  Equation  4.21  are  linearly 
independent,  JB  =  D  =  0.  At  ^  =  -(ho  +  f  ), 


^  P  X  n 

pv, —  -  —f^uv  =  0 


pVj(-AKsmK^  -  -^/^«^(Acosk5  +  0=0 

7t 


(4.22) 


Again,  due  to  linear  independence  of  solution,  C  =  0  and 


-v,K  sinK^ 


^  COSK^ 
1C 


(4.23) 


or,  the  characteristic  equation  for  k  is 


=  1 

fcw^o  tanK(/tj,  +  0 

The  comprehensive  homogeneous  solution  for  v  is 

v(z,f)  =  Ei4,cosTc„(z-C) 

»=1 


(4.24) 


(4.25) 


4.1.3  Particular  solution 

Returning  to  Equation  4.7  to  determine  the  particular  solution  for  v,  we 
first  e?q)ress  the  forcing  term  as  an  eigenfunction  expansion.  We  define 


FJit)  = 


I  dir  It  V 

</ - ^(*0  +  Z  +  7— COSK„(z  -  0  > 

pv,  ^ ^ _  (4.26) 

<  COSK„(z  -  C),  C0SK^(Z  -  C)  ) 


where  the  triangular  brackets  represent  the  inner  product  of  the  two  enclosed 
functions  (Greenberg  1988),  then 


N 


r(0  =  I  Fn(f>  C0SK(Z  -  0 

B=1 


(4.27) 


We  determine  the  particular  solution  by  the  method  of  variation  of  parameters. 
Starting  from 
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Substituting  into  Equation  4.7  yields 


dA 


dt 


-COSTC„(z  -  -  V,Ki^„COSK„(z  -  +  V^lC^i4^COSK„(z  - 


=  F^cosk^Cz  -  0 


(4.29) 


which  simplifies  to 


=  F  (4.30) 

dt 


or 


(4.31) 


Substituting  Equation  4.31  into  4.28  provides  the  particular  solution 


S  /  cosK^(z  -  Q 

n~l  0 


(4.32) 


The  complete  solution  for  v  is  the  sum  the  homogeneous  and  particular  solutions 


4.1.4  Full  solution  for  velocity  profile 

Returning  to  Equation  4.6,  the  solution  for  the  cross-shore  velocity 
profile  is  the  sum  of  Equations  4.5  and  4.33 


U  =  — 

pv,l 


fcw\, 


N 


+  +  /F,e'’'*"^*'*^^iT]cosK^(z-C) 

11=1  0 


-^.(4.34) 


4.1.5  /  independent  of  z 

It  is  convenient  to  consider  the  case  of  /,  the  forcing  function  on  the 
right  hand  side  of  Equation  4.1,  independent  of  z.  The  first  term  on  the  right  hand 
side  of  Equation  4.1  is  independent  of  z,  and  the  second  term  has  a  weak 
dependence  in  shallow  water.  The  assumption  of  /independent  of  z  is  used  often 
in  the  literature  (Svendsen  and  Lorenz  1989,  Svendsen  and  Putrevu  1990,  Sanchez- 
ArcUla  et  al.  1992,  and  Putrevu  and  Svendsen  1993,1995)  and  simplifies  the 
analytical  integration  in  Equation  4.26.  The  solution  for  F„,  assuming /is  not  a 
function  of  z,  is  given  by  the  solution  of  Equation  4.26 


if-^‘ 


n 


^  pfcw^c 


sinK„(*(,  +  C)  + 


1 

'cOSK^iflQ  +  O  -  1 

pv,  dt 

_ _ , 

(4.35) 


The  solution  for  is  given  by 


.rr  ^,(0)  COSK„(/»0  +  Q  '  1 

(17„  -  -^- - )sinK„(Ao  +  C)  +  ^ - - - 


Pfcw^c 


(4.36) 


y(Ao  +  0  +  ■Jsin2K„(;io  +  0 


The  velocity  distribution  over  the  vertical,  given  by  Equation  4.34, 
contains  a  time  integration.  An  efficient  method  is  required  to  evaluate  the 
integral.  Using  integration  by  parts,  the  integral  becomes: 


=  F.(r)£:^7:  -  (4.37) 


)o 


where  the  prime  denotes  a  time  derivative.  The  integral  in  Equation  4.37,  can  be 
broken  into  two  integrals 

/F„(T)c'''’^^"'‘^dT  =  7  (4.38) 


Pulling  out  the  term  expC-i',  K„^At),  gives 


If  we  define  the  integral  /„  as 


t  2 


i,m  =  /.-’■“■‘"'V.w* 


(4.40) 


then,  Equation  4.39  becomes 
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(4.41) 


»-Af 


Partial  integration  can  be  used  repeatedly  to  evaluate  the  second  term  on  the  right- 
hand  side  of  Equation  4.41,  giving  the  following  series 


r-Ai 


2 


Fnit)  ^  F'jO)  ^ 

(4.42) 


where  k  is  0(1)  or  larger  and  v,  is  0(0.01).  Assuming  that  F„  varies  slowly  with 


time  (F„’  «  0),  Equation  4.42  is  substituted  for  the  integral  in  Equation  4.41  to 


give 


/,«)  =  «  AO  4 


2 


(4.43) 


The  analytic  profile  model  is  given  by  Equations  4.34,  4.35,  4.36,  and 
4.43.  The  profile  model  is  also  used  within  the  depth-integrated  model  to  define 
the  bottom  velocity,  V(-ho},  for  calculation  of  the  bottom  friction.  It  is  important 
that  the  forcing/resisting  functions  (radiation  stress  and  bottom  friction)  are  the 
same  in  both  models.  If  they  are  not,  then  f  from  the  depth-integrated  model  will 
not  be  in  equilibrium  with  the  radiation  stress  forcing  and  bottom  friction  in  the 
profile  model.  Disequilibrium  between  models  causes  erroneous  profile  solutions 
and  comparable  solutions  could  not  be  obtained  in  the  two  models. 


4.2  Input/Output  Parameters 

*  The  input  forcing  functions  for  the  profile  model  are  water  surface 
gradient  from  the  depth-integrated  model,  radiation  stress,  and  bottom  friction. 
Additional  input  includes  "grid"-like  parameters  defining  the  accuracy  and 
resolution  of  the  solution.  The  eddy  viscosity  coefficient  is  also  required  as  input. 
Unlike  the  depth-integrated  model,  the  input  is  not  required  globally,  but  is  only 
required  in  the  region  where  the  profile  solution  is  desired.  The  model  output  is 
the  cross-shore  velocity  at  specified  locations  over  the  vertical. 

4.2.1  Input  parameters 
Water  surface  gradient 

The  infragravity  wave  water  surface,  7,  is  obtained  from  the  output  of 
the  depth-integrated  model.  The  gradient  of  f  is  calculated  using  second-order, 
central  finite  differences.  At  the  offshore  and  onshore  boundaries  of  the  depth- 
integrated  model,  forward  and  backward  differences  are  used. 

Radiation  stress 

Radiation  stress  is  divided  into  two  components  to  apply  to  the  profile 
model:  (a)  a  surface  stress  and  (b)  a  depth-uniform  component.  The  radiation 
stress  is  defined  by 
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(4.44) 


c 

=  / ip4  *  P)dz  -  +  Cf 

(Mei  1989).  Substituting  P  =  Pr,  +  pg(^  -  z)  and  i?  =  f  f  into  Equation  4.44 
yields 

=  /  (P“w  +  Pd)*  + 

where  =  -wj,  assuming  linear  theory,  or 

=  /(p»;  - 

Applying  the  linear  theory  definitions  of  u„,  w^,  and  rj,  integrating  over  depth,  and 
taking  the  time  average  reduces  Equation  4.46  to  Equation  3,31  (the  radiation 
stress  applied  to  the  depth-integrated  model).  The  cross-shore  gradient  of  the  first 
term  on  the  right-hand  side  of  Equation  4.46  is  equivalent  to  the  second  term  on 
the  right-hand  side  of  Equation  4.1.  Assuming  shallow  water,  this  term  is  uniform 
over  depth  and  is  calculated  in  the  model  as 

^  ££  ^  (4.47) 

dx  Sh  dx 

assuming  a  mild  slope.  The  surface  stress  term  is  calculated  as 
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Using  Equations  4.47  and  4.48,  one-third  of  the  stress  is  applied  as  the  surface 
stress  and  two-thirds  is  applied  uniformly  over  depth  below  the  mean  water 
surface.  Wave  height  gradients  are  calculated  with  second-order,  central  finite 
differences  (forward  and  backward  differences  at  the  boundaries). 

Bottom  friction 

The  bottom  friction  is  specified  as  it  was  for  the  depth-integrated 
model.  The  formulation  is  given  in  Equation  2.37.  A  constant  bottom  friction 
coefficient  of  0.01  is  used.  Choice  of  the  bottom  friction  coefficient  was  discussed 
in  Section  3.2.1.  The  same  coefficient  is  used  in  both  the  depth-integrated  and 
profile  models. 

"Grid"-like  parameters 

The  number  of  solution  points  over  depth  must  be  specified.  This 
determines  the  resolution  of  the  solution.  The  x-y  locations  from  the  depth- 
integrated  grid  where  current  profiles  are  to  be  calculated  must  also  be  specified. 
The  number  of  eigenfunctions,  N,  to  be  used  in  the  solution  is  an  additional  input 
parameter.  The  number  of  eigenfunctions  determines  the  accuracy  of  the  solution. 
The  resolution  (number  of  points  over  depth)  and  accuracy  (number  of 
eigenfunctions)  are  unrelated  in  the  analytic  solution,  unlike  a  numerical  solution. 


Eddy  viscosity 


‘  The  eddy  viscosity  has  been  assumed  constant  in  the  derivation  of  the 
profile  model.  A  depth-varying  eddy  viscosity  would  require  a  numerical  solution 
over  the  vertical  instead  of  the  present  analytic  solution,  Svendsen  and  Putrevu’s 
(1994)  eddy  viscosity  parameterization  is  used  for  all  cases  studied  in  this  thesis 


'’r  = 


0.01*0^/^ 


0.8 


(hA-^ 


\h) 


+  0.2 


•tb 


inside  surf  zone 
outside  surf  zone 


(4.49) 


where  h,,  is  the  depth  at  incipient  brealdng  and  =  0.01  hj,  (g  The 
breaking  depth  is  estimated  as  where  is  the  incident  zero- 

moment  wave  height. 


4.2.2  Output  parameters 

The  model  output  is  the  time  series  of  U  at  specified  elevations  and  x-y 

locations. 

4.3  Example 

As  an  example  model  application,  steady  short-wave  forcing  conditions 
are  applied.  The  steady  forcing  case  includes  the  processes  of  setdown  seaward  of 
the  breakpoint,  setup  in  the  surf  zone,  and  undertow.  This  example  is  selected 
because  it  is  important  to  show  that  the  depth-integrated  and  profile  models  are,  at 
least  approximately,  in  equilibrium,  i.e.,  produce  the  same  mean  velocities  for 
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steady  forcing.  The  profile  model  forcing  includes  the  water  surface  slope  from 
the  depth-integrated  model  and  the  radiation  stress  forcing.  If  the  profile  model 
forcing  differs  from  the  depth-integrated  model  forcing,  then  the  surface  slope 
gradient  will  not  balance  the  radiation  stress  forcing.  In  a  steady  case,  this  lack  of 
equilibrium  is  highlighted  by  a  profile  solution  that  diverges  from  the  depth- 
integrated  solution.  Also,  for  comparison,  a  simpler  steady  analytic  solution  for 
the  velocity  profile  is  derived  for  comparison. 

The  steady-state  analytic  solution  is  derived  from  Equation  4.1, 
neglecting  the  first  term  on  the  left-hand  side  of  the  equation,  dV/bt.  Integrating 
the  equation  twice  gives 

U  .  *  C,z  *  C,  (4.50) 

where /is  the  forcing  function  given  by  the  right-hand  side  of  Equation  4.1.  Q 
and  Cj  are  integration  constants.  Applying  the  boundary  conditions  given  in 
Equations  4.2  and  4.3,  the  constants  are  found  to  be 


py  V, 


(4.52) 


Integrating  Equation  4.50  over  depth  gives  the  mean  velocity 


V  =  :!i  ilA  +  ^  +  /  ILL^  +  _!_(*  .  0  (4.53) 

P  2v,  f^u\  3v, 


The  depth-integrated  model  gives  the  same  steady-state  solution,  neglecting  the 
terms  containing  the  eddy  viscosity. 

The  example  case  is  a  beach  configured  with  a  1-m  deep,  flat  portion 
extending  3.66  m  and  then  a  1:20  plane  slope.  The  depth-integrated,  Predictor- 
Corrector  model  was  set  up  with  m  =  31,  ny  =  5,  Ax  =  Ay  =  0.6096  m,  and  At 
=  0. 125  sec.  The  minimum  water  depth  was  0. 13  m.  With  steady  forcing,  the 
cross-shore  depth-integrated  velocity  is  zero.  The  short-wave  transformation,  used 
to  calculate  the  radiation  stress  forcing,  was  calculated  for  a  monochromatic  wave, 
using  linear  shoaling  and  the  Dally  et  al.  (1985)  breaker/reformation  model  with  an 
incipient  breaker  index  of  0.8.  The  forcing  was  ramped  up  over  400  time  steps 
(50  sec).  The  profile  model  was  applied  using  20  eigenvectors.  Solutions  were 
generated  at  21  points  over  the  vertical.  The  still- water  level  was  used  in  place  of 
the  mean  water  level  in  the  calculations. 

The  results  of  the  depth-integrated  model  are  given  in  Figure  4.1.  The 
time  series  of  f  and  U  are  given  for  three  grid  locations:  the  offshore  boundary,  x 
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=  0.0  m,  mid  surf  zone,  x  =  11.0  m,  and  near  the  onshore  boundary,  x  = 

21.9  m.  The  effects  of  the  cold  start  and  the  50-sec  ramping  are  apparent  in  the 
first  60-70  sec,  then  the  variables  reach  steady  conditions. 

The  profile  model  results  are  given  in  Figure  4.2a  and  b.  The  vertical 
axis,  Z,  is  dimensionless  depth,  (ho  +  z)/ho.  The  plots  show  that  the  profile  model 
solution  (dotted  line)  oscillates  around  the  steady  profile  (solid  line)  given  by 
Equations  4.50-4.52  until  about  a  time  of  90  sec.  Small  amplitude  oscillations 
continue  until  about  400  sec.  The  rms  errors  for  selected  grid  cells  for  t  —  400- 
750  sec  are  given  in  Table  3.2.  The  overall  rms  error  was  0.031  m/sec  for  t  = 
200-750  sec  and  0.026  m/sec  for  t  =  400-750  sec. 


Table  4.1  rms  error  in  plane-beach  velocity  profile  computations 

Grid  cell 

1 

7 

13 

19 

25 

31 

37 

Eu.rms  (m/sec) 

0.067 

0.0085 

0.051 

0.030 

0.012 

0.0066 

0.0058 

figure  4.1.  Predictor-Corr^tor  depth-integra^  results  for  plane  beach  with 
steady  forcing  (solid  line  —  dashed  line  --  U)  discussed  in  Section  4.3. 
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figure  4.2a.  Hgenfunction  velocity  profiles  for  a  plane  beach  with  steady  forcing 
(solid  line  —  analytic  steady  solution,  dashed  line  -  depth-integrated  model  results, 
dotted  line  -  Eigenfunction  solution)  discussed  in  Section  4.3 
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Figure  4.2b.  Eigenfunction  velocity  profiles  for  a  plane  beach  with  steady  forcing 
(solid  line  --  analytic  steady  solution,  dashed  line  --  depth-integrated  model  results, 
dotted  line  --  Eigenfunction  solution)  discussed  in  Section  4.3 


Chapter  5 

MODEL  APPLICATIONS  WITH  ANALYTICAL  SURF  BEAT  FORCING 

The  purpose  of  this  chapter  is  to  investigate  the  models  developed  in 
Chapters  3  and  4  using  surf  beat  forcing  from  the  analytic  model  of  Schaffer  and 
Svendsen  (1988)  and  Schaffer  (1990,  1993).  The  depth-integrated  model  is 
compared  directly  to  the  Schaffer-Svendsen  model.  The  profile  model  is  also  run 
for  the  Schaffer-Svendsen  surf  beat  forcing.  Profile  model  results  are  compared  to 
the  periodic  infragravity  wave  profile  model  of  Putrevu  and  Svendsen  (1995). 

5.1  Schaffer  and  Svendsen  Surf  Beat  Model 

Schaffer  and  Svendsen  (1988)  and  Schaffer  (1990)  developed  a 
theoretical  two-dimensional  model  for  surf  beat  generation  by  incident  wave  groups 
on  a  plane  beach.  The  surf  beat  is  forced  by  the  modulation  of  radiation  stress. 
Short  waves  are  modulated  both  inside  and  outside  the  surf  zone,  i.e.,  wave 
breaking  does  not  destroy  the  wave  groups.  The  break  point  in  this  model  was 
held  constant.  Schaffer  (1993)  extended  the  model  to  include  both  a  constant  break 
point  with  wave  modulation  in  the  surf  zone  and  a  varying  break  point  position 
with  no  modulation  in  the  surf  zone,  as  suggested  by  Symonds  et  al.  (1982). 


5.1.1  Surf  beat  model  description 

>  The  governing  equations  for  the  two-dimensional  Schaffer-Svendsen 
surf  beat  model  are  Equations  3.1  and  3.2.  The  model  neglects  the  bottom  friction 
term  and  the  depth  nonlinearity  in  the  momentum  equation.  The  two  equations 
were  cross  differentiated  (d/dx  was  taken  of  Equation  3.1  and  d/dt  was  taken  of 
Equation  3.2)  and  combined  to  eliminate  Q^.  The  resulting  surf  beat  governing 
equation  is 


_3_ 

ar 


.  ar) 


P  dx^ 


(5.1) 


where  the  water  surface  elevation  is  divided  into  a  time-mean  component,  f  and 
a  time-varying  component,  f;, 

C  =  To(x)  +  T[(x^)  (5-2) 

and  likewise  the  radiation  stress  is  divided  into  a  time-mean  component,  S„  o>  and  a 
time-varying  component, 

‘  (5^) 

Mean  values  of  the  variables  have  been  neglected  in  Equation  5.1.  The  onshore 
boundary  condition  is  full  reflection  of  the  infragravity  wave.  The  solution  is 
e^ressed  as 


(5.4) 


Z,  =  aZ?  +  PCf  + 

where  the  superscripts  1  and  2  denote  two  linearly  independent  homogeneous 
solutions,  and  the  superscript  p  denotes  a  particular  solution  which  is  chosen  to 
vanish  at  the  shoreline.  Full  reflection  is  specified  by  a  =  jS.  The  offshore 
boundary  condition  is  specified  as  a  incoming  bound  wave 

,  (5.5) 

P(SA.  -  cl) 

where  is  the  shelf  depth  and  Cg,  is  the  group  celerity  in  the  offshore,  shelf 
region.  At  the  break  point,  matching  conditions  of  continuity  of  mass 

[g:  =  0  (5.6) 

and  momentum 


+ 

1 

dx 

Pgh 

dx  . 

(5.7) 


are  applied,  where  the  square  brackets  indicate  the  difference  in  variables  evaluated 
on  the  +  and  -  x  sides  of  the  matching  location. 

The  short-wave  modulation  is  formed  by  two  monochromatic  waves 
with  slightly  different  frequencies,  =  ca/l  +  e)  and  coj  =  (a/l  -  e).  The  mean 
frequency  is  w,  and  e  is  a  small  perturbation.  The  difference  frequency  is  defined 
as  w  =  0);  -  0)2  and  is  the  frequency  of  the  infragravity  wave.  The  amplitude  of 
the  first  component  is  defined  as  a,  and  the  amplitude  of  the  second  component  is 
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defined  as  ba.  It  is  assumed  that  the  short-wave  height  modulation  is  small, 

5  <  <  1.  The  mean  short-wave  wavenumber  is  k,  and  the  wavenumbers  for  the 
component  waves  =  k  (1  +  e  C/C^ )  and  kj^  k  (1  -  e  C/C^ ),  where  C  and 
Cg  are  the  wave  celerity  and  group  celerity,  respectively.  The  tinie  varying 
radiation  stress  is  calculated  from  linear  theory,  which  using  Equation  5.3  yields 


|AiP 


(5.8) 


Neglecting  terms  of  0(5^  and  higher,  the  time  varying  wave  amplitude  is  given  by 


iA,f  = 


-  k)  2  5  exp 
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if-^dx 
J  C 


fl^26exp 
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if-^dx 
j  C 


(x  i  x^ 
(x  i  Xj) 


(5.9) 


where  x  is  measured  from  the  shoreline,  Xj  is  the  location  of  the  break  point  (which 
may  be  time  varying),  y  is  the  breaker  index  (ratio  of  wave  amplitude  to  water 
depth  at  breaking),  and  k  is  the  parameter  which  determines  wave  groupiness  in  the 
surf  zone.  The  integration  in  Equation  5.9  is  calculated  using  trapezoidal 
integration. 

For  a  value  of  k  =  1,  the  time- varying  short-wave  forcing  in  the  surf 
zone  is  zero,  as  seen  from  Equation  5.9.  This  case  is  analogous  to  the  model  of 
Symonds  et  al.  (1982).  The  position  of  the  break  point  varies  with  time,  as 
specified  by  the  limiting  condition  \A\^  =  y^  ho,  where  \A\  includes  both  the 


steady  and  time-varying  short-wave  components.  Within  the  surf  zone,  wave 
amplitude  is  constant  for  each  depth.  For  k  =  0,  the  groupiness  in  the  surf  zone  is 
maintained,  but  the  amplitude  is  reduced  in  proportion  to  water  depth.  This  case  is 
analogous  to  Schaffer  and  Svendsen  (1988)  and  Schaffer  (1990).  Dq)th  at  the 
break  point  is  constant  and  defined  by  h,,  =  a/y.  Visualizations  of  the  short-wave 
modulation  for  jc  =  1  and  0  are  given  in  Figure  5.1.  The  three  regions,  marked  I, 
n,  and  in  in  the  figure,  represent  the  surf  zone,  the  shoaling  region,  and  a 
constant  depth  "shelf  region,  respectively.  The  region  B,  for  ic  =  1,  is  the  zone 
where  incipient  breaking  takes  place. 

The  solution  for  the  infragravity  wave  amplitude  is  given  by 


^(J)  = 
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dx^ 
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Ho’(y) 


2  PS*. 


* 

p«  .  r  Hi(y)^dx 
i  dx^ 


(5.10) 


where  is  the  beach  slope,  7  =  1  and  11  for  the  region  inside  and  outside  the  surf 
zone,  respectively,  i  =  V- 1 ,  Tf/-’  and  are  Hankel  functions  of  the  first  and 
second  kind,  and  the  argument  of  the  Hankel  functions  is  given  by 


and  1  in  regions  I,  H,  and  HI  for  Schaffer-Svendsen  surf  beat  model  (from  Schaffer 


1993) 
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The  integrals  in  Equation  5.10  are  evaluated  using  Simpson  Rule  integration.  The 
radiation  stress  derivatives  are  calculated  using  second-order  central  finite 
differences.  The  amplitude  of  the  infiagravity  wave  in  region  IE,  the  constant 
depth  shelf  region,  is  given  by 


aim 

•’xt.l 


pfe*.  -  ci) 


(5.12) 


The  last  term  in  Equation  5.12  is  the  input  forced  long  wave  (Longuet-Higgins  and 
Stewart  1962,  1964).  This  assumes  dynamic  equilibrium  between  the  radiation 

•  stress  forcing  and  the  forced  long  wave,  i.e.,  a  long  propagation  distance  in  the 

constant  depth  region.  In  very  shallow  water,  the  group  velocity  approaches  the 
free  wave  velocity  and  a  near-resonance  condition  exists. 

The  complex  constants  o;  and  j3  in  Equations  5.10  and  5.12  are 
determined  by  matching  the  solutions  between  regions  I  and  n  and  regions  n  and 

•  m.  At  the  shoreline,  jS  is  set  equal  to  a,  which  equates  the  reflected  and  incident 

infragravity  wave  amplitudes  at  the  shoreline,  as  stated  earlier  this  is  the  condition 
of  full  reflection.  The  time-varying  solution  for  the  infragravity  wave  surface 

^  elevation  is  given  by 


Cj  =  +  ♦) 


(5.13) 
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where  *  denotes  the  complex  conjugate  of  the  preceding  term  and  f  is  given  by 
Equation  5. 10  or  5. 12. 

The  mean  (steady)  radiation  stress  forcing  was  discussed  by  Schaffer 
(1993),  but  was  not  included  in  the  model  solution  (Equations  5.10  and  5.12).  The 
mean  radiation  stress  is  given  by 
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_  0 
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(5.14) 


and  the  mean  water  surface  elevation  is  given  by 


Cn  = 


a^k 


2sinh2Mn 


3y' 


+  1 


iK  -  - 


X  t  X, 


X  X. 


(5.15) 


where  the  subscript  b  represents  variables  evaluated  at  the  break  point.  Equations 
5.14  and  5.15  have  been  added  to  the  model  to  calculate  mean  radiation  stress  and 
mean  water  surface  elevation. 


5.1.2  Model  input  and  output 

The  input  required  to  run  the  surf  beat  model  includes  descriptions  of 
the  two  short-wave  trains,  beach  slope,  shelf  depth,  the  parameter  k,  and 


parameters  defining  the  solution  grid.  The  output  parameters  for  each  grid  cell  are 
water- depth,  infragravity  wave  surface  elevation  and  depth-integrated  velocity,  and 
radiation  stress. 

Input  parameters 

The  input  parameters  defining  the  short-wave  components  are  the  wave 
amplitudes,  a  and  8a,  in  meters  and  the  wave  frequencies,  and  co^,  in  sec'^ 

The  beach  configuration  is  defined  by  specifying  the  shelf  depth  in  meters,  ,  and 
the  beach  slope,  The  computation  grid  is  specified  by  the  nondimensional 
depth  change  between  grid  steps.  Ah,  which  is  defined  as  the  depth  change  Ah  in 
meters  multiplied  by  the  deepwater  wavenumber,  w/Zg.  The  horizontal  grid  step 
and  number  of  grid  points  are  internally  calculated  based  on  beach  slope  and  shelf 
depth.  The  dimensionless  minimum  and  maximum  depths  for  model  output  must 
also  be  specified.  These  depths  are  nondimensionalized  using  the  deepwater 
wavenumber.  Additional  input  parameters  are  the  breaker  index  y,  gravitational 
acceleration  g  (SI  units  will  be  used  throughout  this  study),  and  the  parameter  k. 
Schaffer  (1993)  suggests  that  the  value  of  k  is  near  unity.  Values  of  x  <  1 
represent  continued  wave  grouping  into  the  surf  zone  (see  Equation  5.15).  Values 
of  X  >  1  represent  a  reversal  of  wave  grouping  into  the  surf  zone,  i.e.,  setdown 
under  what  was,  prior  to  breaking,  low  waves  and  setup  under  what  was  high 


waves. 


Output  parameters 


For  each  grid  cell,  the  model  ou^ut  is  the  dimensionless  still-water 
dq)th  ho  the  dimensionless  mean  water  surface  the  dimensionless 

time-varying  water  surface  amplitude  the  dimensionless  mean  radiation 

stress  amplitude  S„  (/(aJb),  the  dimensionless  time-varying  radiation  stress 
amplitude  \S„i\/(aJb),  and  the  dimensionless  time-varying  infragravity  wave 
velocity  amplitude  \  lJi\/( i  a„((ajg)).  a„  is  the  deepwater  wave  amplitude  that  is 
equivalent  to  a. 

Sample  surf  beat  model  output  is  shown  in  Figure  5.2.  This  figure 
reproduces  Figure  5  given  by  Schaffer  (1993).  The  figure  shows  the  spatial 
variation  of  the  time  varying  nondimensional  surf  beat  envelope  (-h/-  the 
amplitude)  as  a  function  of  nondimensional  depth,  ho^i^/g.  The  nondimensional 
surf  beat  for  times  r  =  0  sec  and  -0.25  T  sec,  where  T  is  the  surf  beat  period,  are 
also  shown  in  the  figure.  The  input  parameters  for  this  simulation  are  listed  as 


Case  1  in  Table  5.1. 


Figure  5.2.  Case  1:  nondimensional  surf  beat  envelope  for  Schaffer-Svendsen 
model,  ±  l^^/fSaoo)  i  ”  solid  line,  ^j/(8a„)  at  r  =  0  and  -T/4  —  dashed  and  chain- 
dot  lines 
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Table  5.1.  Surf  beat  model  input 

Case 

a,  m 

ba,  m 

wj,  sec'* 

0)2,  sec'* 

K  m 

K 

7 

■ 

1 

0.62 

0.062 

1.375 

1.125 

18.835 

0.05 

0.4 

1.0 

2 

0.80 

0.08 

2.18 

2.02 

3.0 

0.08 

0.4 

0.0 

3 

0.80 

0.08 

2.25 

1.95 

3.0 

0.08 

0.4 

0.0 

5.1.3  Depth-integrated  model  surf  beat  results 

The  Schaffer-Svendsen  surf  beat  model,  with  the  additional  calculation 
of  mean  radiation  stress  and  water  surface  elevation,  provides  the  input  required  to 
run  the  depth-integrated  model  described  in  Chapter  3.  The  surf  beat  model  also 
provides  an  analytic  solution  to  which  the  depth-integrated  model  results  can  be 
compared.  In  this  section  and  in  Section  5.2,  the  surf  beat  model  output  will  be 
used  to  drive  the  depth-integrated  and  profile  models.  The  periodic  nature  of  the 
solution  provides  the  opportunity  to  investigate  model  input  parameters  and  output 
trends. 

As  discussed  in  Section  3.1.1  of  Chapter  3,  the  Predictor-Corrector 
model  gave  results  which  oscillated  from  grid  point  to  grid  point  in  the  cross-shore 
direction  when  executed  with  surf  beat  input  (see  Figure  3.1).  The  reason  for  this 
oscillation  is  a  spatial  splitting  of  the  solution  at  alternating  grid  points  when 


bottom  friction  and  depth  nonlinearity  are  neglected,  as  they  are  in  the  surf  beat 
model  described  above.  For  the  following  comparisons,  the  Lax-Wendroff  version 
of  the  depth-integrated  model  is  used. 

Three  surf  beat  applications  will  be  used  for  model  comparisons: 

Case  1  is  the  application  shown  in  Figure  5.2,  and  Cases  2  and  3  were  chosen  to 
approximate  natural  beach  conditions  similar  to  those  occurring  during  the 
SUPERTANK  laboratory  project  (see  Chapter  6).  The  surf  beat  ou^ut  of  mean 
and  time-varying  radiation  stress  throughout  the  grid  is  used  as  input  to  the  depth- 
integrated  model  along  with  the  incident  infragravity  wave  surface  elevation  at  the 
offshore  boundary  given  by  Equations  5.12  and  5.15 

^  ^ ^ _  (5.16) 

2sinh2fc/t^  p(g/i,  - 


Case  1 

The  surf  beat  input  conditions  for  Case  1  are  listed  in  Table  5.1.  The 
grid  spacing  used  in  the  surf  beat  and  Lax-Wendroff  models  was  Ax  =  1.88  m 
with  nx  =  200.  The  time  step  used  in  the  Lax-Wendroff  model  was  Ar  = 

0.126  sec.  The  maximum  Courant  number  was  0.91.  The  Lax-Wendroff  model 
was  started  from  a  cold  start  (^  =  0  m  and  fit  =  0  m/sec  over  the  entire  grid). 
The  forcing  functions  (f,  and  were  linearly  ramped  from  zero  to  the  full  values 
specified  by  the  surf  beat  model  over  500  time  steps.  Lax-Wendroff  model  results 


are  given  in  two  forms:  "snapshots"  which  illustrate  the  entire  domain  at  a  given 
time  step  and  time  series  plots  for  given  x  locations.  The  orientations  of  the  x  axis 
in  the  surf  beat  and  Lax-Wendroff  models  are  opposite.  The  example  output  plots 
use  the  x-axis  orientation  of  the  Lax-Wendroff  model  which  is  x  =  0  at  the 
offshore  boundary  and  x  increasing  toward  the  shore. 

Figure  5.2  gives  an  overview  of  the  spatial  variation  of  the  surf  beat 
forcing  for  Case  1.  The  majority  of  the  model  domain  is  outside  the  surf  zone  (the 
breaker  point  is  shown  with  a  vertical  line).  The  figure  shows  six  partial  nodes  in 
the  surface  elevation,  with  the  nodes  located  inside  and  just  outside  the  surf  zone 
having  almost  zero  amplitude.  Antinodes  exist  near  the  breaker  line  and  at  the 
shoreline.  Also,  note  the  discontinuity  in  elevation  at  the  breaker  point.  This  is 
the  narrow  incipient  breaker  zone  formed  by  the  time-varying  break  point  for 
K  =  1.  The  period  of  the  infragravity  wave  is  25  sec.  The  short-wave 
modulation,  5,  is  10  percent  (see  Table  5.1) 

What  is  happening  in  the  surf  beat  solution  is  that  the  short-wave  height 
varies  and  creates  a  long  wave  (with  a  trough  under  high  wave  groups  and  a  crest 
under  low  wave  groups).  The  variation  is  driven  by  two  short-wave  trains  with 
different  amplitudes  and  frequencies.  The  infragravity  wave  frequency  is  equal  to 
the  difference  in  frequency  of  the  short-wave  trains. 

Figures  5.3a  -  e  give  snapshots  of  the  infiragravity  water  surface 
elevations  at  5  sec  intervals  over  one  wave  period.  The  solid  lines  are  the  surf 
beat  analytic  solution  and  the  dashed  lines  are  the  Lax-Wendroff  model  results. 


Figure  5.3.  Case  1  snapshots  of  water  surface  elevation  from  depth-integrated  and 
analytic  Schaffer-Svendsen  models  (solid  line  —  analytic  surf  beat  solution,  dashed 
line  —  numerical  solution) 
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The  times  given  at  the  top  of  the  figures  are  in  seconds,  and  start  well  after  the 
ramp-up  period  (approximately  60  sec).  The  figures  show  small  surface 
oscillations  in  the  shoaling  region,  a  strong  mean  setdown  at  the  break  point,  and 
signiticant  time-varying  oscillation  at  the  break  point  and  shoreline  (anti-nodes  in 
the  surf  beat  solution).  The  dimensionless  mean  and  root-mean-square  (rms)  error 
for  selected  locations  are  given  in  Table  5.2.  The  dimensionless  errors  are  defined 
as 
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where  SS  represents  the  Schaffer-Svendsen  surf  beat  solution,  LW  represents  the 
Lax-Wendroff  solution,  and  Hq  is  the  local  mean  wave  height.  The  error 
calculations  begin  after  1000  time  steps  to  avoid  the  effect  of  the  numerical  ramp- 
up  time.  The  numerical  Lax-Wendroff  results  match  the  analytic  solution  very 
well.  The  mean  errors  increase  in  magnitude  from  less  than  0.(X)1  %  of  Hp 
(0.0001  m)  in  the  offshore  region  to  0.01  m  near  the  shoreline  (where  Hg  goes  to 
zero)  where  both  the  mean  and  time-varying  water  surface  amplitudes  are  largest. 
The  largest  rms  errors  occur  where  the  gradient  in  the  time-varying  envelope  are 
greatest,  e.g.,  in  the  mid  surf  zone,  which  suggests  that  the  errors  are  caused  by 
inaccuracies  in  time  integration,  which  cause  small  phase  errors. 


Table  5.2.  Lax-Wendroff  error  for  Case  1 

X,  m 

82.9 

167.6 

252.4 

344.7 

363.5 

0.001 

0.009 

0.001 

o 

o 

1 

-0.97 

0.25 

0.49 

0.36 

0.98 

0.27 

% 

0.001 

0.001 

0.001 

-0.01 

-0.002 

Ev,ms}  % 

0.06 

0.02 

0.07 

0.07 

0.34 

Figures  5.4a-e  give  time  series  of  the  water  surface  elevations  at  five 
locations.  These  plots  again  show  good  agreement  of  the  numerical  solution  with 
the  analytic  solution.  There  is  a  very  small  error  in  the  period  of  the  numerical 
solution  resulting  in  a  phase  error  developing  in  time. 

Figures  5.5a-e  give  snapshots  of  the  depth-integrated  velocity  for  the 
same  times  shown  in  Figure  5.3.  Again,  the  solid  lines  are  the  analytic  solution 
and  the  dashed  lines  are  the  numerical  solution.  Errors  in  the  computed  velocities 
are  again  given  in  Table  5.2.  The  error  is  defined  by  Equation  5.17.  Time  series 
of  the  depth-integrated  velocity  are  given  in  Figures  5.6a-e.  As  shown  in  the 
water  surface  elevation  time  series,  a  small  phase  error  is  developing. 


Rgure  5.4.  Case  1  time  series  of  water  surface  from  depth-integrated  and  analytic 
Schaffer-Svendsen  models  for  selected  locations  (solid  line  ~  analytic  surf  beat 
solution,  dashed  line  --  numerical  solution) 
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Figure  5.5.  Case  1  snapshots  of  depth-integrated  velocity  from  depth-integrated 
and  analytic  Schaffer-Svendsen  models  (solid  line  -  analytic  surf  beat  solution, 
dashed  line  —  numerical  solution) 


c.  X  =  167.6  m 


Rgure  5.6.  Case  1  time  series  of  depth-integrated  velocity  from  depth-integrated 
and  analytic  Schaffer-Svendsen  models  for  selected  locations  (solid  line  ~  analytic 
surf  beat  solution,  dashed  line  --  numerical  solution) 


•  The  configuration  of  Case  2  is  somewhat  similar  to  the  wave  channel  at 
SUPERTANK,  with  an  offshore  depth  of  approximately  3  m  and  a  surf  zone  width 
on  the  order  of  25  m  (for  no  nearshore  bar).  The  surf  beat  input  parameters  for 
Case  2  are  also  listed  in  Table  5.1.  The  surf  beat  output  is  shown  in  Figure  5.7. 
The  grid  spacing  used  in  both  the  surf  beat  and  Lax-Wendroff  models  was  Ax  = 
0.25  m  and  /IX  =  150.  The  time  step  used  in  the  Lax-Wendroff  model  was  Ar  = 
0.033  sec.  The  maximum  Courant  number  was  0.71. 

The  infragravity  periods  at  SUPERTANK  were  typically  in  the  range  of 
20  to  40  sec.  The  period  of  the  infragravity  wave  for  Case  2  is  39  sec.  A  value 
of  AC  =  0  was  selected  because  wave  groupiness  was  observed  in  the  surf  zone  at 
SUPERTANK.  Remember  that  for  ic  =  1,  the  time-varying  short-wave  amplitude 
in  the  surf  zone  is  zero  and  wave  groupiness  is  destroyed  (Equation  5.9).  For 
<c  =  0,  the  time-varying  amplitude  is  reduced  as  function  of  depth,  but  maintains 
its  grouped  nature.  Note  that  for  k  =  0  (Figure  5.7),  that  the  discontinuity  in  the 
surface  elevation  envelope  is  only  in  the  slope  and  not  in  magnitude,  due  to  the 
constant  break  point.  The  main  feature  of  the  surf  beat  solution  is  a  partial  node  in 
the  mid  surf  zone.  Since  the  envelope  is  fairly  flat,  the  solution  does  not  vary 
greatly  across  the  grid. 

Figures  5.8a-e  give  snapshots  of  the  infragravity  surface  elevations  over 
one  wave  period.  As  expected  from  the  fact  that  the  short-wave  frequency 
perturbation  e  is  small  and  thus  the  infragravity  wave  is  long  compared  to  the 


Figure  5.7.  Case  2:  nondimensional  surf  beat  envelope  for  Schaffer-Svendsen 
model,  ±  l?j/f3aoo)l  —  solid  line,  ^i/(baj)  at  f  =  0  and  -TIA  —  dashed  and  chain- 
dot  lines 
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e.  t  =  64.8  sec 


Figure  5.8.  Case  2  snapshots  of  water  surface  elevation  from  depth-integrated  and 
analytic  Schaffer-Svendsen  models  (solid  line  -  analytic  surf  beat  solution,  dashed 
line  —  numerical  solution) 
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domain,  the  solution  is  fairly  bland,  with  small  oscillations  around  the  mean 
setdown/setup  solution.  The  water  surface  is  rocking  like  a  teeter-totter  around  the 
central  node.  There  is  a  small  mean  bias  in  the  numerical  solution  which  increases 
near  the  shoreline.  Figure  5.9  shows  time  series  for  five  locations.  In  this  case, 
the  surf  beat  amplitude  at  the  shoreline  is  relatively  small.  The  phase  difference 
between  the  solutions  is  also  small,  but  increases  near  the  shoreline.  Mean  and 
rms  errors  in  elevation  are  reported  in  Table  5.3. 

Figures  5.10  and  5.11  give  snapshot  and  time  series  plots  for  the  depth- 
integrated  velocity  for  Case  2.  The  velocities  are  low  and  also  rock  around  the 
central  node  in  the  solution.  Again,  the  largest  errors  occur  near  the  shoreline. 
Errors  are  given  in  Table  5.3. 


Table  5.3.  Lax-Wendroff  error  for  Case  2 

jc,  m 

17.5 

22.5 

28.8 

35.0 

-0.004 

-0.009 

0.01 

0.04 

0.19 

0.82 

0.53 

0.59 

0.75 

1.00 
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-0.01 

-0.007 
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0.87 
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figure  5.9.  Case  2  time  series  of  water  surface  from  depth-integrated  and  analytic 
Schaffer-Svendsen  models  for  selected  locations  (solid  line  -  analytic  surf  beat 
solution,  dashed  line  ~  numerical  solution) 
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e.  t  =  64.8  sec 


figure  5.10.  Case  2  snapshots  depth-integrated  velocity  from  depth-integrated  and 
analytic  Schaffer-Svendsen  models  (solid  line  ~  analytic  surf  beat  solution,  dashed 
line  —  numerical  solution) 
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e.  X  =  28.5  m 


Figure  5.11.  Case  2  time  series  of  depth-integrated  velocity  from  depth-integrated 
and  analytic  Schaffer-Svendsen  models  for  selected  locations  (solid  line  --  analytic 
surf  beat  solution,  dashed  line  —  numerical  solution) 
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Case  3 

The  surf  beat  input  parameters  for  Case  3  are  listed  in  Table  5.1.  The 
surf  beat  output  is  shown  in  Figure  5.12.  The  grid  spacing  used  in  both  the  surf 
beat  and  Lax-Wendroff  models  was  Ax  =  0.25  m  and  nx  =  150.  The  time  step 
used  in  the  Lax-Wendroff  model  was  Ar  =  0.035  sec.  The  maximum  Courant 
number  was  0.76.  The  configuration  of  Case  3  is  similar  to  Case  2,  but  the  period 
of  the  infragravity  wave  is  21  sec,  versus  39  sec  for  Case  2.  Consequently,  the 
infragravity  wave  is  shorter  and  the  surf  beat  solution,  shown  in  Figure  5.12  now 
has  both  a  node  and  anti-node  in  the  surf  zone. 

Figures  5.13a-e  give  snapshots  of  the  infiragravity  surface  elevations 
over  one  wave  period.  The  solution  is  similar  to  Case  2,  although  the  gradients  in 
water  surface  are  greater.  Figure  5.14  shows  time  series  for  five  locations.  The 
time  series  show  good  agreement  in  phase  and  magnitude.  The  errors  in  elevation 
are  reported  in  Table  5.4.  Figures  5.15  and  5.16  give  snapshot  and  time  series 
plots  for  the  depth-integrated  velocity  for  Case  3.  Errors  are  given  in  Table  5.4. 
As  seen  in  Case  2,  errors  in  velocity  increase  near  the  shoreline.  Generally,  the 
agreement  is  good. 
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Elgure  5.12.  Case  3:  nondimensional  surf  beat  envelope  for  Schaffer-Svendsen 
model,  ±  I  |/f3aoo)  |  --  solid  line,  at  r  =  0  and  -r/4  --  dashed  and  chain- 

dot  lines 


e.  t  =  69.1  sec 


Figure  5.13.  Case  3  snapshots  of  water  surface  elevation  from  depth-integrated 
and  analytic  Schaffer-Svendsen  models  (solid  line  -  analytic  surf  beat  solution, 
dashed  line  —  numerical  solution) 
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e.  X  =  28.5  m 


Figure  5.14.  Case  3  time  series  of  water  surface  from  depth-integrated  and 
analytic  Schaffer-Svendsen  models  for  selected  locations  (solid  line  —  analytic  surf 
beat  solution,  dashed  line  -  numerical  solution) 
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e.  t  =  69.1  sec 


Kgure  5.15.  Case  3  snapshots  of  depth-integrated  velocity  from  depth-integrated 
and  analytic  SchMfer-Svendsen  models  (solid  line  —  analytic  surf  beat  solution, 
dashed  line  -  numerical  solution) 
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e.  X  =  28.5  m 


Figure  5.16.  Case  3  time  series  of  depth-integrated  velocity  from  depth-integrated 
and  analytic  Sch^fer-Svendsen  models  for  selected  locations  (solid  line  -  analytic 
surf  beat  solution,  dashed  line  —  numerical  solution) 


94 


Table  5.4.  Lax-Wendroff  error  for  Case  3 

jc,  m 

Wk 

17.5 

22.5 

28.8 

35.0 

0.02 

-0.03 

0.01 

0.03 

0.05 

0.78 

0.39 
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0.43 

0.63 

Ev,.,  % 
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0.02 

0.001 

-0.01 

-0.03 

Ev.n^,  % 

0.62 

0.46 

0.45 

0.53 

0.78 

5.1.4  ProHle  model  surf  beat  results 

The  profile  model,  described  in  Chapter  4,  was  also  run  for  Cases  2 
and  3.  The  model  was  driven  with  output  from  the  Lax-Wendroff  model.  The 
number  of  eigenvectors  used  was  20,  and  the  number  of  points  solved  over  the 
depth  was  21.  Results  for  Cases  2  and  3  are  given  in  Figures  5.17  and  5.18, 
respectively.  Both  figures  show  similar  trends.  The  figures  show  velocity  profiles 
over  an  entire  infragravity  wave  cycle.  Both  mean  and  time-varying  forcing  have 
been  applied. 

In  the  region  outside  the  surf  zone  (x  =  7.25  m),  the  velocity  profiles 
"lean"  offshore  at  the  surface  and  are  fairly  uniform  over  the  lower  two-thirds  of 
the  depth.  This  is  due  to  an  increasing  amplitude  of  the  short  waves,  and  thus  a 


Figure  5.17.  Case  2  surf  beat  velocity  profiles  from  eigenfunction  profile  model 
at  locations  outside  and  through  the  surf  zone 


Figure  5.18.  Case  3  surf  beat  velocity  profiles  from  eigenfunction  profile  model 
at  locations  outside  and  through  the  surf  zone 
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positive  gradient  in  the  radiation  stress.  Since  the  water  depth  is  relatively  deep  in 
this  region,  the  surface  stress  does  not  greatly  impact  the  velocity  distribution  near 
the  bottom. 

In  the  region  from  the  outer  to  mid  surf  zone  (;c  =  17.25  m  and  x  = 
28.5  m),  the  velocity  profiles  exhibit  the  classical  parabolic  undertow  shape  with 
offshore  flow  near  the  bottom  and  onshore  flow  near  the  surface.  The  shape  is 
driven  by  the  breaking-induced  decrease  in  the  radiation  stress.  The  mean 
radiation  stress  gradient  is  dominant.  The  mean  flows  are  zero  since  this  is  a 
closed  system  with  no  inflow  or  outflow.  The  surface  stress  continues  to  dominate 
the  velocity  profiles  in  this  region.  The  surface  stress  is  strong  and  depths  are 
shallower,  so  the  surface  effects  reach  near  the  bottom. 

In  the  very  nearshore  region  (x  =  34.75  m),  the  behavior  of  the 
profiles  changes  again.  In  this  region,  the  time-varying  radiation  stress  forcing  is 
reduced  because  the  short  waves  are  small.  The  dominant  forcing  is  depth  uniform 
and  due  to  the  gradients  in  the  infragravity  wave  elevation  (strong  nodes  and 
antinodes  in  the  inner  surf  zone).  The  resultant  velocity  profiles  are  less  curved. 
The  top  and  bottom  of  the  profiles  also  tend  to  be  in  phase. 

In  the  next  section,  the  mean  forcing  will  be  neglected.  The  time- 
varying  results  will  be  discussed  in  greater  detail. 
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5.2  Putrevu  and  Svendsen  Periodic  Profile  Model 


Putrevu  and  Svendsen  (1995)  developed  an  infiragravity  velocity  profile 
model  for  the  surf  zone  which  assumes  periodicity  in  the  forcing  functions.  The 
velocity  profiles  are  driven  by  radiation  stress  gradients  and  free  and  forced 
infragravity  waves.  The  model  assumptions  and  governing  equations  are  similar  to 
the  profile  model  described  in  Chapter  3.  The  solution  method  is  different  in  that 
it  assumes  periodicity  in  all  driving  forces  (i.e.,  periodic  wave  groups).  Although 
the  periodicity  is  the  same  in  all  driving  forces,  the  relative  phase  of  the  surface 
stress  forcing  and  the  depth  uniform  forcing  can  vary  and  must  be  specified. 

5.2.1  Putrevu-Svendsen  model  description 

The  governing  equation  for  the  Putrevu-Svendsen  profile  model  is 
Equation  4.1.  The  model  was  developed  assuming  a  gently  sloping  beach,  a  long 
straight  coast  (longshore  homogeneity)  with  normal  wave  incidence,  shallow  water 
(pressure  variation  in  the  infragravity  waves  is  hydrostatic  under  the  free  surface), 
and  Reynolds’  stresses  are  parameterized  by  an  eddy  viscosity  v, .  Only  the  time- 
varying  solution,  111(1),  is  developed  by  Putrevu  and  Svendsen,  mean  water  surface 
elevations  and  velocities  are  neglected.  The  steady  solution  could  be  linearly 
added.  The  boundary  conditions  in  z  are  given  by  Equations  4.2  and  4.3,  with  a 
slight  difference  in  the  bottom  friction  formulation  (tj,  =  0.5  The 

initial  condition.  Equation  4.4,  is  replaced  with  the  assumption  of  periodicity  in 


All  components  of  the  infragravity  wave  problem,  Vj,  f;  ,t,  ,  uj,  and 
/(see  Chapter  4  for  definitions)  are  assumed  to  have  the  form 

ijf  =  ■jjf^  (5.18) 

where  w  is  the  infragravity  wave  frequency,  ^  is  a  dummy  variable  which 
represents  any  of  the  components  listed  above,  and  is  the  component  amplitude, 
which  may  be  a  function  of  x  and  z.  The  radiation  stress  forcing  terms  and  uj 
are  in  phase  with  the  short-wave  temporal  variation.  The  surface  elevation 
gradient  term  will  not  in  general  be  in  phase  with  the  short-wave  variation.  The 
surface  elevation  includes  both  forced-  and  free-wave  components.  The  component 
phasing  is  selected  so  that/  is  maximum  at  r  =  0  and  the  surface  stress  is  defined 

where  A  is  the  phase  of  t,  relative  to  /  and  /  is  the  surface  stress  amplitude. 

The  infragravity  velocity  profile  amplitude,  U^,  is  split  into  three  parts 

and  the  time-varying  velocity  profile  is  given  by  Uj  =  ReiU„  e  The  first  part 
of  the  solution,  14;  ,  represents  the  depth  uniform  response  to  the  depth  uniform 
forcing  and  is  given  by 


This  is  a  particular  solution  to  the  governing  equation. 


The  second  and  third  components  of  the  solution,  U„2  and  ,  are 
homogeneous  solutions  which  are  subject  to  the  inhomogeneous  boundary 
conditions  at  the  bottom  and  surface  (Equations  4.2  and  4.3),  respectively.  The 
solution  for  is  given  by 

- ^  *  ♦=>  (5.22) 

pv,PV2 


where 


P  = 


(5.23) 


^  _  ^cos^(p(z+/i(^)  +  smh^(p(z+*j^) 

ysiiih^(p*(^  +  siii^(pA(^ 


(5.24) 


+  tan"* 


/ sinhp^j,  cospfe 


^coshp^Q  sinp^, 


o)J 


(5.25) 


4)2  =  tan”^ 


/  — 


sinhp(z+^j^  sinPCz  +  fto)'! 


y  coshp(z+^o)  cosP(z+Ao) 


(5.26) 


The  solution  for  is  given  by 
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and  where  ^0  water  surface  elevation,  and  Uib  and  are  determined 

by  the  solutions  for  Uai  and  17^ 
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The  shape  of  the  velocity  profiles  is  determined  by  the  relative 
magnitudes  of  the  three  velocity  terms  ^^2,  and  ,  the  depth  variation  given 
by  F2  and  Fj  ,  and  the  relative  phase  between  terms  determined  by  A. 

5.2.2  Application  of  Putrevu-Svendsen  and  eigenfunction  proHle  models  to 
surf  beat 

Putrevu  and  Svendsen  (1995)  present  three  dimensionless  parameters 
that  determine  the  relative  importance  of  the  three  terms  in  Equation  5.32.  The 
first  parameter  is^  /g,  which  is  the  forcing  that  is  uniform  over  depth.  The  depth- 
uniform  forcing  is  composed  of  gradients  in  the  infragravity  water  surface 
elevations  and  the  radiation  stress  below  the  mean  water  level  (gradient  of  the  first 
term  on  the  right-hand  side  of  Equation  4.46  or  by  Equation  4.47),  The  second 
parameter  is  F;  =  Vl  /S/(p  g),  which  is  the  surface  stress  forcing.  The  surface 
stress  forcing  is  the  shear  at  the  mean  water  level  due  to  the  radiation  stress  above 
the  mean  water  level  and  is  the  gradient  of  the  second  term  on  the  right-hand  side 
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of  Equation  4.46  or  by  Equation  4.48.  The  third  parameter  is  P2  =  Xw  ^</(^  v,  ^), 
which  is  the  bottom  stress  forcing  (Equation  2.37)  which  resists  the  flow.  In  this 
section  we  will  look  at  the  relative  importance  of  these  three  parameters  for  the 
surf  beat  Cases  2  and  3.  Also,  the  results  from  the  Putrevu-Svendsen  model  will 
be  compared  with  the  eigenfunction  model  derived  in  Chapter  4. 

Case  2 

Figure  5.19  shows  the  three  parameters  discussed  above  as  a  function 
of  X  for  surf  beat  Case  2.  The  top  figure  shows  the  relationship  of  the  two  major 
parameters,  /g  (solid  line)  and  P,  (dashed  line).  For  orientation,  the  bottom 
plot  shows  the  dimensionless  envelope  of  the  infragravity  wave,  BA  =  \^i\ /(ha). 

A  discontinuity  occurs  in  the  parameters  at  the  breaker  point.  The  surface  stress 
term  dominates  the  depth-uniform  forcing  outside  the  surf  zone  meaning  that  there 
is  significant  depth-variation  in  the  velocity  profiles.  As  the  short  waves  decay  in 
the  surf  zone,  the  surface  stress  term  loses  dominance  because  the  surface  stress  is 
decaying  in  proportion  to  H  dH/hx.  In  the  inner  surf  zone  there  is  a  crossover 
point,  and  from  there  to  the  shoreline  the  depth  uniform  forcing  dominates.  This 
is  because  the  short  waves  have  decayed  in  this  region,  and  the  infragravity  wave 
is  approaching  an  anti-node  at  the  shoreline  (the  gradient  in  the  infragravity  wave 
surface  elevation  is  the  main  depth-uniform  forcing). 

Figure  5.20  shows  the  ratio  of  the  surface  forcing  to  the  depth  uniform 
forcing  (solid  line).  The  cross-over  point  from  surface  forcing  dominance  to  the 


figure  5.19.  Case  2:  velocity  profile  parameters,  top:  solid  line  —fjg  (depth- 
uniform  forcing),  dashed  line  --  Pi  =  Vl  P/(pg)  (surface  stress  forcing), 
bottom:  solid  line  -  BA  (surf  beat  envelope),  dashed  line  -  uj(2  v,  13)  (bottom 
friction),  chain  dot  line  —  A/x  (relative  phase  of  surface  stress  and  depth-uniform 
forcing) 
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depth-unifonn  forcing  dominance  occurs  at  a  ratio  equal  to  1/  Vl,  as  seen  from 
taking  the  ratio  of  the  amplitudes  of  Equations  5.21  and  5.22. 

The  lower  plot  in  Figure  5.19  shows  the  parameter  Pj  >  which  is 
approximately  constant  through  the  surf  zone.  The  relative  phase,  A,  between^ 
and  P;  is  also  shown  in  Figure  5.19.  The  typical  value  of  A  is  approximately  t, 
meaning  that  the  depth-uniform  and  surface  stress  forcing  are  out  of  phase.  This  is 
expected  for  a  case  with  k  =  0,  since  the  infragravity  wave  would  setdown  under 
high  wave  groups  and  vice  versa.  The  relative  phase  changes  at  the  breaker  line 
and  at  the  shoreline. 

The  dashed  line  in  Figure  5.20  shows  the  variation  of  with  x.  The 
parameter  /3  (Equation  5.23)  is  analogous  to  the  parameter  used  in  the  solution  for 
wave  boundary  layers  (Putrevu  and  Svendsen  1995,  Longuet-Higgins  1953).  The 
parameter  is  inversely  related  to  the  depth  to  which  surface  and  bottom  stresses 
influence  the  velocity  profiles.  Small  values  of  indicate  surface  and  bottom 
stress  impact  the  entire  water  column,  and  large  values  indicate  the  impacts  are 
localized  near  the  boundary. 

To  compare  the  Putrevu-Svendsen  and  eigenfunction  models,  the  surf 
beat  model  described  in  Section  5.1  was  used  to  determine  the  variation  of  short- 
and  long-wave  amplitudes  and  phases  in  the  profile  models.  Only  the  time-varying 
forcing  was  applied,  so  mean  forcing  was  neglected.  The  results  of  the  Putrevu- 
Svendsen  model  are  given  in  Figure  5.21,  and  the  results  of  the  eigenfunction 
model  are  given  in  Figure  5.22.  In  the  plots,  the  velocities  have  been 


Figure  5.21.  Putrevu-Svendsen  time-varying  velocity  profiles  for  Case  2  surf  beat 
forcing  for  locations  outside  and  through  the  surf  zone 


Figure  5.22.  Eigenfunction  time-varying  velocity  profiles  for  Case  2  surf  beat 
forcing  for  locations  outside  and  through  the  surf  zone 
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nondimensionalized  by  the  group  velocity  and  the  elevation,  measured  from 
bottom,  has  been  nondimensionalized  by  the  still-water  depth. 

Both  models  produce  expected  results.  Outside  the  surf  zone,  the 
profiles  are  fairly  uniform  over  depth  with  greater  variation  near  the  surface.  The 
surface  stress  forcing  is  dominant,  but  the  effect  reaches  only  about  one  quarter  of 
the  water  colunm.  Progressing  through  the  surf  zone,  the  depth-uniform  forcing 
becomes  more  and  more  important,  and  the  profile  variation  becomes  more 
uniform  over  depth. 

The  results  firom  the  two  models  are  the  same.  The  major  difference 
between  the  models  is  that  the  Putrevu-Svendsen  model  assumes  that  all  forcing 
components  have  the  same  periodicity,  while  the  eigenfunction  model  does  not.  For 
the  Schaffer-Svendsen  surf  beat  input,  the  forcing  components  do  have  the  same 
periodicity. 

Case  3 

The  parameters  for  Case  3  are  shown  in  Figures  5.23  and  5.24.  The 
difference  between  Cases  2  and  3  is  the  long-wave  period  for  Case  3  is 
approximately  half  the  period  for  Case  2.  This  means  that  the  wave  length  for 
Case  3  is  shorter,  and  thus,  nodes  and  antinodes  for  both  the  long-wave  surface 
elevation  and  the  short-wave  groups  are  captured  in  the  surf  zone.  This  gives 
larger  gradients  in  the  long-wave  elevations  and  short-wave  groups  and  greater 
forcing.  The  short-wave  breaking  characteristics  are  similar  in  both  cases.  For 
Case  3,  the  dominance  of  the  surface  stress  term  is  more  pronounced  outside  the 


no 


Figure  5.23.  Case  3;  velocity  profile  parameters,  top:  solid  line  ~fjg  (depth- 
uniform  forcing),  dashed  line  -  P;  =  Vl  ^/(pg)  (surface  stress  forcing), 
bottom:  solid  line  —  BA  (surf  beat  envelope),  dashed  line  —  uj(2  v,^)  bottom 
friction),  chain  dot  line  —  A/x  (relative  phase  of  surface  stress  and  depth-uniform 
forcing) 
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Figure  5.24.  Relative  strength  of  surface  stress  to  depth-uniform  surf  beat  forcing 
for  Case  3  (dashed  line  —  ^h,  solid  line  --  (T^  ^)/(p  fj) 
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surf  zone  and  in  the  outer  surf  zone.  Also,  the  switch  to  dominance  of  the  depth- 
uniform  forcing  is  also  more  pronounced.  The  overall  trends  for  Case  3  are 
similar  to  Case  2. 

The  plots  of  the  velocity  profiles  (Figures  5.25  and  5.26)  show  the 
progression  from  surface  stress  dominance  to  depth-uniform  forcing  dominance 


from  outside  the  surf  zone  to  the  shoreline. 


Figure  5.25.  Putrevu-Svendsen  time-varying  velocity  profiles  for  Case  3  surf  beat 
forcing  for  locations  outside  and  through  the  surf  zone 


Figure  5.26.  Eigenfunction  time-varying  velocity  profiles  for  Case  3  surf  beat 
forcing  for  locations  outside  and  through  the  surf  zone 
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Chapter  6 

SUPERTANK  LABORATORY  DATA  COLLECTION  PROJECT 

The  SUPERTANK  Laboratory  Data  Collection  Project  was  conducted 
by  the  U.S.  Army  Engineer  Waterways  Experiment  Station,  Coastal  Engineering 
Research  Center  (CERC),  in  the  large  wave  channel  at  O.H.  Hinsdale  Wave 
Research  Laboratory,  Oregon  State  University  (OSU),  during  August  and 
September  1991  (Kraus  et  al.  1992,  Kraus  and  Smith  1994,  Smith  and  Kraus 
1995).  The  purpose  of  SUPERTANK  was  to  collect  data  unbiased  by  scaling 
distortions  to  verify  and  refine  existing  predictive  technology  and  provide  data  and 
insight  on  detailed  physical  proces^s  for  development  of  the  next  generation  of 
coastal  processes  numerical  and  physical  models.  The  specific  objectives  of 
SUPERTANK  were  to:  (a)  collect  data  to  verify  and  improve  existing  macro-scale 
beach  profile  change  numerical  simulation  models;  (b)  collect  data  to  develop 
advanced  hydrodynamic,  cross-shore  sand  transport,  and  meso-scale  beach  profile 
change  numerical  simulation  models;  (c)  collect  data  to  quantify  performance  of 
sandbars  constructed  offshore  as  a  beneficial  use  of  dredged  material;  (d)  test  and 
compare  sediment-sensing  acoustic  instruments  in  a  controlled,  field-scale 
environment  in  support  of  dredging  research;  and  (e)  collect  data  to  improve 


understanding  of  micro-scale  fluid  and  sand  motion.  SUPERTANK  was  conducted 
as  a  multi-institutional  effort,  pooling  expertise,  instrumentation,  and  a  wide  range 
of  research  interests. 

This  chapter  gives  an  overview  of  SUPERTANK,  and  then  describes 
hydrodynamic  data  collection  and  data  analysis.  SUPERTANK  data  were  used  to 
test  the  model  components  developed  in  Chapters  3  and  4.  Model  comparisons 
using  the  data  are  made  in  Chapter  7. 

6.1  Experiment  Overview 

SUPERTANK  was  conducted  in  the  large  wave  channel  at  OSU.  The 
channel  is  104  m  long,  3.7  m  wide,  and  4.6  m  deep  and  is  enclosed  in  a  building. 
For  SUPERTANK,  sand  was  trucked  from  the  Oregon  coast  and  placed  in  the 
channel.  The  channel  was  first  partially  filled  with  water  so  that  the  sand  would 
fall  through  the  water  and  become  wetted  and  somewhat  washed.  The  sand  was 
very  well-graded  with  a  median  diameter  of  0.22  mm  and  a  fall  speed  of 
3.3  cm/sec.  Some  tree  bark  and  other  organic  matter  were  contained  in  the  sand. 
The  sand  was  washed  by  running  waves  on  the  beach,  then  draining  the  tank,  and 
repeating  the  procedure.  Organic  matter  on  the  sand  surface  was  removed.  The 
beach  was  composed  of  approximately  600  m^  of  sand  that  formed  a  76-m-long 
beach.  Figure  6.1  shows  a  view  of  the  wave  channel  during  SUPERTANK. 
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Figure  6.1.  View  of  the  OSU  wave  channel  during  SUPERTANK 


Waves  were  generated  in  the  wave  channel  with  a  hydraulically-driven, 
hinged-flap  wave  board.  The  generator  has  direct  digital  control  and  is  equipped 
to  actively  absorb  reflected  waves  at  the  peak  spectral  period.  Broad-  and  narrow- 
band  random  and  monochromatic  waves  were  run  with  zero-moment  wave  heights 
in  the  range  0.2  to  1.0  m  and  with  peak  spectral  periods  in  the  range  of  3  to 
10  sec.  The  TMA  spectral  shape  (Bouws  et  al.  1985),  applicable  to  finite  water 
depths,  was  used  to  design  random-wave  tests,  with  spectral  width  parameter  7 
between  1  (broad-banded)  and  100  (narrow-banded). 

Waves  were  typically  run  in  segments  of  10-,  20-,  40-  and  70-min 
duration.  The  choice  of  these  run  lengths  was  determined  by  the  competing  factors 
of  long  runs  for  efficiency  and  statistical  confidence  versus  short  runs  to  capture 
the  beach  profile  evolution,  prevent  long-wave  buildup,  and  produce  reasonable 
data  file  sizes.  Profile  change  proceeded  rapidly  at  first  and  gradually  slowed  as 
the  profiles  assumed  an  equilibrium  shape,  therefore,  the  wave  run  intervals  in  a 
test  series  were  short  at  first  and  gradually  increased  to  record  the  change  with 
profile  surveys.  Kriebel  and  Smith  (1994)  showed  that  long-period  oscillations  did 
occur  during  runs,  but  did  not  grow  as  a  function  of  run  length.  The  magnitudes 
of  low-firequency  wave  heights  were  similar  to  those  observed  in  the  field. 

Over  a  six-week  period,  228  data  runs  were  completed  for  a  total  of 
approximately  130  hours  of  wave  action.  The  major  test  series  were  erosional 
waves,  accretional  waves,  dune  erosion,  seawall  impacts,  suspended  sediment 
transport,  berm  flooding,  and  wave  attenuation  by  offshore  mounds.  At  the  peak 


of  data  collection  activities,  the  wave  channel  was  instrumented  with  16  resistance 
wave  gauges,  10  capacitance  wave  gauges,  18  two-component  electromagnetic 
current  meters,  34  optical  backscatter  sensors,  10  pore-pressure  gauges,  3  acoustic 
sediment  concentration  profilers,  1  acoustic  Doppler  current  profiler,  1  four-ring 
acoustic  benthic  stress  gauge,  1  laser  Doppler  velocimeter,  5  video  cameras,  and  2 
underwater  video  cameras.  Synchronous  sampling  by  separate  data  acquisition 
systems  was  accomplished  by  digital  input  of  a  WWV  time  code  (broadcast  by  the 
National  Institute  of  Standards  and  Technology,  Fort  Collins,  Colorado)  to  all 
computer  clocks. 

The  beach  profile  was  surveyed  with  an  auto-tracking,  infrared 
Geodimeter  positioned  at  the  shoreward  end  of  the  wave  channel.  The  Geodimeter 
targeted  a  prism  attached  to  the  top  of  a  survey  rod  mounted  on  a  manually 
operated  carriage.  The  survey  rod,  which  moved  freely  in  a  supporting  roller 
guide,  made  contact  with  the  bed  through  a  pair  of  wide-tread  rubber  wheels,  x, 
y,  z  survey  coordinates  were  collected  automatically  on  a  personal  computer 
connected  to  the  Geodimeter.  The  nominal  horizontal  survey  spacing  was  0.3  m, 
with  finer  resolution  of  steep  features  such  as  bars  and  dunes.  The  beach  profile 
was  surveyed  along  the  center  line  of  the  channel.  Occasionally,  additional  lines 
located  0.9  m  from  each  of  the  channel  walls  were  surveyed  to  assess  uniformity 
of  the  profile  across  the  channel.  Beach  profile  data  were  visually  checked  for 
obvious  errors  (e.g.,  occasionally,  the  survey  rod  would  stick  and  the  wheels  were 


not  touching  the  bottom).  These  points  were  removed  and  replaced  through 
nearest  neighbor  interpolation. 

6.2  Hydrodynamic  Data  Collection 

At  SUPERTANK,  waves  and  currents  were  measured  along  the  wave 
channel  with  16  resistance  wave  gauges  and  18  electromagnetic  current  meters. 
Wave  gauges  were  deployed  at  3.66  m  intervals  from  the  mid  surf  zone  to  the 
wave  paddle.  The  wave  gauges  measured  the  water  surface  elevation  as  a  function 
of  time  to  give  wave  height,  wave  period,  and  mean  and  time-varying  water  level 
throughout  the  channel.  Current  meters  were  deployed  in  vertical  arrays  inside 
and  just  outside  the  surf  zone.  The  current  meters  measured  the  cross-shore  and 
vertical  water  velocities  as  a  function  of  time  to  give  mean  and  time-varying 
current  velocities. 

6.2.1  Listrument  description 
Resistance  wave  gauges 

Sixteen  resistance  wave  gauges  measured  the  elevation  of  the  water 
surface  to  determine  water  level  and  wave  parameters.  The  resistance  wave  gauges 
were  fabricated  by  the  staff  of  the  O.H.  Hinsdale  Wave  Research  Laboratory 
pibble  and  Sollitt  1989).  Each  gauge  was  2.1  m  tall  and  consisted  of  two 
0.69-mm  seven-strand  stainless  steel  cables,  spaced  6.4  cm  apart.  The  resistance 
gauges  use  the  conductivity  of  the  water  to  determine  the  water  elevation.  The 
resistance  between  the  two  wires  is  inversely  proportional  to  the  depth  of  water. 


By  monitoring  the  change  in  resistance,  the  water  surface  elevation  was 
determined.  Each  gauge  was  driven  by  a  constant  amplitude  AC  voltage  source 
and  the  current  drawn  by  the  gauge  was  measured.  By  Ohm’s  law,  with  the 
voltage  held  constant,  the  current  is  inversely  proportional  to  the  resistance  or 
directly  proportional  to  the  water  elevation.  Current  sensing  is  accomplished  by 
measuring  the  voltage  drop  across  a  1-ohm  resistor  placed  in  series  with  the  gauge. 
The  conversion  from  volts  to  engineering  units  is  determined  by  calibration  of  the 
gauges,  as  described  in  Section  6.2.5. 

Electromagnetic  current  meters 

Cross-shore  and  vertical  water  velocities  were  measured  with  18 
Marsh-McBimey  Model  512  electromagnetic  current  meters.  Each  current  meter 
sensor  consists  of  a  3.8-cm  diameter  sphere,  which  contains  an  electromagnet. 

Four  electrodes  are  equally  spaces  around  the  sphere.  The  electromagnet  generates 
a  magnetic  field,  and,  as  the  water  (a  conductor)  flows  through  the  magnetic  field, 
it  produces  a  voltage  that  is  proportional  to  its  velocity  (Faraday  principle).  The 
voltage  is  detected  with  the  electrodes.  A  rod  containing  the  instrument  cable  is 
attached  to  the  sphere,  which  in  turn  is  attached  to  a  flat  metal  mounting  bar. 

6.2.2  Instrument  location 
Resistance  wave  gauges 

Wave  gauges  were  deployed  at  3.66-m  intervals  on  the  west  wall  of  the 
wave  channel.  Gauges  were  attached  to  the  tank  wall  using  existing  2.54-cm  wall 


inserts  (threaded  holes).  Gauge  mounts  consisted  of  2.54-cm  aluminum  rods, 
threaded  on  one  end  to  attach  to  the  wall  inserts.  The  other  end  of  each  rod  was 
machined  to  accept  a  Delrin  plate,  to  which  the  wave  gauge  wires  were  attached. 
The  gauge  mounts  extended  horizontally  out  from  the  wall  approximately  0.46  m 
and  the  gauges  extended  vertically  from  z  =  -0.3  m  to  -2.4  m,  where  z  is 
measured  vertically  from  the  top  of  the  wave  tank.  The  typical  still-water  level 
was  at  z  =  -1.52  m.  The  top  mount  incorporated  an  adjustable  eyebolt  to  set  wire 
tension. 

Electromagnetic  current  meters 

Fourteen  current  meters  were  deployed  on  the  east  wall  of  the  wave 
channel  in  vertical  arrays  of  one  to  four  gauges.  The  arrangement  of  the  current 
meters  varied  from  week  to  week,  with  some  changes  from  run  to  run,  to  meet 
measurement  objectives.  Four  additional  current  meters  were  deployed  in  a 
vertical  array  (approximate  0.3-m  vertical  spacing)  off  a  mobile  carriage.  The 
position  of  the  carriage  and  elevation  of  the  meters  on  the  carriage  were  varied 
from  run  to  run. 

The  current  meter  mounts  consisted  of  aluminum  plates  and  rods.  The 
plates  were  bolted  into  the  wall  inserts  with  2.54-cm  bolts.  Spacers  welded  to  the 
back  of  the  plates  set  the  plates  approximately  5  cm  out  from  the  wall.  A  vertical 
array  of  rectangular  slots  was  machined  into  the  plates  at  16.2-cm  spacings. 
Aluminum  rods  of  1.9-cm  diameter,  threaded  and  keyed  on  one  end,  slid  into  the 
slots  in  the  plate  and  were  fastened  with  a  nut  on  the  backside  of  the  plate  (the 


space  between  the  wall  and  plate  allowed  access  behind  the  plate).  The  elevation 
of  the  rods  could  be  easily  adjusted  between  rans  by  removing  the  nut  behind  the 
plate  and  refastening  the  rod  into  another  slot.  The  opposite  end  of  the  rod  was 
marhinprf  flat  on  one  side,  and  the  mounting  bar  of  the  current  meter  was  attached 
to  the  flat  portion  of  the  rod  with  two  small  hose  clamps.  Current  meter  cables 
were  secured  to  the  plate  with  cable  ties  attached  through  small  holes  punched  into 
the  plate  for  that  purpose.  The  meters  extended  approximately  0.6  m  out  from  the 
rhannrf  wall.  Figure  6.2  shows  an  example  of  the  current  meter  mounting 
configuration.  The  current  meters  deployed  off  the  carriage  were  hose  clamped  to 
angles  bolted  to  an  aluminum  wing-shaped  section  hanging  down  from  the  carriage. 
The  wing  elevation  was  adjusted  with  a  winch  and  clamped  into  place  to  provide  a 
steady  platform  for  the  meters. 

'6.2.3  Instrument  characteristics 
Resistance  wave  gauges 

The  largest  source  of  error  in  the  resistance  wave  gauge  is  expected  to 
be  the  change  in  resistivity  (change  in  conductivity  of  the  water,  contamination  of 
the  probes,  or  residual  errors  in  the  calibration  procedure),  typically  2  percent. 
The  range  of  the  gauges  was  fixed  by  the  vertical  range  of  the  wires,  z  =  -0.3  to 
-2.4  m.  Occasionally,  waves  exceeded  the  top  of  a  gauge  (produced  by  a 
combination  of  large  wave  height  and  high  water  level)  and  saturated  the  signal. 
Gauge  calibrations  were  linear. 


125 


Electromagnetic  current  meters 

The  manufacturer’s  specifications  state  that  the  current  meters  have  an 
accuracy  of  ±2  percent.  The  range  of  the  meter  is  ±5  V,  which  translates  to 
±  5  m/sec  or  3  m/sec  (depending  on  factory  calibrations).  The  meters  were 
dq)loyed  in  groups  of  two  or  four  using  the  close  proximity  option.  This  reduced 
interference  between  closely  spaced  meters  by  driving  all  meters  in  a  group  with 
one  timing  pulse. 

6.2.4  Experiment  procedures 

The  sequence  of  a  typical  run  was: 

a.  Specify  wave  conditions  and  data  collection  parameters.  The  wave 
conditions  were  either  monochromatic  (specified  wave  height  and  period)  or 
random  (specified  target  zero-moment  wave  height  peak  period  Tp,  and 
peakedness  parameter  7  to  generate  a  TMA  spectrum).  Data  collection  parameters 
included  sampling  frequency  and  run  length. 

b.  Check  gauges  and  water  level.  Checks  were  made  of  the  gauges 
(checking  gauge  voltages  to  ensure  they  were  operating  properly)  and  the  water 
level. 

c.  Position  the  carriage.  The  carriage  with  the  current  meters  was 
positioned  at  a  strategic  location  (e.g.,  in  the  expected  incipient  breaking  zone, 
adjacent  to  wall-mounted  meters,  in  the  trough). 


d.  Begin  wave  action  and  data  collection.  One  computer  program 
controlled  both  the  wave  generation  and  the  collection  of  the  resistance  wave  gauge 
and  electromagnetic  current  meter  data.  The  data  collection  began  as  soon  as  the 
wave  generation  program  began  (although  wave  action  did  not  begin  immediately) 
and  ended  when  the  wave  action  ceased.  Typically,  wave  action  lasted  10  to 

70  min. 

e.  Lab  notes.  Notes  were  kept  documenting  significant  features  of  the 
run,  e.g.,  breaker  type  and  approximate  position,  instrument  status,  and  visual 
observation  of  sand  ripples. 

f.  beach  profile  survey.  At  the  beginning  of  the  day  and  following 
each  run,  the  beach  profile  was  surveyed  (the  post-run  survey  typically  served  as 
the  pre-run  survey  for  the  next  run).  Divers  also  measured  the  elevation  of  the 
lowest  current  meter  in  each  vertical  array  and  made  any  necessary  meter 
adjustments.  Current  meter  positions  were  surveyed  once  a  week  and  after  major 
changes  in  configuration. 

g.  Preliminary  analysis  of  data  for  quality  control.  For  at  least  one 
run  a  day  (typically  the  first  run  of  the  day),  plots  were  made  of  the  time  series 
and  spectra  for  all  gauges  to  ensure  that  the  gauges  were  operating  properly. 


6.2.5  Calibration  and  zero  reference 


Resistance  wave  gauges 

Wave  gauges  were  calibrated  once  a  week  during  the  project.  They 
were  calibrated  by  slowly  raising  and  lowering  the  water  level  (approximately 
0.3  m/hr)  and  recording  the  water  level  change  with  the  gauges.  Simultaneously,  a 
video  camera  recorded  the  water  level  relative  to  a  0.3-m-increment  scale  painted 
on  the  tank  wall.  Wave  gauge  data  were  linearly  regressed  with  the  video 
measurements  to  estimate  gauge  gains  and  offsets.  Calibration  gains  for  each  day 
of  the  experiment  were  linearly  interpolated  between  the  weekly  calibrations. 

Gauge  offsets  were  calculated  in  the  analysis  software  as  the  average  water  surface 
elevation  during  the  quiescent  period  at  the  beginning  of  a  run  (average  over  the 
first  200  points  in  random  wave  runs,  average  over  the  first  20  points  in 
monochromatic  wave  runs).  The  run-by-run  offsets  were  calculated  to  adjust  for 
small  changes  in  water  level  between  runs  (leaking  or  filling  of  the  tank)  and 
changes  in  the  sand  level  on  the  most  shoreward  gauges  (due  to  beach  profile 
change). 

Electromagnetic  current  meters 

Current  meters  were  calibrated  between  September  24  and  27,  1991  at 
the  U.S.  Geological  Survey  Indoor  Hydraulic  Laboratory  Facility  in  Bay  St.  Louis, 
Mississippi.  Although  most  of  the  meters  had  been  calibrated  within  one  year 
prior  to  the  project,  the  post-project  calibrations  were  used  for  all  analyses.  As 
with  the  wave  gauges,  the  current  meter  offsets  for  submerged  meters  were 


calculated  from  the  beginning  of  the  data  record.  For  gauges  that  were  out  of  the 
water^at  the  beginning  of  the  run,  the  calibration  offsets  were  applied. 

6.2.6  Sampling 

Analog  signals  from  the  wave  gauges  and  electromagnetic  current 
meters  were  conditioned  to  a  0-  to  10-V  DC  dynamic  range.  Each  signal  was  low- 
pass  filtered  through  a  fifth-order,  10-Hz  Bessel  filter  at  the  input  to  the  data 
acquisitidn  system  to  eliminate  noise  and  avoid  aliasing.  Analog  signals  were 
digitized  at  16  Hz,  optically  encoded,  and  transmitted  to  a  VAX  3400  server. 
Digitized  records  were  serially  recorded  in  256  time  unit  blocks  to  random  access 
memory  and  transferred  to  a  hard  drive. 

6.3  Data  Analysis 

The  analysis  routines  for  the  resistance  wave  gauge  and  electromagnetic 
cmrent  meter  data  were  nearly  the  same,  with  additional  preprocessing  required  for 
the  current  meter  data.  The  data  were  filtered  to  separate  short-wave  motions  from 
low-frequency  motions,  and  then  spectral  and  zero  up-crossing  analyses  were 
performed  on  the  raw,  low-pass,  and  high-pass  signals.  Sample  time  series, 
spectra,  and  mean  data  plots  are  presented.  Analysis  of  wave  reflection  and  wave 
groups  are  also  given. 


6.3.1  Analysis  procedure 
Preprocessing 

The  electromagnetic  current  meter  data  required  two  types  of 
processing  prior  to  analysis.  First,  a  low-pass  filter  with  a  2-Hz  cutoff  was  applied 
to  remove  high-frequency  noise.  High-frequency  noise  was  present  in  some  data 
records  for  deployments  with  the  current  meters  spaced  less  than  30  cm  apart  due 
to  interference  between  the  gauges.  Second,  the  elevation  of  each  current  meter 
was  conapared  with  the  water  surface  elevation  (estimated  from  a  resistance  wave 
gauge  at  the  same  cross-shore  location)  to  determine  if  the  current  meter  was 
sufficiently  submerged.  The  criterion  for  submergence  was  that  the  water  level 
was  5  cm  or  more  above  the  surveyed  current  meter  elevation.  Current  meter 
readings  were  set  to  zero  for  the  portion  of  the  record  when  the  gauge  was  out  of 
the  water.  The  current  meters  deployed  on  the  carriage  and  current  meters 
shoreward  of  the  most  shoreward  resistance  gauge  could  not  be  checked  for 
submergence  because  there  was  not  a  co-located  resistance  wave  gauge.  These 
gauge  outputs  were  analyzed  assuming  100  percent  submergence;  thus,  the  results 
are  viewed  with  healthy  skepticism. 

Filtering 

To  separate  short-wave  motions  from  low-frequency  motions,  a  non¬ 
recursive,  low-pass  filter  was  applied.  The  period  cutoff  for  the  filter  was  set  to 
twice  the  peak  period  of  the  incident  short  waves  (the  peak  period  of  the  long- 
period  motion  was  generally  20  sec  or  longer  as  compared  to  typically  3  to  5  sec 


for  the  short  waves).  The  period  response  of  the  filter  is  shown  in  Figure  6.3  for 
a  cutoff  period  of  6  sec.  The  filter  was  applied  beginning  with  the  1009th  data 
point.  The  first  1008  points  were  skipped  to  ensure  that  the  waves  had  reached  the 
most  shoreward  gauge.  Two  hundred  points  were  sacrificed  at  the  beginning  (in 
addition  to  the  1008  points  that  were  skipped)  and  end  of  the  data  record  in 
constructing  the  filter.  The  low-pass  time  series  was  output  directly  from  the 
filter.  The  high-pass  time  series  was  obtained  by  subtracting  the  low-pass  time 
series  from  the  preprocessed  data. 

Spectral  analysis 

Spectral  analysis  was  performed  using  a  standard  Fast  Fourier 
Transform  (Press  et  al.  1989).  Data  were  divided  into  blocks  of  4096  points  and 
the  spectrum  was  calculated  for  each  block.  The  total  number  of  blocks  per  run 
was  determined  by  the  run  length  (20-min  run  =  4  blocks,  40-min  run  =  8  blocks, 
and  70-min  run  =  15  blocks).  The  average  spectrum  for  each  run  was  calculated 
by  averaging  the  spectra  over  all  blocks  within  the  run.  The  confidence  bands  for 
each  spectrum  become  tighter  as  the  number  of  averages  increases.  The  degrees  of 
fi-eedom  (DOF)  for  the  spectra  are  approximately  twice  the  number  of  averages 
(20-min  run  =  8  DOF,  40-min  run  =  16  DOF,  and  70-min  run  =  30  DOF). 

Mean  water  elevation,  standard  deviation,  root-mean-square  wave 
height,  and  peak  period  were  determined  from  the  calculated  spectra  for  the 
wave  gauges.  Similarly,  mean  velocities  and  standard  deviations  were  estimated 
from  the  calculated  spectra  for  the  current  meters. 


Figure  6.3.  Period  response  of  the  low-pass  filter  for  cutoff  period  6  sec 


Time  series  analysis 


Time  series  analysis  was  performed  with  a  zero-upcrossing  definition  of 
a  wave.  A  wave  is  defined  between  two  successive  upward  crossings  of  the  water 
surface  about  the  mean  water  elevation.  The  upcrossing  definition  was  chosen 
because  the  saw-tooth  shape  of  the  nearshore  wave  profile  gives  distinct  upward 
crossings,  whereas  downcrossings  vary  from  wave  to  wave.  The  wave  height  H  of 
an  individual  wave  is  defined  as  the  difference  between  the  highest  and  lowest 
water  elevation  between  two  zero  upcrossings.  The  wave  period  T  of  an  individual 
wave  is  defined  as  the  time  between  two  zero  upcrossings.  Very  small  waves  (J? 
<0.6  cm),  associated  with  small  surface  variations  near  the  mean  water  elevation, 
were  not  identified  as  individual  waves.  These  very  small  waves  were  included  as 
part  of  the  preceding  wave. 

The  mean,  standard  deviation,  skewness,  and  kurtosis  of  the  water 
surface  and  the  mean,  root-mean-square,  significant  (H,),  one-tenth,  and  maximum 
(fl^  wave  height  and  period  were  calculated  from  the  wave  gauges.  Mean  and 
maximum  velocities,  standard  deviation,  skewness,  kurtosis,  and  percent  of 
measurements  with  the  gauge  submerged  were  calculated  for  the  current  meters. 

6.3.2  Long-wave  motion  in  two  SUPERTANK  runs 

Data  for  SUPERTANK  Runs  A2007B  and  A0509A  are  analyzed  with 
particular  emphasis  on  determining  the  long-wave  motion.  These  runs  will  be 
studied  in  greater  detail  in  Chapter  7.  The  data  provide  insight  into  the 
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SUPERTANK  results  as  well  as  background  for  the  modeling  results  discussed  in 
Chapter  7.  The  beach  profiles  for  Runs  A2007B  and  A0509A  are  shown  in 
Figures  6.4  and  6.5,  respectively 
Water  surface  elevations 

A  sample  of  the  raw  resistance  wave  gauge  data  is  shown  in  Figure  6.6 
for  Gauge  3,  Run  A2007B.  The  vertical  axis  is  the  water  surface  elevation  f  and 
the  horizontal  axis  is  time.  Every  fourth  data  point  of  the  time  series  is  plotted. 
The  low-pass  signal  is  also  included  in  Figure  6.6.  The  target  wave  conditions  for 
Run  A2(X)7B  were  7],  =  5.0  sec  and  =  0.7  m.  The  wave  spectrum  was  very 
narrow  with  7  =  100.  Wave  groups  are  evident  in  the  time  series,  as  would  be 
expected  with  the  relatively  narrow  incident  spectrum.  Gauge  3  is  located  in 
approximately  0.7  m  water  depth  on  the  shoreward  side  of  a  bar,  within  the  surf 
zone  (Figure  6.4).  The  mean  water  surface  elevation  for  the  record  was  0.02  m. 

The  smoothed  spectra  for  the  raw,  low-pass,  and  high-pass  signals  for 
Gauge  3,  Run  A2007B,  are  given  in  Figure  6.7.  For  this  example,  the  spectral 
peak  is  located  at  0.2  Hz  (5  sec)  and  higher  harmonics  are  apparent  at  0.4,  0.6, 
and  0.8  Hz.  Both  the  wave  kurtosis  and  skewness  increase  in  shallower  depths,  as 
harmonics  grow  and  waves  become  asymmetric.  A  low-frequency  peak  at 
0.023  Hz  (43  sec)  is  also  apparent.  At  Gauge  3,  91  percent  of  the  wave  energy  is 
contained  in  the  high-pass  signal  and  9  percent  is  contained  in  the  low-pass  signal. 
The  energy  in  the  low-pass  signal  tends  to  increase  close  to  shore. 


figure  6.4.  Summary  of  wave  transformation  and  mean  cross-shore  currents  for 
Run  A2007B 
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Figure  6.7,  Wave  spectra  for  the  raw,  low-pass,  and  high-pass  water  surface 
elevation  for  Gauge  3,  Run  A2007B 


Velocities 


Figure  6.8  shows  the  spectrum  for  the  raw  signal  and  the  raw  and  low- 
pass  time  series  for  Channel  30  (cross-shore  velocity  component,  located 
approximately  0.3  m  below  the  still-water  surface)  for  Run  A2007B.  This  current 
meter  is  located  at  the  same  cross-shore  location  as  wave  Gauge  3  shown  in 
Figures  6.6  and  6.7.  The  general  appearance  of  the  raw  time  series  is  similar  to 
the  wave  time  series  in  Figure  6.6.  The  positive  velocities  are  onshore  directed 
(corresponding  to  the  passage  of  a  wave  crest)  and  the  negative  velocities  are 
offshore  (corresponding  to  the  passage  of  a  wave  trough).  The  mean  velocity  at 
this  gauge  was  -0.12  m/sec.  The  energy  spectrum  of  the  velocity  is  similar  to  the 
wave  spectrum  in  Figure  6.7.  The  primary  peak  is  at  0.2  Hz  with  harmonics  at 
0.4  and  0.6  Hz.  Approximately  28  percent  of  the  energy  is  contained  in  the  low- 
pass  signal  and  the  low-ffequency  peak  is  at  0.02  Hz  (50  sec). 

Figure  6.4  and  6.5  give  example  summaries  of  the  wave  and  current 
data  for  Runs  A2007B  and  A0509A,  respectively.  The  incident  wave  conditions 
for  Run  A0509A  were  Tj,  =  3  sec,  =  0.8  m,  and  7  =  20.  Run  A0509A  was 
the  first  SUPERTANK  run,  and  it  started  from  an  equilibrium  profile  shape.  The 
beach  profile  for  Run  A2007B  was  formed  from  a  week  of  erosional  wave 
conditions  (generating  multiple  bars  moving  offshore)  followed  by  a  week  of 
accretional  wave  conditions  (in  which  the  bars  moved  onshore).  The  top  of  the 
figures  shows  the  wave  transformation  as  measured  by  the  resistance  wave  gauges 
and  calculated  from  the  zero-upcrossing  analysis.  is  given  by  the  dashed  line 
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Figure  6.8.  Current  meter  raw  spectrum  and  time  series  of  the  raw  and  low-pass 
signal  for  Channel  30,  Run  A2007B 


and  is  given  by  the  solid  line.  Note  that  the  wave  height  scale  is  on  the  right- 
hand  side  of  the  plot.  For  Run  A2007B,  waves  are  breaking  on  the  outer  bar  and 
decaying  as  they  travel  shoreward.  For  Run  A0509A,  the  incipient  wave  breaking 
zone  was  located  at  the  two  most  shoreward  wave  gauges.  Figure  6.9  shows  the 
decay  of  the  high-pass  wave  height  (solid  line)  and  the  growth  of  the  low-pass 
wave  height  (dashed  line)  for  Run  A0509A. 

The  mean  cross-shore  current  velocities  are  shown  in  Figures  6.4  and 
6.5  as  vectors.  The  plots  illustrate  the  tendency  for  offshore  flow  below  trough 
level  and  onshore  flow  above  trough  level.  The  maximum  mean  cross-shore 
velocities  at  SUPERTANK  were  0.5  m/sec,  shoreward  of  the  bar.  Velocities  in 
the  offshore  region  of  the  channel  were  typically  less  than  0.05  m/sec,  in  the  range 
0.2  to  0.5  (an  estimate  of  the  flow  required  to  balance  mass 

transport).  The  velocity  profiles  in  the  offshore  region  are  linear  in  shape 
(constant  over  depth  or  slightly  increasing  offshore-directed  velocity  from  the 
bottom  to  the  surface).  In  the  nearshore  region,  the  mean  velocity  profile  tended 
to  be  parabolic  in  shape,  with  greatest  velocities  near  the  bottom  and  lower 
velocities  near  trough  level.  The  SUPERTANK  cross-shore  flow  measurements 
show  trends  similar  to  field  measurements  at  the  DELILAH  experiment  (Smith  et 
al.  1992),  conducted  at  the  CERC  Field  Research  Facility  in  North  Carolina.  The 
velocity  profile  shape  offshore  of  breaking  and  across  the  bar  crest  are  similar  in 
the  laboratory  and  field  measurements.  However,  measured  undertow  velocities  in 
the  trough  at  DELILAH  are  generally  smaller  than  those  measured  at 
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Figure  6.9.  High-pass  and  low-pass  significant  wave  heights  for  Run  A0509A 
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SUPERTANK.  This  may  be  because  the  trough  was  deeper  and  wider  and  the 
foreshore  steeper  at  DELILAH  than  at  SUPERTANK  or  due  to  three-dimensional 
circulation  at  DELILAH. 


6.3.3  Short-  and  long- wave  reflection 

Wave  reflection  from  the  nearshore  (beach  and  bars)  was  calculated  for 
selected  wave  runs  at  several  locations  in  the  channel.  Reflection  estimates  were 
calculated  with  a  frequency-domain  linear  theory  analysis  technique  using  synoptic 
time  series  of  cross-shore  velocity  and  water  surface  elevation  (Hughes  1993). 

This  method  assumes  normally-incident  waves. 

Hughes  represents  the  water  surface  elevation  and  cross-shore  velocity 
as  the  sum  of  incident  and  reflected  irregular  linear  waves 
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N  is  the  number  of  frequencies,  k  is  the  wavenumber,  a  is  the  angular  wave 
frequency,  e  is  the  phase  angle,  a  is  the  wave  amplitude,  and  the  subscripts  I  and  R 
denote  incident  and  reflected  wave  components.  Using  trigonometric  identities. 
Equations  6.1  and  6.2  can  be  rewritten 
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Fourier  analysis  is  used  to  calculate  A,  B,  C,  and  D  from  the  time  series.  Then, 
Equations  6.8-11  are  solved  to  calculate  the  incident  and  reflected  wave  amplitudes 
and  phases 
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The  reflection  coefficient  is  the  ratio  of  the  reflected  wave  amplitude  to  the 
incident  wave  amplitude.  The  incident  and  reflected  wave  time  series  are  given  by 
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An  example  of  the  reflection  analysis  output  is  shown  in  Figure  6. 10. 
The  figure  shows  time  series  of  the  total,  incident,  and  reflected  long-wave  water 
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Figure  6.10.  Total,  incident,  and  reflected  low-pass  water  surface  elevations  for 
Run  A2(X)7B  (total  —  solid  line,  incident  --  chain-dot  line,  reflected  -  dashed  line) 


surface  elevation  for  Run  A2007B.  Separation  of  incident  and  reflected  long 
waves  is  required  to  drive  the  depth-integrated  model  discussed  in  Chapter  3. 

Reflection  coefficients  were  generally  0,20  to  0.45  (all  frequencies)  for 
the  SUPERTANK  runs.  Coefficients  for  the  low-pass  frequencies  were  in  the 
range  0.4  to  1.0.  The  reflection  coefficient  was  strongly  related  to  the  incident 
short-wave  peak  period,  with  greater  reflection  for  longer  peak  periods.  The  wave 
generator  used  at  the  wave  channel  was  equipped  to  absorb  reflected  waves  at  the 
peak  frequency,  so  re-reflection  of  wave  energy  from  the  wavemaker  was 
suppressed.  Figure  6.11  shows  total  wave  height  calculated  in  6-min  segments  for 
the  high-pass  and  low-pass  signals  over  a  70-min  run.  The  solid  lines  are  surf 
zone  wave  gauges  (depth  less  than  1  m),  the  chain-dot  lines  are  gauges  in  the 
deepest  portion  of  the  channel  (depth  approximately  3  m),  and  the  dotted  and 
dashed  lines  are  at  intermediate  depths.  There  is  a  variation  in  long-wave  height 
with  time,  but  there  is  not  a  consistent  increase  in  height  with  time  from  re¬ 
reflection  of  waves  from  the  wavemaker,  as  might  be  expected. 

6.3.4  Wave  groups 

The  incident  wave  spectra  at  SUPERTANK  included  narrow-banded 
q)ectral  shapes  (TMA  spectra  with  7  =  20  and  100).  The  narrow  spectra 
produced  wave  groups  which  induced  time  variations  in  the  low-pass  surface 
elevations  and  velocities.  The  time  series  of  high-pass  wave  energy  (£)  and  low- 
pass  surface  elevation  (f)  for  Run  A2007B  are  given  in  Figure  6. 12  for  three 


Figure  6.11.  High-pass  and  low-pass  wave  heights  as  a  function  of  time  (surf 
zone  —  solid  line,  deep  —  chain-dot  lines,  intermediate  depth  —  dashed  and  dotted 
lines) 
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Figure  6,12.  Incident  wave  energy  (E  —  dotted  line)  and  low-pass  surface 
elevation  (f-  solid  line)  offshore  (x  ==  0  m,  top),  seaward  of  incipient  breaking  (x 
=  40  m,  middle),  and  surf  zone  (x  =  51  m,  bottom) 
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locations  in  the  tank  (x  =  0,  40,  and  51  m  in  Figure  6.4).  The  high-pass  wave 
energy  was  calculated  by  squaring  the  high-pass  surface  elevation  and  then  low- 
pass  filtering  the  squared  elevation  to  get  the  envelope  of  the  energy.  The  top 
panel  of  Figure  6.12  is  a  location  near  the  wave  generator,  the  middle  panel  is  a 
location  just  seaward  of  the  incipient  breaking  zone,  and  the  bottom  panel  is  in  the 
mid  surf  zone.  The  middle  panel  shows  strong  correlation  between  the  high-pass 
wave  energy  and  the  low-pass  surface  elevation.  The  bottom  panel,  however, 
shows  that  in  the  mid  surf  zone,  most  of  the  short-wave  energy  has  been  dissipated 
due  to  wave  breaking.  The  long-wave  surface  elevation  maintains  its  structure  and 
represents  a  system  of  free  waves  propagating  towards  the  shoreline  superimposed 
on  similar  outgoing  waves  which  are  the  long-wave  motion  reflected  from  the 
shore. 

As  shown  in  the  middle  panel  of  Figure  6.12,  the  strongest  correlation 
between  the  short-wave  energy  and  the  long-wave  water  surface  elevation  occurs 
just  seaward  of  the  breaker  zone.  The  correlation  of  these  two  signals  is  shown  in 
Figure  6.13  as  a  function  of  lag  time.  A  strong  negative  correlation  occurs  at  a 
time  lag  of  approximately  zero  (short-wave  energy  leads  slightly  in  time),  as  would 
be  expected  if  the  long-wave  motion  was  generated  by  radiation  stress  forcing 
(Longuet-Higgins  and  Stewart  1964).  The  negative  correlation  means  that  high 
wave  groups  force  a  setdown  of  the  low-pass  water  surface  elevation,  and  low 
wave  groups  force  a  setup. 
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Figure  6.13.  Correlation  of  the  incident  wave  energy  and  the  low-pass  surface 
elevation  at  a  location  just  seaward  of  the  incipient  breaker  zone 
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Figure  6.14  shows  a  repeat  of  the  time  series  of  short-wave  energy  (£) 
and  long-wave  water  surface  elevation  (f)  from  the  center  panel  of  Figure  6. 12 
with  the  low-pass  velocity  measured  with  a  co-located  current  meter  (just  outside 
the  incipient  breaker  zone).  Offshore  flow  is  correlated  with  the  high  wave  groups 
and  depressions  in  the  low-pass  surface  elevation.  The  time  series  show  there  is  a 
mean  setdown  in  the  water  surface  in  this  location  and  a  mean  offshore  flow  at  this 
mid-depth  current  measurement  location.  Figure  6.15  shows  the  correlation 
between  the  long-wave  water  surface  elevation  and  the  long-wave  cross-shore 
velocity  corresponding  to  Figure  6.14.  The  positive  correlation  is  fairly  strong, 
with  the  water  surface  elevation  lagging  the  velocity  by  a  few  seconds.  The 
similar  correlation  of  the  long-wave  water  surface  elevation  and  velocity  in  the  mid 
surf  zone  is  shown  in  Figure  6.16.  In  the  surf  zone,  the  correlation  is  negative, 
and  the  water  surface  leads  the  cross-shore  velocity  by  approximately  one  fourth  of 
the  low-pass  wave  period. 

The  very  narrow  spectra  generated  at  SUPERTANK  correspond  to  time 
series  with  significant  wave  groups.  The  wave  groups  generate  long  waves  similar 
to  the  surf  beat  model  of  Schaffer  and  Svendsen  discussed  in  Chapter  5. 


Figure  6.15.  Correlation  of  the  long-wave  surface  elevation  and  the  long-wave 
velocity  at  a  location  just  seaward  of  the  incipient  breaker  zone 


Figure  6.16.  Correlation  of  the  long-wave  surface  elevation  and  the  long-wave 
velocity  at  a  location  in  the  mid  surf  zone 
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Chapter  7 

MODELING  TIME-  AND  DEPTH-VARYING  CURRENTS  AT  SUPERTANK 

Chapters  3  and  4  described  the  depth-integrated  and  profile  models, 
respectively,  and  Chapter  6  described  the  SUPERTANK  Laboratory  Data 
Collection  Project  data  set.  This  chapter  puUs  together  both  the  modeling 
components  and  the  laboratory  data.  Measurements  made  at  SUPERTANK  are 
used  to  drive  the  numerical  models  for  predicting  time-varying  velocity  profiles 
and  to  evaluate  model  results.  To  drive  the  numerical  models,  SUPERTANK 
measurements  provided  the  incident  infragravity  wave  time  series  at  the  offshore 
model  boundary,  time  series  of  radiation  stress  throughout  the  grid  domain,  and  the 
still-water  depth.  Low-pass  measurements  of  velocity  and  water  surface  elevation 
were  used  to  evaluate  the  models. 

7.1  SUPERTANK  Runs  A0509A  and  A2007B 

SUPERTANK  Runs  A0509A  and  A2007B  were  selected  for  simulation 
with  the  models.  The  bathymetry  for  Run  A0509A  has  an  equilibrium  shape,  as 
shown  in  Figure  7.1.  The  positions  of  current  meters  are  shown  with  triangles  in 
the  figure.  The  short-wave  parameters  for  Run  A0509A  were  zero-moment  wave 
height  of  0.8  m,  peak  period  of  3  sec,  and  a  spectral  peakedness  parameter  of  20, 


Figure  7.2.  Bathymetry  and  current  meter  positions  for  SUPERTANK 


Run  A2007B 
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a  relatively  narrow  spectrum.  The  incident  wave  conditions  are  summarized  in 
Table- 7.1  The  significant  infragravity  wave  period  was  approximately  26  sec. 

This  run  was  selected  because  it  includes  measurements  both  outside  and  inside  the 
surf  zone  and  illustrates  several  typical  features  of  model  results.  Also,  this  is  the 
only  SUPERTANK  run  with  a  simple,  non-barred  bathymetry. 

The  bathymetry  for  Run  A2(X)7B  has  a  complex  barred  shape,  as  shown 
in  Figure  7.2.  The  multiple-bar  shape  was  formed  under  a  series  of  accretionary 
wave  runs  (low  steepness  waves).  The  positions  of  current  meters  are  shown  with 
triangles  in  the  figure.  The  short-wave  parameters  for  Run  A2007B  were  zero- 
moment  wave  height  of  0.7  m,  peak  period  of  5  sec,  and  a  spectral  peakedness 
parameter  of  100,  a  very  narrow  spectrum.  The  significant  inffagravity  wave 
period  was  approximately  37  sec.  This  run  was  selected  because  of  the  very 
narrow  spectrum  and  the  resulting  strong  infragravity  modulation.  The  current 
measurements  were  made  just  outside  the  surf  zone  and  in  the  surf  zone. 

The  following  sections  describe  the  application  of  the  depth-integrated 
and  profile  models  to  SUPERTANK  Runs  A0509A  and  A2007B  and  evaluation  of 
model  results. 
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Table  7.1.  Wave  and  grid  parameters  for  Runs  A0509A  and  A2007B 

Run 

Tp,  sec 

7 

m 

ny 

Ax,  m 

At,  sec 

A0509A 

0.8 

3.0 

20. 

57 

5 

0.91 

0.125 

A2007B 

0.7 

5.0 

100. 

37 

5 

0.61 

0.125 

7.2  Depth-Integrated  Model 

The  two-dimensional-horizontal,  dq)th-integrated  Predictor-Corrector 
model  described  in  Chapter  3  is  used  to  calculate  the  long-wave  water  surface 
elevations  for  SUPERTANK.  The  governing  equations  are  continuity  and 
momentum  (Equations  3.1-3).  The  offshore  boundary  condition  is  an  incident 
long-wave  time  series,  and  the  boundary  condition  at  the  shoreline  is  full  reflection 
of  the  long  wave.  Radiation  stress  forcing,  calculated  from  the  SUPERTANK 
short-wave  measurements,  is  applied  throughout  the  nearshore.  The  depth- 
integrated  model  provides  the  long-wave  water  surface  elevation  throughout  the 
nearshore  to  drive  the  time-varying  velocity  profile  model  (Section  7.3). 

7.2.1  Model  input  and  set  up 

As  input,  the  depth-integrated  Predictor-Corrector  model  requires  the 
incident  infragravity  wave  surface  elevation  time  series  at  the  offshore  boundary. 


157 


the  time  series  of  radiation  stress  throughout  the  domain,  bottom  friction 
information,  initial  conditions,  and  specification  of  the  grid  domain. 

Offshore  boundary  condition 

As  described  in  Section  6.3.1,  the  SUPERTANK  water  surface  and 
current  measurements  were  filtered  to  sq)arate  short-wave  motions  from 
infragravity-wave  motions  using  a  low-pass  filter.  The  low-passed  signal  includes 
wave-group  forced  long  waves  and  free  long  waves  that  are  reflected  from  the 
beach  and  wave  generator  (the  generator  absorbed  reflected  waves  only  at  the 
incident  short-wave  peak  period).  The  short-wave  measurements  are  used  to 
ralmlatfi  radiation  stress,  and  the  infragravity  wave  measurements  are  used  to  drive 
the  model  at  the  offshore  boundary  and  to  evaluate  model  output  within  the 
domain. 

Incident  infragravity-wave  time  series  at  the  offshore  boundary  were 
calculated  using  a  linear  technique  to  separate  the  incident  and  reflected  waves 
using  co-located  measurements  of  cross-shore  velocity  and  water  surface  elevation 

I 

(Hughes  1993,  Section  6.3.3).  This  required  that  the  offshore  model  boundary  be 
located  at  a  position  with  co-located  wave  and  current  gauges.  The  offshore 
boundary  for  Run  A0509A  is  at  =  0.0  m  in  Figure  7.1,  in  a  water  depth  of 
2.42  m.  The  offshore  boundary  for  Run  A2007B  is  at  x  =  32.9  m  in  Figure  7.2, 
in  a  water  depth  of  0.78  m.  This  boundary  location  was  selected  for  Run  A2007B 
because  it  was  the  deepest  location  with  a  current  meter.  The  separated  incident 
and  reflected  water  surface  elevations  for  Run  A2007B  are  given  in  Figure  6. 10. 
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Note  that  both  forced  and  free  infragravity  waves  existing  in  the  tank  are  included 
in  the  offshore  water  surface  boundary  condition.  The  offshore  boundary  time 
series  was  applied  uniformly  across  the  offshore  boundary  (in  the  longshore 
direction). 

Radiation  stress  forcing 

The  radiation  stress  forcing  throughout  the  domain  is  the  most  difficult 
input  to  supply.  Using  the  SUPERTANK  water  surface  elevation  measurements, 
radiation  stress  is  calculated  from  the  measured  envelope  of  short-wave  energy 
using  linear  wave  theory  (Equation  3.31),  and  assuming  shallow-water  conditions 
for  the  short-wave  motion.  The  envelope  of  the  short-wave  energy  (rj^  was 
calculated  by  squaring  the  measured  short-wave  water  surface  elevation  (tj),  low- 
pass  filtering  using  the  same  filter  used  to  separate  the  short-wave  and  infragravity 
wave  motions,  multiplying  by  two,  and  applying  Equation  3.31,  which  results  in 

=  I  P«  7«)  <’•« 

The  spacing  of  the  SUPERTANK  water  surface  measurements  was 
3.66  m,  which  is  coarser  than  the  grid  spacing  used  in  model  applications.  To  fill 
the  spatial  gaps,  the  short-wave  energy  envelope  was  linearly  interpolated  between 
gauges.  The  most  shoreward  wave  gauge  was  located  at  x  =  47.59  for  Run 
A0509A  and  at  x  =  43.9  m  for  Run  A2(X)7B,  so  the  wave  energy  envelope  was 
linearly  interpolated  from  the  measured  value  at  the  most  shoreward  gauge  to  zero 


at  the  shoreline.  The  gauge  spacing  does  not  resolve  the  incipient  breaker  region 
well.  I  Numerical  or  analytical  models  of  the  transformation  and  breaking  of  the 
short-wave  groups  could  be  used  to  provide  the  short-wave  energy  based  on  the 
energy  equation.  Additional  empirical  assumptions  would  be  required  to  specify 
wave  breaking  for  the  wave  groups.  The  approach  taken  here  is  to  use  the 
laboratory  measurements.  Thus,  the  short-wave  breaker  region  is  not  well 
represented  due  to  a  lack  of  resolution  in  the  measurements.  The  radiation  stress 
forcing  was  applied  uniformly  across  the  grid  in  the  longshore  direction.  The 
envelope  of  the  short-wave  energy  is  used  in  Equation  7. 1  to  calculate  the  time- 
varying  radiation  stress. 

Bottom  firiction 

Bottom  friction  is  specified  using  Equation  2.37.  A  constant  bottom 
friction  coefficient  of^  =  0.01  is  used  for  all  simulations.  The  model  results 
were  insensitive  to  the  choice  of  in  the  typical  range  of  values  between  0.01- 
0.05.  The  short-wave  velocity  amplitude,  u„,  is  calculated  using  linear  theory 


_  a  a 

h  srnhm 


(7.2) 


where  a  is  the  short-wave  amplitude,  a  is  the  short-wave  radial  frequency,  and  k  is 
the  short-wave  wave  number.  To  simplify  the  calculation  of  the  velocity 
amplitude,  the  local  wave  amplitude  was  assumed  equal  to  half  the  incident 
outside  the  surf  zone.  To  approximate  breaking,  a  limiting  amplitude  of  0.4  ho 


was  applied.  Progressing  from  the  offshore  boundary  to  the  shore,  amplitudes 
clipped  by  0.4  hg  were  passed  shoreward  as  the  base  amplitude  in  place  of 
0.5  The  short-wave  amplitude  was  assumed  constant  in  time  for  the 
calculation  of  bottom  friction. 

Initial  conditions 

The  initial  conditions  for  the  depth-integrated  model  (Equation  3.4-12) 
were  specified  as  g,  =  0  m^/sec  and  f  =  0  m  throughout  the  domain.  The 
infragravity  wave  forcing  at  the  offshore  boundary  and  the  radiation  stress  forcing 
throughout  the  domain  were  also  initiated  at  zero  and  linearly  increased  to  their  full 
measured  values  over  the  first  500  time  steps.  This  ramping  of  the  forcing  reduces 
artificial  oscillations  generated  by  the  initial  conditions. 

Grid  domain 

For  Run  A0509A,  the  region  modeled  extends  from  a  depth  of  2.42  m 
to  0.20  m.  The  cross-shore  and  longshore  grid  spacing  was  0.91  m,  and  the  time 
step  used  was  0.125  sec.  With  Ax  =  0.91  m,  there  were  57  grid  cells  in  the 
cross-shore  (5  cells  in  the  longshore).  For  Run  A2007B,  the  region  modeled 
extends  from  a  depth  of  0.78  m  to  0. 12  m.  The  cross-shore  and  longshore  grid 
spacing  was  0.61  m,  resulting  in  37  cells  in  the  cross-shore  and  5  cells  in  the 
longshore.  For  both  SUPERTANK  runs,  the  time  step  used  was  0.125  sec.  The 
grid  spacings  were  selected  so  measurement  points  aligned  with  grid  cells  and  the 
Courant  criterion  of  Ar  <  Ax/C  was  met.  The  time  step  was  selected  to  equal 
twice  the  measurement  sampling  rate.  The  grid  information  for  each  run  is 


summarized  in  Table  7.1.  Water  depths  were  specified  from  the  SUPERTANK 
bathymetry  surveys. 

7.2.2  Model  results 

Results  for  the  infragravity  wave  water  surface  elevation  and  depth- 
integrated  velocities  (U  =  QJh)  are  given  in  Figures  7. 3-7. 8  for  Run  A0509A  and 
Figures  7.9-7.14  for  Run  A2007B.  The  figures  show  time  series  results  for  four  x 
locations.  The  solid  lines  are  the  low-pass  filtered  wave  and  current  measurements 
and  the  dashed  lines  are  the  model  results.  Infragravity  surface  elevations  are 
relative  to  the  still-water  level,  with  negative  values  below  still-water  level  and 
positive  values  above  still-water  level.  Positive  velocities  are  directed  onshore  and 
negative  velocities  are  offshore.  The  measured  velocities  are  determined  as  the 
average  of  the  long-wave  velocities  over  all  gauges  below  the  still-water  level,  not 
taldng  into  account  that  the  sensors  are  unevenly  distributed  over  depth.  Error 
statistics  and  mean  values  at  all  measurement  locations  are  given  in  Tables  7.2  and 
7.3  for  Runs  A0509A  and  A2007B,  respectively.  The  dimensionless  mean 
measured  and  calculated  infragravity-wave  surface  elevations  are  given  by 


Figure  7.3.  Infragravity  surface  elevations  (solid  --  measured, 
calculated)  for  Run  A0509A,  t  =  0-250  sec 


dashed  - 


Figure  7.4.  Depth-integrated  velocities  (solid  line  --  measured,  dashed  line 
calculated)  for  Run  A0509A,  t  =  0-250  sec 


163 


calculated)  for  Run  A0509A,  t  =  250-500  sec 


figure  7.6.  Dq)th-integrated  velocities  (solid  line  --  measured,  dashed  line  - 


calculated)  for  Run  A0509A,  t  =  250-500  sec 


164 


Figure  7.7.  Infragravity  surface  elevations  (solid  --  measured,  dashed 
calculated)  for  Run  A0509A,  t  =  500-750  sec 
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calculated)  for  Run  A2007B,  t  =  250-500  sec 
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Figure  7.13.  Infragravity  surface  elevations  (solid  --  measured,  dashed 
calculated)  for  Run  A2007B,  t  =  500-750  sec 


Figure  7.14.  Dq)th-integrated  velocities  (solid  line  —  measured,  dashed  line 
calculated)  for  Run  A2007B,  t  =  500-750  sec 
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Table  7.2.  Depth-integrated  model  means  and  rms  errors  for  Run  A0509A 


% 

-0.94 

0.14 

3.75 

-1.08 

-2.13 

1.15 

-0.98 

-2.01 

4.64 

-1.05 

-1.99 

4.66 

-1.06 

-1.89 

4.65 

-1.02 

-1.99 

4.72 

-1.02 

-1.93 

4.92 

-1.47 

-1.66 

5.14 

-1.30 

-1.73 

5.16 

-1.34 

-1.49 

5.24 

^  ^V.rrnn  ^ 


0.0 


3.66 


7.32 


10.97 


14.63 


18.29 


21.95 


25.60 


29.26 


32.92 


34.75 


36.58 


40.23 


43.89 


47.59 


-0.30  0.02 


5.14  -0.43  0.00 


-1.75 

-0.01 

-1.30 

-0.01 

-1.20 

-0.01 

-2.35 

-0.003 

Table  I 

f.3.  Depth-integrated  model  means  and  rms  errors  for  Run  A2007B 

jc,  m 

1^,  % 

1^,  % 

% 

u^,  % 

^U,rms9  % 

32.92 

-1.60 

-0.43 

4.41 

-1.80 

-0.05 

2.98 

36.58 

3.26 

7.89 

7.56 

-5.00 

-0.04 

6.44 

37.19 

— 

— 

— 

-2.80 

-0.04 

4.31 

40.23 

0.57 

6.76 

7.79 

-3.72 

-0.05 

5.94 

43.89 

2.82 

20.68 

21.50 

-1.41 

-0.03 

5.38 
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where  is  the  depth-averaged  measured  velocity,  is  the  velocity  calc^ated 
using  the  depth-integrated  model,  N  is  the  number  of  measurements  in  the  time 
series,  and  is  the  local  rms  wave  height.  The  dimensionless  rms  error  in 
long-wave  water  surface  elevation  is  defined 


i=l 


X  100 
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The  dimensionless  mean  measured  and  calculated  depth-integrated  velocities  are 
given  by 
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N 

1^1 

N  \/gK 


X  100 


(7.6) 


X  100 


(7.7) 


where  U„  is  the  depth-averaged  measured  velocity  and  is  the  velocity  calculated 
using  the  depth-integrated  model.  The  dimensionless  rms  error  in  velocity  is 
defined 


'U^ 


N 


X  100 


(7.8) 


Dashes  in  the  columns  of  Tables  7.2  and  7.3  indicate  that  a  measurement  was  not 
made  at  the  given  x  location. 

Surface  elevations 

Comparison  of  the  mean  calculated  infragravity  surface  elevations  to 
the  measured  elevations  illustrates  clearly  the  inaccuracies  in  calculating  the  mean 
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radiation  stress.  The  mean  setdown  and  setup  are  driven  by  the  gradients  in  the 
mean,  radiation  stress  (caused  by  wave  shoaling  and  breaking).  As  discussed 
previously,  the  wave  height  variation  around  breaking  is  not  well  captured  in  the 
measurements  due  to  the  fairly  coarse  spacing  of  the  wave  gauges  (3.66  m).  This 
results  in  a  similar  inaccuracy  in  the  estimation  of  the  radiation  stress  forcing. 
Linear  theory  used  to  calculate  radiation  stress  is  known  to  be  inaccurate. 
Furthermore,  the  simple  method  used  to  Calculated  radiation  stress  does  not 
account  for  the  transition  region  between  incipient  breaking  and  wave  setup 
(Svendsen  1984a)  where  changes  in  wave  shape  counteract  the  rapidly  decreasing 
wave  height  to  produce  essentially  constant  radiation  stress.  This  is  illustrated  by 
looking  at  the  most  shoreward  points  in  Tables  7.2  and  7.3.  For  Run  A0509A,  the 
calculations  show  significant  wave  setup  at  a:  =  47.59  m,  and  the  maximum 
setdown  occurring  at  x  =  29.26  m.  The  measurements  show  the  maximum 
setdown  is  much  further  shoreward  at  x  =  43.89  m  and  only  a  small  setup  at  x  = 
47.59  m.  The  over  estimation  of  setup  at  the  shoreline  in  the  numerical  model 
results  in  exaggerated  setdown  in  the  offshore  region  (in  the  closed  calculation 
domain,  overprediction  of  the  mean  wave  level  in  one  area  will  cause  an 
undeiprediction  in  other  areas).  The  mean  setdown  is  approximately  twice  the 
measured  values,  although  the  maximum  calculated  setdown  is  about  equal  to  the 
measured  value.  For  Run  A2007B  where  all  the  gauges  except  the  most  seaward 
is  in  the  surf  zone,  the  setup  calculated  at  the  gauge  locations  is  as  much  as  an 
order  of  magnitude  larger  than  the  measurements. 


Next,  looking  at  the  rms  error  in  the  infragravity  wave  calculations,  the 
average  rms  error  in  surface  elevation  is  5.8%  of  (0.040  m)  for  Run  A0509A 
and  10.3%  of  (0.044  m)  for  Run  A2007B.  The  error  increases  in  the 
shallowest  depths.  The  reason  for  the  increase  in  error  at  shallower  depths  is 
three-fold:  (a)  again,  the  linear  radiation  stress  calculations  are  less  accurate  in 
shallower  depths  as  the  waves  become  more  nonlinear,  (b)  the  measurements  are 
less  accurate  in  shallower  depths  as  the  sand  depth  at  the  gauge  changes  during  the 
run,  potentially  causing  a  change  in  the  gauge  offset  calculated  at  the  beginning  of 
the  run,  and  (c)  the  error  in  mean  water  level  increases  substantially  in  shallow 
water.  Up  to  55  percent  of  the  rms  error  accounts  for  the  mean  errors  previously 
discussed.  The  error  in  the  mean  setup  accounts  for  a  large  percentage  of  the  rms 
error  in  all  locations  with  a  positive  mean  water  level,  where  large  mean  errors 
occurred.  Removing  the  mean  setup  from  the  total  water  surface  elevations 
reduced  the  average  rms  errors  to  5.1%  (0.035  m)  and  7.3%  (0.031  m)  for  Runs 
A0509A  and  A2007B,  respectively. 

Velocities 

The  calculated  time  means  of  the  depth-integrated  cross-shore  velocities 
are  approximately  zero,  as  would  be  expected  from  a  depth-integrated  model  with 
no  net  inflow  or  outflow.  However,  the  depth-averaged  measured  mean  velocities 
do  show  a  net  offshore  flow  or  undertow.  This  is  due  to  the  fact  that  all 
measurements  used  were  below  the  still-water  level  and  the  wave  mass  flux,  which 
is  included  in  the  calculated  volume  flux,  occurs  above  that  level.  Data  from 
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gauges  shown  in  Figure  7.2  that  were  at  or  above  the  still-water  level  tended  to  be 
unreliable  and  were  neglected.  The  measurements  are  capturing,  to  varying 
degrees  depending  on  the  measurement  elevation,  the  return  flow  (undertow)  from 
the  wave-driven  mass  transport.  The  return  flow  is  generally  expressed  as 

=  C_  ^/gh  —  (7.9) 

mass  mass  *  ®  ^ 

where  is  the  return  flow  velocity  and  is  a  variable  coefficient  (see 
Svendsen  1984b).  The  values  of  C^,  based  on  the  depth-averaged  measured 
velocities  and  significant  wave  heights,  range  firom  0.02  to  0.13,  with  the  larger 
values  in  the  surf  zone.  The  model  cannot  represent  the  mean  measured  undertow 
without  explicitly  specifying  the  mass  transport,  which  is  not  done  here. 

The  average  rms  error  in  depth-averaged  velocity  is  3.7%  ofV gh 
(0.13  m/sec)  for  Run  A0509A  and  5.0%  ofVgT  (0.13  m/sec)  for  Run  A2007B. 
Removing  the  mean  from  the  velocities  reduced  the  average  rms  error  for  Run 
A0509A  slightly  (6  percent)  and  the  error  for  Run  A2007B  was  reduced  22  percent 
to  3.9%  (0.10  m/sec).  The  errors  were  generally  fairly  constant  throughout  the 


domain. 


7.3  Profile  Model 


7.3.1  Model  input  and  set  up 

The  profile  model  requires  the  water  surface  gradient,  radiation  stress, 
bottom  friction  coefficient,  eddy  viscosity,  grid-like  parameters,  and  initial 
conditions  as  input. 

Water  surface  gradient,  radiation  stress,  and  bottom  friction 

The  input  parameters  of  the  water  surface  gradient  and  the  radiation 
stress  are  closely  tied  to  the  depth-integrated  model.  As  described  in  Section 
4.2.1,  the  infragravity  water  surface  elevations  are  obtained  from  the  output  of  the 
depth-integrated  model.  The  radiation  stress  calculated  within  the  depth-integrated 
model  is  saved  along  with  the  water  surface  elevations.  Input  for  the  profile  model 
are  gradients  of  T  and  radiation  stress  which  are  calculated  with  second-order 
finite  differences  within  the  profile  model.  Note  that  the  Ax  for  the  derivatives  is 
the  same  as  used  in  the  depth-integrated  model,  although  the  velocity  profiles  can 
be  calculated  anywhere  within  the  grid  used  for  the  depth-integrated  model  without 
determining  the  solution  for  the  entire  domain.  The  bottom  friction  is  calculated 
using  Equation  2.37  with  a  constant  bottom  friction  coefficient  of^w  =  0.01,  as  in 
the  depth-integrated  model.  It  is  very  important  that  forcing  in  the  depth-integrated 
and  profile  models  be  consistent.  Otherwise,  the  water  surface  gradient  calculated 
in  the  depth-integrated  model  will  not  be  in  equilibrium  with  the  other  profile 
model  forcing. 


Eddy  viscosity 


‘  The  eddy  viscosity  is  calculated  using  Equation  4,49.  For  complicated 
cross-shore  bathymetry,  such  as  SUPERTANK  Run  A2007B,  the  "inside  surf 
zone"  formulation  for  eddy  viscosity  is  applied  for  all  points  shoreward  of  the  most 
offshore  point  where  —  0.8. 

Grid-like  parameters  and  initial  conditions 

The  profile  model  was  applied  at  all  jc  locations  in  the  depth-integrated 
grid,  using  the  central  point  in  y.  Profiles  were  determined  using  10 
eigenfunctions  and  specifying  output  at  21  locations  over  the  vertical.  The  initial 
condition  was  U  =  0  m/sec  for  all  x  and  z. 

7.3.2  Model  results 

Example  profile  model  results  are  shown  in  Figures  7.15-7.24.  In 
Figures  7.15-7,22,  the  solid  line  is  the  profile  model  ouq)ut,  the  dashed  line  is  the 
depth-integrated  model  output,  and  the  symbols  are  measurements.  The  vertical 
axis  is  the  dimensionless  depth  Z,  where 

z  =  hll  (7.10) 

K 

The  six  panels  in  each  figure  represent  six  snapshots  in  time  at  3-sec  intervals. 
Figures  7.15-7.18  are  results  from  Run  A0509A,  and  Figures  7.18-7.22  are  results 
from  Run  A2007B.  Figures  7.23  and  7.24  are  overlays  of  calculated  velocity 
profiles  at  3-sec  time  increments  from  t  =  500  to  750  sec,  for  Runs  A0509A  and 


Figure  7.15.  Velocity  profiles  at  x  =  25.6  m  (solid  --  profile  model,  dashed  -- 
depth-integrated  model,  symbols  -  measurements)  for  Run  A0509A 


U  (ny^)  U  (n/sec)  U  (rr/sec) 


Figure  7.16.  Velocity  profiles  atx  =  36.6  m  (solid  -  profile  model,  dashed  -- 
depth-integrated  model,  symbols  —  measurements)  for  Run  A0509A 
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Figure  7.17.  Velocity  profiles  atx  =  40.2  m  (solid  —  profile  model,  dashed  - 
depth-integrated  model,  symbols  —  measurements)  for  Run  A0509A 


Figure  7.18.  Velocity  profiles  at  x  =  43.9  m  (solid  ~  profile  model,  dashed  - 
depth-integrated  model,  symbols  —  measurements)  for  Run  A0509A 
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Figure  7.21.  Velocity  profiles  at  x  =  40.2  m  (solid  --  profile  model,  dashed  - 
dqith-integrated  model,  symbols  —  measurements)  for  Run  A2007B 


Figure  7.22.  Velocity  profiles  at  a:  =  43.9  m  (solid  --  profile  model,  dashed  - 
depth-integrated  model,  symbols  —  measurements)  for  Run  A2007B 
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A2007B,  respectively.  Error  statistics  and  mean  velocities  at  the  measurement 
locations  are  given  in  Tables  7.4  and  7.5  for  Runs  A0509A  and  A2007B, 
respectively.  The  tables  give  the  x  locations  (see  Figures  7. 1-7.2),  water  depths, 
and  dimensionless  vertical  positions  Z  (Equation  7.10)  at  which  the  measurements 
were  made.  The  dimensionless  mean  measured  and  calculated  velocities  are 
defined 


where  U„  is  the  measured  velocity,  is  the  velocity  calculated  using  the  profile 
model,  and  N  is  the  number  of  measurements  in  the  time  series.  The 
dimensionless  rms  error  in  velocity  is  defined 
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Table  7.4.  Profile  model  mean  values  and  rms  errors  for  Run  A0509A 


x/m 

h,  m 

Z 

% 

% 

^U,rmsy  ^ 

Ev.^,  % 

25.60 

1.67 

0.14 

-0.82 

-0.71 

2.24 

2.24 

25.60 

1.67 

0.41 

-1.18 

-0.39 

2.28 

2.14 

25.60 

1.67 

0.69 

-1.60 

0.49 

3.28 

2.54 

25.60 

1.67 

0.87 

1.77 

1.15 

3.20 

3.14 

34.75 

1.38 

0.20 

-1.85 

-1.08 

3.56 

3.50 

34.75 

1.38 

0.41 

-1.85 

-0.62 

3.08 

2.83 

34.75 

1.38 

0.61 

-1.66 

0.19 

3.60 

3.09 

34.75 

1.38 

0.84 

-0.05 

1.73 

4.56 

4.20 

36.58 

1.33 

0.14 

-1.75 

-0.51 

3.28 

3.03 

36.58 

1.33 

0.37 

-1.56 

-0.32 

3.57 

3.35 

36.58 

1.33 

0.59 

-1.58 

0.16 

4.01 

3.62 

36.58 

1.33 

0.82 

-0.38 

0.75 

3.37 

3.17 

40.23 

1.13 

0.12 

-1.71 

-0.56 

3.98 

3.81 

40.23 

1.13 

0.53 

-2.60 

-0.003 

5.09 

4.38 

40.23 

1.13 

0.81 

0.65 

0.70 

3.80 

3.80 

43.89 

0.88 

0.20 

-4.03 

-1.62 

6.87 

6.44 

43.89 

0.88 

0.73 

-0.69 

1.41 

5.68 

5.28 

*rms  error  with  mean  removed 
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Table  7.5.  Profile  model  mean  values  and  rms  errors  for  Run  A2007B 

h,  m 

Z 

u^. 

Van, 

U^,  m/sec 

UJ,  m/sec 

36.58 

0.74 

0.35 

-6.02 

-1.55 

lAl 

5.62 

36.58 

0.74 

0.87 

-4.07 

-2.82 

9.87 

7.07 

37.19 

0.88 

0.23 

-4.68 

0.61 

IM 

5.52 

37.19 

0.88 

0.56 

-3.19 

-0.07 

5.98 

5.11 

37.19 

0.88 

0.88 

-0.58 

-1.53 

5.31 

5.22 

40.23 

0.50 

0.20 

00 

od 

t 

-1.23 

9.90 

7.07 

40.23 

0.50 

0.82 

0.71 

1.18 

8.85 

8.84 

43.89 

0.59 

0.31 

-5.79 

-1.59 

6.97 

5.57 

43.89 

0.59 

0.84 

2.95 

2.64 

6.58 

6.58 

*rms  error  with  mean  removed 


Run  A0509A 

Profiles  for  Run  A0509A  are  given  in  Figures  7.15  -  7.18  and  range 
from  locations  well  outside  the  surf  zone  at  ;r  =  25.6  m  to  the  outer  surf  zone  at  x 
=  43.9  m.  These  figures  show  a  short  sequence  of  profiles  at  a  time  when  the 
depth-integrated  model  performed  well  in  terms  of  predicting  the  infragravity 


motion  (see  Figures  7.7  and  7.8).  Figure  7.23  gives  an  overview  of  profile  shapes 
outside  and  inside  the  surf  zone. 

Outside  the  surf  zone  (x  <  40  m),  typical  velocity  profiles  exhibit  three 
characteristics:  (a)  the  greatest  velocity  variation  generally  occurs  near  the 
surface,  (b)  the  velocities  near  the  surface  and  bottom  are  out  of  phase  (note  the 
"braided"  appearance  of  the  profiles  near  the  surface),  and  (c)  the  profile  variation 
over  depth  is  relatively  symmetric  about  U  =  0  m/sec  (time  mean  of  the  velocity 
at  any  elevation  is  near  zero).  These  characteristics  are  illustrated  in  Figures 
7.15  -  7.17  and  7.23.  The  velocity  profiles  are  uniform  over  the  lower  two-thirds 
of  the  depth.  Near  the  surface,  the  velocity  magnitude  is  greater  and  the  direction 
is  often  opposite  that  of  the  lower  part  of  the  water  column.  Typically,  in  the 
offshore  region  the  time-varying  surface  stress  is  the  dominant  forcing.  Since  the 
mean  forcing  is  small,  the  profiles  are  approximately  symmetric  about  U  = 

0  m/sec.  The  water  depth  is  relatively  large,  so  the  surface  stress  causes  higher 
velocities  near  the  surface. 

Inside  the  surf  zone  (x  >  40  m),  velocity  profiles  have  the  following 
characteristics:  (a)  the  mean  undertow  profile  is  the  dominant  featurej  and  (b)  the 
time-variation  is  fairly  constant  in  magnitude  and  phase  over  depth  (the  profiles  are 
approximately  parallel,  instead  of  the  braided  pattern  outside  the  surf  zone).  The 
mean  velocity  profile  inside  the  surf  zone  is  the  classic  parabolic  shape  (see 
Figures  7.18  and  7.23),  unlike  outside  the  surf  zone  where  the  mean  profile  is 
approximately  vertical.  This  implies  that  the  dominant  forcing  in  the  outer  surf 


zone  is  the  mean  surface  stress  (mean  surface  stress  portion  of  the  radiation  stress 
gradient  due  to  wave  breaking).  The  typical  time  variation  of  the  profile  in  the 
surf  zone  is  an  oscillation  about  the  mean  undertow  profile,  although  there  is  a 
small  phase  difference  between  the  surface  the  bottom.  Figures  7.18  and  7.23 
show  that  the  velocity  profiles  are  relatively  vertical  for  positive  (onshore) 
velocities  and  more  parabolic  for  negative  (offshore)  velocities. 

The  average  rms  error  in  the  calculated  velocities  is  3.9%  oiV gh 
(0.14  m/sec),  which  is  approximately  equal  to  the  error  in  the  depth-integrated 
model.  Considering  only  the  time-varying  solution,  the  rms  error  decreases 
slightly.  Mean  velocities  are  in  general  not  well  represented  in  the  model.  The 
reasons  for  the  disagreement  are  the  same  as  discussed  in  Section  7.2.2  for  the 
depth-integrated  model.  The  error  is  uncorrelated  with  elevation  in  the  water 
column,  but  error  does  increase  in  the  surf  zone.  The  profile  model  reproduces 
the  characteristics  and  the  profile  shapes  measured  in  the  SUPERTANK  runs. 

Run  A2007B 

Profiles  for  Run  A2007B  are  given  in  Figures  7.19-7.22  and  7.24,  and 
range  from  locations  just  outside  the  surf  zone  through  the  outer  and  mid  surf 
zone.  The  profile  characteristics  discussed  in  the  surf  zone  for  Run  A0509A  apply 
here  as  well,  i.e.,  dominance  of  parabolic  profile  shape,  and  time-varying  phase 
and  magnitude  are  constant  over  depth.  In  addition,  in  the  mid  to  inner  surf  zone 
(j:  >  40  m),  model  results  and  measurements  show  less  velocity  variation  over  the 
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vertical  as  the  time-varying,  depth-uniform  forcing  (surface  elevation  gradient) 
becomes  dominant.  This  is  shown  clearly  in  Figure  7.24  at  ;c  =  50.0  m. 

The  average  rms  error  in  the  calculated  velocities  for  Run  A2007B  is 
7.6%  ofV~gh  (0.20  m/sec),  which  is  about  50  percent  greater  that  the  depth- 
integrated  model  results.  Considering  only  the  time-varying  solution,  the  rms  error 
decreased  to  6.3%  (0.16  m/sec).  Both  Runs  A0509A  and  A2007B  show  an 
increase  in  calculated  velocity  errors  in  and  near  the  surf  zone. 

The  main  source  of  error  in  the  calculations  is  inaccuracies  in  the 
forcing  (linear  radiation  stress  calculation  and  omission  of  the  transition  zone  in  the 
region  of  incipient  wave  breaking,  see  Section  7.2.2).  Additional  errors  stem  from 
the  linear  separation  of  the  incident  and  reflected  long  wave  used  to  drive  the 
dq)th-integrated  model,  the  linear  separation  of  the  short-  and  long-wave 
components,  and  the  lack  of  resolution  in  the  water  surface  measurements  for 
calculating  radiation  stress.  Laboratory  measurement  errors  are  also  a  factor. 

7.3.3  Profile  parameters 

As  discussed  in  Section  5.2.2,  Putrevu  and  Svendsen  (1995)  identified 
dimensionless  parameters  that  together  describe  infragravity  wave  profiles: 
which  determines  the  depth  to  which  the  effect  of  the  surface  and  bed  shear 
stresses  are  felt  (|3  =  (<ji/(2  w  is  the  infragravity  wave  frequency,  and  h  is 
water  depth);  A  which  is  the  relative  phase  between  the  surface  stress  and  depth- 
uniform  forcing;  and  /3/(p  fj  which  is  the  relative  importance  of  the  surface 


stress  and  depth-uniform  forcing,  where  T,  and^  are  the  amplitudes  of  the  surface 
stress  and  depth-uniform  forcing,  respectively,  A  summary  of  estimated  parameter 
values  for  the  SUPERTANK  Runs  A0509A  and  A2007B  are  shown  in  Figures 
7.25  and  7.26.  The  parameters  are  estimated  using  peak  infragravity  wave 
frequencies  of  0,039  Hz  for  Run  A0509A  and  0.027  Hz  for  Run  A2007B.  The 
peak  frequencies  for  the  surface  stress  and  depth-uniform  forcing  were  not  exactly 
the  same.  Also,  the  peak  frequencies  varied  somewhat  across  the  grid.  The 
interpretation  of  parameters  values  illustrates  the  mechanisms  behind  the  trends 
discussed  previously,  including  penetration  of  surface  stresses,  the  relative  phasing 
of  the  forcing,  and  the  dominant  forcing. 

Parameter  Bh 

The  larger  value  of  j8/i  outside  the  surf  zone  implies  that  surface  stress 
forcing  will  impact  a  smaller  percentage  of  the  water  column  than  at  locations 
inside  the  surf  zone,  and  thus,  the  greatest  vertical  variation  will  be  near  the 
surface.  For  Run  A0509A,  jSA  varies  from  10  to  1,  which  is  reflected  in  the 
difference  between  the  profiles  with  greater  variations  near  the  surface  outside  the 
surf  zone  (large  /Sft)  and  profiles  with  a  uniform  variation  in  magnitude  over  depth 
in  the  surf  zone  (small  jSft).  For  Run  A2007B,  is  relatively  constant  (1.5  to 
1.0),  and  the  profile  variation  over  depth  is  similar  to  the  surf  zone  profiles  of  Run 
A0509A  which  have  similar  values  of  ^h. 
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Figure  7.26.  Relative  strength  of  surface  stress  and  depth-uniform  forcing  (solid 
line),  (dashed  line),  and  relative  phase  of  surface  stress  and  depth-uniform 
forcing  for  Run  A2007B 
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Parameter  A 


;  The  relative  phase  between  the  surface  stress  and  depth-uniform  forcing 
is  the  parameter  A.  This  parameter  gives  the  simplest  interpretation  of  the  results. 
Physically,  however,  this  parameter  is  difficult  to  predict  because  the  depth- 
uniform  component  of  the  forcing  is  composed  of  both  the  pressure  gradient  from 
the  long-wave  motion  and  the  depth-uniform  part  of  the  variation  of  the  radiation 
stress  due  to  the  modulation  of  the  short-wave  motion.  For  the  cases  of  pure  surf 
beat  analyzed  in  Chapter  5,  the  value  of  A  was  approximately  constant  and  equal  to 
X  outside  the  surf  zone  and  0.5-0. 8  x  inside  the  surf  zone.  The  trend  was 
dependent  on  the  choice  of  x  =  0. 

The  SUPBRTANK  measurements  show  similar  trends  in  A.  For  both 
Runs  A0509A  and  A2007B,  A  «  x  outside  the  surf  zone  (same  as  for  the  surf  beat 
model)  and  in  the  outer  surf  zone.  Inside  the  surf  zone  (x  >  40  m  for  Run 
A0509A  and  the  inner  surf  zone,  x  >  45  m,  for  Run  A2007B)  A  decreases  to  x/2. 
Thus,  inside  the  surf  zone,  the  surface  stress  and  depth-uniform  forcing  are  closer 
to  being  in  phase.  The  decease  in  A  in  the  surf  zone  indicates  that  the  upper  and 
lower  portions  of  the  velocity  profile  should  also  be  less  out  of  phase,  as  seen 
previously. 

Parameter  B/(o  f„) 

The  last  parameter,  /3/(p 7^,  is  a  measure  of  the  ratio  of  the  surface 
stress  (part  of  the  radiation  stress  due  to  short-wave  modulation)  to  the  depth- 
uniform  forcing  (contribution  of  the  pressure  gradient  due  to  the  infragravity 


motion  and  the  modulation  of  the  short-wave  radiation  stress  below  trough  level, 
which  is  assumed  to  be  depth  uniform).  The  parameter  shows  that  the  forcing 
outside  the  surf  zone  tends  to  be  surface  stress  dominated  ^/(p  fj  >  0.7), 
while  in  the  surf  zone  the  forcing  is  balanced  between  surface  stress  and  depth- 
uniform  forcing.  In  Run  A2007B  (Figure  7.26)  the  depth-uniform  forcing 
dominates  (T^  j8/fp  fJ  <  0.7)  in  the  inner  surf  zone.  In  the  relatively  long,  flat 
inner  surf  zone  in  Run  A2007B  the  short-wave  energy  is  dissipated  (T^  becomes 
small)  and  the  profile  shape  is  vertical  (depth-uniform  forcing  dominates,  see 
Figure  7.24,  x  =  50  m).  Surface  stress  dominated  profiles  vary  over  depth,  and 
depth-uniform  dominated  profiles  are  linear  and  vertical  in  shape. 

The  model  results  for  SUPERTANK  Runs  A0509A  and  A2007B  exhibit 
the  following  characteristics.  Outside  the  surf  zone,  the  greatest  velocity  variation 
occurred  near  the  surface  (large  and  large  /3/(p  fJ).  Also,  velocities  near  the 
surface  and  the  bottom  were  out  of  phase  (A  =  tt).  The  profile  variation  over 
dqith  was  fairly  symmetric  about  U  =  0  m/sec  because  the  time-mean  forcing  was 
small.  Inside  the  surf  zone,  the  mean  undertow  profile  was  the  dominant  feature 
due  to  the  dominance  of  the  time-mean  forcing.  The  time-variation  of  the  velocity 
was  fairly  constant  over  depth  (small  jS/i  and  small  A).  Small  values  of  A  indicate 
that  the  surface  and  depth-uniform  forcing  are  in  phase,  so  the  velocities  at  the 
surface  and  through  the  water  column  are  also  in  phase.  In  the  inner  surf  zone  or 
a  flat  surf  zone,  the  velocity  profile  was  vertical  due  to  the  small  value  of 
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Chapter  8 

SUMMARY  AND  CONCLUSIONS 


The  purpose  of  this  study  was  to  numerically  model  time-  and  depth- 
varying  cross-shore  currents  forced  by  gravity  and  inftagravity  waves  in  the 
nearshore.  The  study  included  numerical  model  development,  laboratory  data 
collection,  and  model  application  and  evaluation  with  the  laboratory  data. 

8.1  Research  Results 

Both  depth-integrated  and  depth-varying  models  were  developed  and 
linked.  The  depth-integrated  model  includes  the  propagation,  generation,  and 
reflection  of  long  waves  in  a  nearshore  domain.  The  ou^ut  water  surface 
elevations  from  the  depth-integrated  model  were  used  to  drive  a  local  eigenfunction 
solution  of  the  velocity  over  dqpth.  Both  models  are  time  dependent. 

Using  the  depth-integrated  equations  of  mass  and  momentum,  two 
depth-integrated,  finite-difference  flow  models  were  developed  based  on  the 
Predictor-Corrector  (two-dimensional)  and  Lax-Wendroff  (one-dimensional) 
methods.  The  Predictor-Corrector  model  was  a  modification  of  the  SHORECIRC 
model.  Model  forcing  terms  included  radiation  stress,  bottom  friction,  and  a 
specified  incoming  long  wave.  The  Predictor-Corrector  model  exhibited  spurious 


water  elevation  and  velocity  oscillations  in  time  when  it  was  executed  without 
bottom  Motion  and  depth  nonlinearity.  This  problem  is  similar  to  spatial 
oscillations  in  a  Leap-Frog  model. 

The  velocity  profile  model  is  an  analytical  eigenfunction  solution  for 
the  cross-shore  momentum  equation,  assuming  longshore  uniformity.  Eddy 
viscosity  was  assumed  constant  for  a  given  depth,  but  varied  across  shore  as  a 
function  of  local  water  depth  and  breaker  location.  Both  the  depth-integrated  and 
profile  models  shared  the  same  radiation  stress  and  bottom  Motion  formulations,  to 
ensure  equilibrium  between  the  solutions. 

The  depth-integrated  and  profile  models  were  applied  to  the  analytical, 
depth-integrated  surf  beat  solution  of  Schaffer  and  Svendsen  (1988)  and  Schaffer 
(1990,  1993).  The  Schaffer-Svendsen  solution  includes  a  forced  infragravity  wave 
propagating  across  a  plane  beach.  The  infragravity  wave  is  released  through 
breaking  and  reflects  from  the  shoreline.  The  Schaffer-Svendsen  solution  was 
modified  to  include  both  mean  and  time-varying  forcing  (original  model  included 
only  the  time-varying  component).  Three  test  cases  simulated  one  condition 
published  by  Schaffer  (1993)  and  two  conditions  similar  to  the  SUPERTANK 
configuration.  The  Lax-Wendroff  depth-integrated  model  was  applied  since  the 
analytical  surf  beat  solution  was  linear  and  neglected  bottom  Motion.  The  depth- 
integrated  model  results  compared  well  to  the  analytical  solution.  Root-mean- 
square  errors  in  ranged  from  0.25  to  1.0%  of  Hg  (0.(X)2  to  0.01  m)  and  rms 
errors  in  V ranged  from  0.02  to  0.9%  of Vgh  (0.002  to  0.04  m/sec).  The 


greatest  errors  were  typically  near  antinodes  in  the  solution  at  the  shoreline.  The 
solutions  showed  antinodes  in  elevation  at  the  shoreline  and  near  the  breaker  line. 

The  profile  model  was  also  run  with  the  Schaffer-Svendsen  surf  beat 
forcing.  The  velocity  profiles  outside  the  surf  zone  leaned  offshore  at  the  surface 
and  were  fairly  uniform  over  the  lower  two-thirds  of  the  depth.  This  was  due  to 
the  shoaling  of  the  short  waves  and  the  relatively  large  water  depth.  In  the  mid 
surf  zone,  the  velocity  profiles  exhibited  a  parabolic  undertow  shape  with  offshore 
flow  near  the  bottom  and  onshore  flow  near  the  surface.  This  shape  was  due  to 
the  gradient  in  wave  height  caused  by  wave  breaking.  In  the  very  nearshore,  the 
velocity  profiles  varied  uniformly  over  depth  because  the  depth-uniform  forcing 
(water  surface  gradients)  dominated. 

Putrevu  and  Svendsen  (1995)  developed  an  infragravity  wave  velocity 
profile  model,  assuming  a  constant  periodicity  in  all  forcing  functions.  This 
condition  is  satisfied  in  the  Schaffer-Svendsen  surf  beat  model.  Using  only  the 
time-varying  portion  of  tiie  forcing  (neglecting  the  mean  forcing),  the  eigenfunction 
profile  model  results  were  compared  to  the  Putrevu-Svendsen  solution.  The 
profiles  were  the  same,  as  would  be  expected. 

Putrevu  and  Svendsen  (1995)  identify  three  dimensionless  parameters 
that  represent  the  relative  importance  of  the  surface  stress  forcing,  depth-uniform 
forcing,  and  the  bottom  stress.  For  the  surf  beat  examples,  the  parameters 
elucidated  the  following  trends:  (a)  outside  the  surf  zone,  the  surface  stress  term 
dominated  and  profiles  showed  the  greatest  variation  near  the  surface, 


(b)  progressing  through  the  surf  zone,  the  depth-uniform  forcing  became  more 
important  and  the  profile  became  uniform  over  depth,  and  (c)  bottom  stress  was 
fairly  constant  through  the  nearshore  and  was  less  important  than  the  surface  stress 
and  depth-uniform  forcing  terms. 

Two  additional  factors  introduced  by  Putrevu  and  Svendsen,  and  A, 
indicate  the  depth  to  which  surface  stress  forcing  penetrates  and  the  phase 
difference  between  the  depth  uniform  and  surface  stress  terms,  respectively.  The 
parameter  decreases  from  offshore  to  nearshore  indicating  greater  penetration  of 
the  surface  stress  forcing  near  the  shoreline.  This  reinforces  the  trend  toward 
depth  uniform  profiles  in  shallow  depths  and  greater  vertical  nonuniformity  in 
deeper  depths.  For  the  surf  beat  results,  A  was  approximately  equal  to  x.  In  the 
inner  surf  zone,  the  value  of  A  decreased  to  0.5  -  0.8  x.  Since  in  the  inner  surf 
zone  the  short-wave  energy  is  dissipated  (surface  stress  forcing  near  zero),  the 
change  in  A  did  not  appear  to  impact  the  profile  shape. 

Water  surface  elevation  and  cross-shore  current  data  were  collected 
during  the  large-scale  SUPERTANK  Laboratory  Project  to  provide  measurements 
to  evaluate  the  depth-integrated  and  profile  models.  SUPERTANK  was  performed 
in  the  large  wave  tank  at  Oregon  State  University  with  a  sand  beach.  A  range  of 
near-prototype  wave  conditions  were  run,  and  waves,  water  levels,  currents,  and 
beach  profile  evolution  were  measured.  SUPERTANK  hydrodynamic  data  were 
filtered  with  a  cut  off  of  half  the  peak  frequency  to  separate  short-  and  long-wave 
componraits.  Details  of  the  analysis  and  sample  data  were  presented  in  Chapter  6. 


Two  STJPERTANK  runs  were  identified  to  simulate  with  the  models, 
Runs  A0509A  and  A2007B.  The  runs  differed  from  each  other  in  peak  period  (3 
versus  5  sec),  bathymetry  (non-barred  versus  barred),  spectral  peakedness  (7  =  20 
versus  100),  and  gauge  positions  (outside  the  surf  zone  and  outer  surf  zone  versus 
outer  to  mid  surf  zone). 

The  depth-integrated  Predictor-Corrector  model  grid  was  set  up  so  that 
gauge  locations  aligned  with  grid  cells.  Still-water  depths  were  specified  using  the 
pre-run  profile  surveys.  The  depth-integrated  model  was  forced  with  measured 
incident  long-wave  time  series  at  the  offshore  boundary  and  radiation  stress 
throughout  the  domain  calculated  from  measured  or  interpolated  short-wave  time 
series.  The  radiation  stress  was  calculated  using  linear  theory.  A  constant  bottom 
friction  coefficient  off^  =  0.01  was  applied. 

The  average  rms  errors  for  the  depth-integrated  model  were  1  to  21% 
of  (0.04  m)  for  Y  and  2  to  6%  ofVYh  (0.13  m/sec)  for  U.  The  rms  errors 
in  surface  elevation  increased  in  shallow  depths.  This  increase  in  shallow  water 
was  due  to  the  linear  calculation  of  radiation  stress  (due  to  increased  wave 
linearity)  and  an  increase  in  mean  gauge  errors  at  shallow  depths  (due  to  sand  bed 
changes).  Up  to  22  percent  of  the  rms  error  in  U  is  due  to  mean  errors.  The 
measured  depth-averaged  velocities  are  weighted  toward  the  lower  water  column 
and  thus  tended  to  show  a  net  offshore  flow,  instead  of  the  net  zero  flow  expected 
for  a  closed  system. 


The  profile  model  was  executed  with  10  eigenfunctions  for  21  locations 
over  depth.  The  water  elevations  and  radiation  stress  were  supplied  by  the  depth- 
integrated  model.  The  same  bottom  friction  coefficient  was  used  as  used  in  the 
dq)th-integrated  model.  The  rms  error  in  the  calculated  velocities  was  4%  of 
V^h  (0.14  m/sec)  for  Run  A0509A  (most  gauges  outside  the  surf  zone)  and  8%  of 
V^h  (0.20  m/sec)  for  Run  A2007B  (most  gauges  inside  the  surf  zone).  Errors 
increased  in  and  near  the  surf  zone  for  both  runs.  Mean  errors  again  accounted  for 
as  much  as  20  percent  of  the  rms  error.  Velocity  errors  were  larger  for  the  profile 
model  than  the  depth-integrated  model.  Mean  velocities  are  generally  not 
rq)resented  well  by  the  model.  Errors  were  uncorrelated  with  elevation  in  the 
water  column.  The  profile  model  results  were  quite  good  during  times  when  the 
d^th-integrated  results  were  good 

The  overall  trends  in  the  profiles  were  similar  to  those  seen  in  the  surf 
beat  simulations.  The  profile  model  results  for  SUPERTANK  exhibited  the 
following  characteristics: 
a)  outside  the  surf  zone 

♦  the  greatest  velocity  variation  occurs  near  the  surface  (large 
and  surface  stress  dominating) 

♦  velocities  near  the  surface  and  the  bottom  are  out  of  phase 
(A  =  t) 

♦  profile  variation  over  depth  is  fairly  symmetric  about  U  = 

0  m/sec  (small  mean  forcing) 


b)  inside  the  surf  zone 

i  ♦  mean  undertow  profile  is  the  dominant  feature  (mean  forcing  is 
dominating) 

♦  the  time  variation  of  the  velocity  is  fairly  constant  over  depth 
(small  and  small  A) 

♦  in  the  inner  surf  zone  or  a  flat  surf  zone,  the  velocity  profile  is 
uniform  over  depth  (small  surface  stress  forcing) 

8.2  Future  Work 

The  greatest  limitation  of  the  present  study  was  the  specification  of 
radiation  stress  using  linear  theory  and  coarsely  spaced  water  surface 
measurements.  Svendsen  (1984a)  and  Svendsen  and  Putrevu  (1993)  show  the 
limitations  of  linear  theory  in  the  surf  zone,  but  the  ease  of  application  for  a  broad 
range  of  conditions  has  kept  linear  theory  as  the  tool  of  choice.  The  coarsely 
spaced  wave  gauges  at  SUPERTANK  did  not  capture  the  breaker  position  for 
model  application.  Future  experiments  need  closely  spaced  gauges  in  the  surf  zone 
to  catch  the  important  transition  processes  in  the  breaker  zone.  A  model  of  the 
time-varying  nearshore  wave  transformation  and  breaking  would  also  fill  this  void. 
Until  such  a  model  is  available,  use  of  the  depth-integrated  and  profile  model  will 
be  limited  to  cases  with  dense  wave  measurements  or  use  of  ad-hoc  assumptions 
about  the  relationships  between  short  waves  and  infragravity  waves  (e.g.,  Symonds 
et  al.  1982,  Schaffer  and  Svendsen  1988). 


Another  improvement  to  the  depth-integrated  model  would  be  to 
formulate  a  wetting  and  drying  boundary  condition  at  the  shoreline  to  replace  the 
present  fixed  shoreline  location.  Prediction  of  swash  processes  is  useful,  but  not 
necessarily  important  to  study  velocity  profiles  away  from  the  shoreline. 

The  models  developed  in  this  study  focused  on  cross-shore  processes. 
The  next  step  would  be  to  generalize  the  solution  to  include  both  the  longshore  and 
cross-shore.  The  complications  in  doing  this  include  more  complicated  boimdary 
conditions,  refraction  processes,  and  solution  methods.  Application  of  such  a 
quasi-three-dimensional  model  could  be  used  to  represent  cell  circulation  and  rip 
current  generation,  which  are  still  poorly  understood.  Edge-wave  and  shear-wave 
forcing  could  also  be  applied. 

For  coastal  engineers,  calculation  of  hydrodynamic  processes  is  usually 
an  intermediate  step  in  the  calculation  of  sediment  transport  and  bathymetry 
change.  The  depth-integrated  and  profile  models  are  useful  tools  for  studying 
cross-shore  transport  of  suspended  sediment.  Together  with  the  sediment 
concentration  and  beach  profile  data  from  SUPERTANK,  the  models  could  be 
applied  to  study  sediment  transport  rates  and  quantities  under  controlled  conditions. 
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